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Frequency Dependence of D-Region Scattering at VHF 
J. C. Blair, R. M. Davis, Jr., and R. C. Kirby 


' 
| Contribution from Central Radio Propagation Laboratory, National Bureau of Standards, Boulder, Colo. 
(Received March 30, 1961) 


Results are given of a one-year program of observation of frequency dependence of 
system loss for D-region scattering at VHF. Continuous simultaneous observations were 
made at five frequencies from 30 to 108 Me/s, using narrow beam antennas scaled in dimen- 
sions and height according to wavelength; the path was from Long Branch, IIl., to 
Boulder, Colo., a distance of 1,300 km. Hourly values of system loss closely followed the 
relationship p,/paxf". The value of n varied with time of day and season, with 90 percent 
of all observations lying in the range 7 to 9% for scaled antenna systems. Diurnal and 
seasonal variation of n are attributed to changing relative roles of turbulence and meteoric 

















reflections. Effects of normal ionospheric absorption are not apparent. 


A special analysis 


is made of frequency dependence during weak signal conditions. Results are also given of 
studies of the frequency dependence of fading characteristics, the effect of beamwidth, and 


SID behavior. 
1. Introduction 


Continuous radio-wave scattering at very high 
frequencies, from turbulent irregularities and mete- 
oric ionization in the D-region of the ionosphere, has 
been intensively investigated [Bailey, Bateman and 
Kirby, 1955] since 1951, and has been exploited for 
communication systems in several areas of the world 
[JTAC, 1960]. This paper reports an experimental 
study of the frequency dependence in the range 30 
to 108 Me/s. One of the authors has already pub- 
lished preliminary data [Blair, 1959]. 

Frequency dependence is of considerable engi- 
neering importance in the design of ionospheric- 
scatter communication systems. The higher the 
frequency one can use, the greater the freedom from 
interference, especially that propagated to great dis- 
tances by the F2-layer and sporadic-F ionization. 
To raise the frequency, however, increases propa- 
gation losses until, in spite of the accompanying de- 
crease of the background level of galactic noise, 
excessive transmitter power is required for a needed 
transmission capacity. For design purposes one is 
especially interested in the frequency dependence 
characteristic of weak-signal conditions. 

Interpretation of the physics of the scattering 
process also requires knowledge of the frequency 
dependence. In studying the role of turbulence, 
under the usual assumptions of single isotropic 
scattering at an average angle 6, observations of fre- 
quency dependence serve to evaluate the spectrum of 
sizes of irregularities, S(k). The scattering process 
acts like a narrow band filter, emphasizing the sizes 
of irregularities which have wavenumbers equal to 


k= 2 sin (1) 
where 
\=wavelength 
6=scattering angle, shown in figure 1. 
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Wheelon [1960a] has recently summarized turbu- 
lence theory for D-region scatter. He treats the 
single scattering description of electromagnetic re- 
sponse of electron density irregularities in time and 
space, and discusses several theories for turbulent 
mixing of electron density. 

Neglecting absorption, the total received power in 
VHF scatter propagation is the resultant of compo- 
nents due to turbulent scattering, meteor reflections, 
and noise, which may be expressed as a function of 
frequency: 


Palf) =D), (A) +Pm (f)+ pw (f) (2) 


where 

Pa=available power at terminal of receiving antenna 

ps=component of available power from scattering 
from turbulent irregularities 

Pm=component of available power from meteor 
reflections 

py=component of available power from external 
noise usually mainly galactic 

f=carrier frequency, megavycles per second 


] 
a yan 
Se SE x<\5 
QPEL ES <> 





“7 Y “Me“M#_ g 
A? V 47 S/, # Vy Y/) ; 


Ficure 1. Geometry of D-region scattering. 








Other modes which sometimes occur, such as £,, 
auroral echoes, F2-propagated ground backscatter, 
are recognizable by the abnormal character of the 
recording, and such data are not used in the analysis 
of frequency dependence. Adequate power is used 
in the experiment to keep the signal-to-noise ratio 
large, usually greater than 10, so that py can be 
neglected. The effect of absorption can be indicated: 


pelf)=p.(f) exp Fas ‘ ” atid ai 


+pn(f) exp Fas. ie a sec od | (3) 


where 


Ja is the gyrofrequency, megacycles per second 

a is absorption coefficient, per element thick- 
ness, dh, of absorbing region 

h, is the effective scattering height 

hm is the height of meteoric reflection 


for the oblique transmission through the absorbing 
region enroute to and from the scattering region or 
meteor reflection. 

This is the first published report of comprehensive 
results of simultaneous observations at a series of 
frequencies, rather than a pair at a time, throughout 
the range of interest. The frequency dependence is 
expressed in terms of n in the relation 


P:|Pa Xf” 
where 


P:=power to transmitting antenna’ 
Pa=available power at receiving antenna 
f=frequenc) 


Earlier experiments showed a variation of n. The 
Bailey, Bateman and Kirby [1955] measurements, 
obtained in two nonsimultaneous tests, indicated 
lower values of n in the frequency range below 50 
Mc/s than in the range above 50 Mc/s. This led 
Wheelon [1957] to suggest a theoretical explanation 
for a dichotomy in the scaling law; he emphasized 
at the same time the need of carefully scaled simul- 
taneous measurements over the whole frequency 
range. 

Some interpretation is given in this report of the 
variation in relative contribution of meteoric reflec- 
tions and turbulent scattering, and the role of 
absorption. 


2. Experimental Conditions 


System loss, p,/~a, has been compared at a series 
of frequencies, using a fixed path length. The angle 
6 is nearly constant; diurnal and seasonal changes 
in height of scattering introduce slight variations. 

Identical antenna patterns and gains are preserved 
at all frequencies by scaling antenna dimensions and 
height. Thus, the scattering volume is illuminated 





and viewed uniformly throughout the frequency 
range. As a consequence of scaling antennas, an 
extra power of 2 is contained in the exponent, n, for 
system loss, because the aperture of the receiving 
antenna varies as \”. One must subtract a power 
of 2 from the values of n (divide system loss by f*) 
for scaled antennas given in this paper,’ to make 
comparison with theoretical frequency dependence of 
scattering cross section o, or with S(k) values. 

Figure 2 shows the geographical layout of the 
experimental path. Continuous transmissions were 
maintained at five frequencies from near 30 Mc/s to 
near 108 Me/s. Narrow-beam rhombic antennas 
were used, scaled in dimensions and height; the leg 
length was 25 wavelengths, one-half the obtuse angle 
was 80.2° and the horizontal half-power beamwidth 
was 6°. The elevation of the axes of the main lobes 
was 4.6°; the beams were directed to intersect at 
the path midpoint at an ionospheric height of 85 
km. Losses in transmitting and receiving transmis- 
sion line systems were measured, and observed data 
at each frequency corrected accordingly. Finally, 
aircraft measurements of antenna gains and radia- 
tion patterns were carried out. The average value 
of gain was 18.8 db relative to a halfwave reference 
dipole at the same height; measured values of gain 
for all antennas were found to be within +0.8 db of 
the average value. The sites were nearly ideally flat 
over at least the first Fresnel zone for ground reflec- 
tion, and the aircraft measurements of radiation 
patterns confirmed the uniformity of the main lobes 
at all frequencies. 

Transmitter power supplied to the antennas ranged 
from a nominal value of 2 kw at 30 Me/s to 20 kw 
at 108 Me/s; the power at each frequency was cali- 
brated and records were kept so that the received 
signal intensity could be adjusted to a standard 
reference level. Receiver noise bandwidths were 300 
c/s at all frequencies except 108 Me/s, for which a 
25 c/s band was used. 

! This use of n is at variance with its use in the references, including previous 
publications of the authors; it has been adopted because of the convenience of 
the system loss concept. Furthermore, one can effectively measure or interpret 
data only for the scaled antenna situation. Even in communication practice 


the concept of scaled (or constant gain) antennas is more convenient than that 
of constant-aperture antennas. 








Ficure 2. Experimental transmission path and frequencies. 
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In the long-term recording program a detector- 
AGC circuit with an integration time constant of 
12 sec was used; the recorder response was propor- 
tional to logarithm of the average signal voltage. 
For the short-term signal distributions encountered, 
comparison with fast-response recorder and totalizer 
measurements showed that median values given by 
the long-time constant record, scaled for periods of 
5 min or so, corresponded with true median values 
of the distributions within a decibel. The use of 
integration considerably facilitated manual reduction 
of the records. Studies of short-term statistics were 
carried out separately. 

The average signal voltage thus measured was 
equated to the rms value for computation of received 
power. Earlier tests [Boggs and Hekimian, 1953] 
showed that using the time average of the signal in 
this way may introduce up to 3 db or more error in 
estimation of absolute power, but that little error is 
made in relative power ratios. 


3. Results 
3.1. Short-Term Signal Statistics 


It was considered important to establish whether 
the statistical description of the short-term behaviors 
of the signals were sufficiently alike throughout the 
frequency band to warrant the assumption: (a) that 
the same propagation process was being observed at 





all frequencies and, to be assured (b) that the 
receiver was integrating signals alike at all fre- 
quencies. A number of fast-response recordings 
(receiver bandwidth 300 c/s, recorder response 25 c/s) 
were made of the fading signal received simulta- 
neously at all frequencies. 

Figure 3 shows cumulative distributions of signal 
envelopes at all frequencies, for typical daytime and 
nighttime conditions. The data shown were ob- 
tained from 10-min sampling periods. It is seen 
from the figure that: 


a. The distributions are similar at all frequencies. 

b. The distributions approximate a Rayleigh 
distribution during the midday hours, and are some- 
what peaked at the “high-signal’”’ end during the 
late night and early morning hours. 


Wheelon [1960] has recently discussed the amplitude 
distribution for a received signal containing effects 
of many meteor bursts, and obtains theoretical 
distributions similar to those shown for the night 
hours in figure 3. 

Median crossings per second of positive slope, N, 
were obtained from these records. These fading 
rates were observed to vary over a wide range of 
values. For short periods (of the order of 30 sec) 
the value of N may be from 0.2 ¢/s to 1.0 ¢/s at 30 
Me’/s, with corresponding values of 1.0 ¢/s to 3 ¢/s 
at 108 Me/s. When signals are influenced by 
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Figure 3. Cumulative distribution of fading envelope at five frequencies for ten minute periods; 
typical daytime and nighttime observations. 
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sporadic-# ionization, the fading rate may remain 
very low for several minutes. Figure 4 shows typical 
values of N versus carrier frequency, averaged over 
10-min periods. 

Most of the observations show N proportional to 
carrier frequency to a power between 3/4 and 5/4. 
One would expect to observe linear dependence of 
fading rate on carrier frequency if the drift speed of 
the scattering region were high compared to the 
speed of change of structure of irregularities. Linear 
dependence of N and good correlation of fading 
structure at all the frequencies would be expected if 
a single intense scatterer in motion, such as a meteor, 
were propagating the dominant signal, or if the 
signal were made up of an ensemble of such reflec- 
tions. Departure from linear dependence of N on 
frequency and poor frequency correlation almost cer- 
tainly implies fading dominated by changing structure 
within the scattering region. 


3.2. Hourly System Loss 


Continuous recording, using the receiver integra- 
tion discussed in section 2 above, was carried out 
from October 1, 1958 through September 30, 1959. 
Hourly median values of the average signal trace 
were scaled for each hour of observation to obtain 
system loss values at the five frequencies. 
































2.0 T T T . . STS T T T me ee | 
1655-1700 gg 
JAN 13 4 JUNE 277 
FS _S\\80-1200 /9 
p JUNE 27 
Pies © - ie 4 7 
. al a ae 0000-0010 7 
wore ; ’ coe JUNE 28 4 
a 71200-1210 ie 1 e 
g 05+ 7 FEBS Zi q 4 
fo} 
¢ + i a 4 
803+ + + : 
« 
— 
> 
T0.2F 41 - 4 
DAYTIME NIGHT TIME 
1.0 L a ose we ee L ews 
2030 50.7000 200 202 30 50 670 ~=— 100 200 
FREQUENCY , Mc/s 
FicurE 4. Observations of fading rate versus frequency; 


median crossings (one sign of slope) observed in ten minute 
period; typical daytime and nighttime observations. 


Cumulative distributions of the hourly system loss 
values at all frequencies are plotted in figure 5 for 
4 months of data representing conditions observed 
during summer, winter, and equinoxes. The solid 
curves represent the scatter signal data, with hours 
of identifiable sporadic-F eliminated, and including 


some hourly values interpolated through short 
periods. One notes that the frequency dependence 





is a little lower at the weak-signal end of the dis- 
tribution. This is discussed later. 
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for months representative of winter, equinox and summer. 
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Figure 6 shows the characteristic diurnal variation 
of the scatter signal at all frequencies for the same 
months as above. Monthly median values of the 
system loss for each hour of the day are plotted. One 
notes that the diurnal variation is consistently greater 
at lower frequencies than it is at higher frequencies. 
The effect, especially noticeable in the December and 
June curves, is progressive through the frequency 
range. This corresponds to diurnal variation of the 
frequency dependence which will be shown more 
explicitly later. 


3.3. Frequency Dependence of Hourly System Loss 


Using hourly median values of signal intensity 
observed at the five frequencies, linear least-squares 
regression lines of system loss in decibels versus log 
frequency were computed on an electronic computer 
for each hour’s data for the 1-year observing period. 
The slope of this line is proportional to the hourly 
value of n. Standard error of n was computed for 
each hour to check the goodness of fit of the straight 
line to the data. Median values of 0.30 for March 
1959 and 0.34 for June 1959 are examples of the 
standard error. 

Figure 7 is a representative plot of the observed 
hourly frequency dependence of system loss. Ob- 
served signal intensity at each frequency is plotted 
relative to that observed at the highest frequency, 
108 Me/s. The lines are indicative of individual 
hourly observations but the plotted points in the 
figure are actually median values of all hourly obser- 
vations in July 1959, for hours near midday and 
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FicurE 7. Frequency dependence of system loss; observed 
signal intensity at each frequency plotted relative to that 
observed at 108 Mc/s; values shown in plot represent all 
hourly values near midday and midnight during July 1959. 


near midnight. The spread of hourly observations 
is indicated by the vertical lines through the points, 
the extremities of which represent the values observed 
10 and 90 percent of the time. Because of the good- 
ness of the fit of the points to straight lines at all 
hours, the spread shown is also nearly representative 
of the computed regression lines. 

The actual correspondence of results with this 
“straight line’’ law is the most striking departure of 
the present results from earlier work. Though 
varies from about 6 to 11, with characteristic diurnal 
and seasonal variation, received power is always 
closely represented by the strict power law, p,/p, of", 
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Figure 6. 


Monthly median plots of hourly values of system loss, showing diurnal variation at 


five frequencies in winter, equinox and summer. 
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over the range of frequencies observed. The 
earlier NBS results [Bailey, Bateman and Kirby, 
1955] had been interpreted as giving a lower value 
of n for the lower part of the frequency range (28 
to 50 Me/s) than was obtained for the higher part 
(50 to 108 Me/s). The two pairs of frequencies had 
not been observed simultaneously, and later studies 
of the siting and antenna characteristics also raised 
some questions about the idealness of the radiation 
patterns. It is felt that the careful siting and cali- 
bration of antennas and power in the present experi- 
ment, together with the feature of simultaneous 
comparison at frequencies throughout the whole 
range, fully justifies revision of the earlier conclusion 
and adoption of the idea of a strict power law. It 
thus appears valid to express the turbulent spectrum 
function S(k) as a power relationship k- ‘*~®’. 

The straight-line law also has implications as to 
absorption effects. During the course of the VHF 
measurements, attenuation data over the same path 
were obtained from a continuous 5 Mc/s transmission. 
Sufficient sensitivity was achieved to permit con- 
tinuous day and night observation at the single 
frequency, so that information was obtained on 
diurnal and seasonal variation of nondeviative 
D-region attenuation. The observed values of 
attenuation at 5 Mc/s were extrapolated to VHF 
by the relation 

1 
—  —-—— e 
att © 


This indicated that if the scatter signal suffered full 
D-region attenuation, the noon values at 30 Me/s 
would be reduced by 3 to 4 db and noticeable curva- 
ture would be introduced into the scaling law. Alter- 
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natively, if the observed straight line law included 
the effect of daytime absorption, then the absence 
of absorption at night would introduce curvature. 
But the experimental results gave a greater value 
of n by day than at night and a straight line at all 
times except during abnormal absorption events such 
as are associated with solar flares. It appears that 
under normal conditions the VHF scatter signal does 
not encounter measurable attenuation and therefore 
must not traverse a significant part of the absorbing 
region. 


3.4. Diurnal and Seasonal Variation of Frequency 
Dependence 


Figure 8 shows that the value of n, varying over 
a range from about 6 to 11, is characteristically 
maximum during the daytime. This figure shows 
values of n exceeded 90 percent, 50 percent, and 10 
percent of the days for each of four months repre- 
senting all seasons. Diurnally n usually varies +1, 
though July shows less. 

Figure 9 shows systematic seasonal variation of 
n. The median of all daytime values (hours 0900 
to 1400) and nighttime values (hours 2100 to 0200) 
is plotted for each month; despite the spread of values 
in each month, indicated by the vertical lines, the 
seasonal variation seems clear. A summer minimum 
is distinct and there is a suggestion in the daytime 
values of a secondary winter minimum. 

It is of interest to inquire whether the summer 
minimum represents absorption effects; i.e., whether 
greater summer D-region absorption tends to reduce 
the net frequency exponent of the scatter signals. 
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Fieure 8. Diurnal variation of n during four seasons. 
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However, the seasonal depression of values of 7 is 
about as great at night as during the day, and n is 
greater during the day. These observations again 
tend to rule out absorption as a factor. It seems 
more likely that the diurnal and seasonal variation 
of n represents changing relative roles of meteoric 
reflections and turbulence. 

A frequency exponent of 7 has been estimated for 
scatter signals that are composed entirely of meteoric 
echoes [Villard, Eshleman, Manning, and Peterson, 
1955]. Thus, the value of 7 could be expected to 
prevail at times when the signal is mostly dominated 
by meteoric reflections. Higher values are provided 





by turbulence theory. Comparison of n for broad- 
beam and narrow-beam systems discussed in 3.6 also 
strengthens the idea that variation of the relative 
role of meteor reflections may account for part of 
the diurnal and seasonal variation of n. It is also 
possible that the S(k) for turbulence varies with 
time of day and season, though the theories have 
not suggested this. 

Figure 10 shows cumulative distributions of hourly 
values of n observed during each month of the pro- 
gram. For the year, the median value of all obser- 
vations of n was 8. Ninety-five percent of the 
observations exceeded 6.9 and five percent exceeded 
cae 

3.5. Weak-Signal Conditions 


It is of interest to the design engineer to know 
the frequency dependence of system loss during con- 
ditions of weakest signals, which generally govern 
the limitations of the system capacity. Values of n 
associated with signal levels exceeded large percent- 
ages of the time can be obtained by comparing the 
cumulative distributions of system loss for a year’s 
observations at each frequency. The following 
values are obtained: 








Obtained for system loss cor- 
responding to signal level 
exceeded: 





| 
Value of n | 
| 
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8.1 | 50% of time 
7. 6 90% 
7. 4 | 99 % 
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Ficure 10. Cumulative distribution of values of n observed during each month of the year. 
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A value n=7.4 may thus be taken for engineering 
computations concerned with weak signal conditions. 
If galactic noise is limiting at the receiver, the net 
signal-to-noise ratio may be taken to vary as ~f~°, 
allowing for the frequency dependence of the noise. 


3.6. Effect of Beam Width 


A parallel experimental system, using horizontal 
Yagi transmitting and receiving antennas with 
azimuthal beam widths of 56°, was available over 
the same path as the narrow-beam systems. This 
broad-beam system provided the opportunity for 
simultaneous frequency-dependence observations at 
two frequencies, near 30 and 50 Me/s. 

Figure 11 compares monthly median values of 7 
at all hours of the day for the two systems for 
September 1959. The broad-beam system values 
of n range from } to 1 lower than the narrow-beam 
values. The diurnal variation of n is more pronounced 
for the broad-beam system, and distinct minimums 
are observed during early evening and early morning 
hours. 

These results are interpreted as additional evidence 
that a lower value of n is associated with the meteoric 
contribution more effectively than with the narrow- 
beam system directed along the Great Circle, espe- 
cially during early morning and evening hours, for 
geometrical reasons discussed in the literature. More 
noticeable absorption effects in the broad-beam sys- 
tem might also tend to reduce the effective exponent, 
as transmission off the Great Circle must traverse 
a longer path in the D-region. These absorption 
effects have manifest themselves during at least one 
intense absorption event associated with a solar flare, 
when signals in the broad-beam system were attenu- 
ated much more severely than the narrow-beam 
system signals. 


3.7. Abnormal D Ionization 


A number of observations were made during ab- 
normal events such as daylight occurrences of sudden 
ionospheric disturbances associated with solar flares, 
and magnetic disturbances. 

During SID’s attenuation was characteristically 
observed at the lowest frequencies, and an increase 
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Ficure 11. Comparison of values and diurnal variation of n 
for narrow and broad beam systems (measurements made 
simultaneously). 








of scattered signal intensity was observed at the 
highest frequencies. Attenuation of as much as 30 db 
has been observed at 30 Mc/s, but the signal never 
dropped below the threshold of observability. 
Usually, the attenuation was overtaken by increased 
scattering to produce a net enhancement within 10 
to 30 min. Attenuation is rarely observed at 50 
Mc/s; usually an enhancement is noticeable. 

Figure 12 shows signal behavior at five frequencies 
during an SID of March 29, 1958. Absorption is 
severe at 30 Mc/s, and is noticeable at 41 Me/s. 
Note that a signal enhancement is observed at all 
frequencies, somewhat greater and later at the 
lower frequencies. The initial rise in signal intensity 
at 108 Me/s is essentially simultaneous with the onset 
of the absorption at 30 Me/s. As time progresses 
increased scattering develops at the lower frequen- 
cies, eventually overtaking the attenuation effects 
at 30 Me/s. 

It is of interest that during the same event for 
which the narrow-beam system data is given in figure 
12, approximately twice as much attenuation (in 
decibels) was observed on the broad-beam Yagi 
system at 30 Me/s, and transmission at frequencies 
as high as 74 Mc/s suffered noticeable attenuation. 
This is attributed to the greater path length through 
the absorbing D-region for the energy propagated 
off the Great Circle; it is also likely that dependence 
of the broad beam propagation on meteoric echoes 
at higher ionospheric heights than the absorbing D- 
region tends to increase the effective attenuation. 


4. Summary 


The frequency dependence of system loss has been 
found to follow a strict power law p,/paaf" over 
the frequency range of the experiment, 30 to 108 
Mc/s, at all times except during abnormal ionospheric 
absorption. Thus, it seems the spectrum of turbu- 
lence S(k) may be characterized by an \"~? depend- 
ence. Figure 13 compares the distribution of n with 
theoretical values from Wheelon’s [1960a] summary 
of spectral models. 

The value of n varies diurnally and seasonally ; 90 
percent of all observations lie in the range 7 to 
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Figure 12. Signal behavior at five frequencies during an 
intense ionospheric absorption event associated with a solar 
flare, March 29, 1958. 
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Ficure 13. Cumulative distribution of all hourly values of n 
measured during the year, with theoretical values noted for 
comparison. 


93. The median value of the observations is 8. 
These values, for scaled antenna systems, include 
the f? factor, which is the effect of variation of re- 
ceiving antenna aperture with frequency. Daytime 
exponents are characteristically greater than night- 
time values by 1 or 2; summer values are lower than 
those observed at other seasons, though winter 
daytime values are also lower than equinoctial day- 
time values. For system design, a value of 7.4 
may be taken to characterize weakest signal periods. 
Since galactic noise varies as about /~*”, the re- 
sultant signal/noise ratio varies as f~°. 

Systematic temporal variation of m is accounted 
for by (a) diurnal and seasonal variation in the 
relative contributions of meteoric reflections and 
turbulent scattering, as well as (b) possible diurnal 
and seasonal variation of S(k). The latter pos- 
sibility has not yet been suggested by the theories. 





It is concluded that normal ionospheric absorp- 
tion does not affect the scattered signal, so the 
scattering would appear to occur mostly below the 
absorbing region. During abnormal intense absorp- 
tion events associated with solar flares or other 
disturbances, attenuation is observed at frequencies 
below 40 Mc/s. Enhancement of signal is character- 
istic at all frequencies during such events, but is 
observed at the lowest frequencies only after the 
initial absorption effects no longer dominate. 

A broad-beam system yields lower values of n 
than are obtained with a narrow-beam system. In 
the experiment, the system having an azimuthal 
beam width of 56° gave values of n lower by 3 tol 
than the system using 6° beam width. The broad- 
beam system exhibited more pronounced minimums 
in the diurnal variation of n during the early morning 
and early evening hours, and susceptibility to greater 
attenuation during absorption events. 


A. D. Wheelon and K. L. Bowles contributed valuable 
suggestions on the initiation and conduct of the program. 
K. F. Florman designed the antenna systems and P. P. 
Viezbicke supervised the calibration of antennas and meas- 
urement of radiation patterns. J. L. Green supervised the 
transmitting station; R. P. Fitzpatrick and C. H. Johnson 
were principal observers and data analysts at the receiving 
station. 
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In calculation of the theoretical scattering coefficient for a terrain, previous authors 
° ° ° ° 2 
tentatively assumed the normalized autocovariance function p(r)=e-4™ for the ground 


elevation as a function of distance from a given point. 
were made using maps with contours ranging from one to twenty-five feet. 
in curves of p(r), which are approximated by exp (—|]r|/B). 


Recently autocorrelation studies 
These resulted 
The theoretical scattering 


cross section (go) of many such terrains can be expressed as 





op=—4y2 2 hs (= )e —4k? o 260826 
sin @ 


_ (4k?) "(cos?@) n+l (1) 





— 2 ote B’ sin? @+-n? 3/2” 


where a, A, k, and 6 are standard deviation of the target terrain, wavelength, wavenumber 


(2x/d) and the angle of incidence respectively. 


to k, the above expression becomes 


40? 
o=— ~ (@cot46) 


AB 


For the case where 1/B is small as compared 


for@<0°. (2) 


These expressions, when normalized, are in agreement with experimental results of other 


authors. 


It is also noteworthy that the results obtained with an acoustic simulator model 
compared very well with this theoretical expression. 
that the ground is conducting and has random elevation variations. 


This work is based on the property 
Theoretical results 


calculated on the basis of varying ground impedance rather than its elevation sre also in 


agreement with this expression. 
1. Introduction 


In recent years, the calculation of backscatter 
from a rough surface, with very obvious extension 
and application to the radar return from the moon, 
has attracted considerable attention. Here we 
follow the basic approach of Davies [1954] as modi- 
fied by Moore [1957] and Cooper [1958]. The 
modified Kirchhoff-Huygens’ principle is employed 
in the calculation using modified spherical variables 
of integration. The ground model used assumes 
“‘facets,’’ of variable size and height above mean 
ground level, whose position is described statistically. 

Instead of assuming a correlation model for 
terrain, as is usually done, the present approach 
used contour maps of different terrain samples in 
the United States to calculate the terrain-elevation 
autocovariance and other statistical properties. 
The overall approach to the problem is an approxi- 
mation, but the results so obtained are very reason- 
able indeed, insofar as the comparison with terrain 
return and moon-echo data reported in various 
publications [Briggs, 1960; Hughes, 1960; Nielson, 
1960] is concerned. 


1 This work was sponsored by the Naval Ordnance Test Station, China Lake, 
Calif., under contract No. N123(60530)18138A. 





2. Statistical Properties of Terrain 


Various autocovariance [Ament, 1953; Davies, 
1954; Moore, 1957b; Daniels, 1960] functions for 
rough terrain have so far been assumed. In some 
instances [Cohen, 1948] an effort was made to 
derive an approximate expression for the autocorrela- 
tion from actual data on terrain elevation. A 
general expression [Norton, 1960] for space correla- 
tion function p (7) describing the random variation 
of the refractive index over space is given here: 


o(r)= Ea 5 Ge xi) K,(7): (2.1) 


l,=characteristic scale 
u=—constant 
T (u)=gamma function 


K, (7_)=modified Bessel function of the second kind. 
0 


where 


It seems to cover some of the most commonly as- 
sumed [Wheelon, 1959] expressions for correlation 
functions. 

A search of the literature shows that small-scale 
perturbations of the terrain elevation have not so far 
been used to calculate an experimental autocovari- 
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Location of terrain samples and the surface 
characteristic constant, B 


TABLE 1. 























| | 
Average | Standard| Sampling} Surface 
Sample General location | elevation} deviation| interval | characteristic 
| constant, (1/B) 
ft it it 
ae Lawrence, Kans___-- 853.5 5. 63 30 0. 1193 x 10-3 
ee Tees Nos, Ariz__-- 5254 49 5 1. 3213 x 10-3 
ae | Turtle Mountain, } 
| eee | 1990 8.2 5 0. 7070 x 10-3 
Daicancene Gila River, Ariz_-_.-- 2617 13.6 4 . 9236 x 10-8 
5......-.| White River, Ariz 5028 5.18 5 | 1. 3610 x 10-8 
ES Mountain Park, 
_ ee 7530 1030 62.5 | 0. 2729 x 10-3 
, ro Sandia Park, | | 
ta aetna 6780 250 20 =| . 9860 x 10-3 





ance function. For this reason, contour maps for 
seven different terrains in the United States were 
selected, as listed in table 1. This set included 
relatively flat land, rolling plains and some rugged 
mountain areas. The average elevation (above sea 
Jevel) and standard deviation of these samples varied 
from 850 to 7,530 ft and 5 to 1,030 ft, respectively. 
Three distinctly different but random lines were 
drawn on each of these contour maps. Along these 
lines, the elevation of the terrain was read to within 
one tenth of a foot (by interpolation) on 1-ft contour 
maps, and within 5 ft on 20- and 25-ft contour maps. 
Using a reasonable compromise between resolution 
and confidence level [Blackman and Tukey, 1958], 
and due to the unavailability of less-than-one-foot 
contour maps, horizontal sampling intervals of 5, 20, 
30, and 62.5 ft were used. The number of points 
for a subsample varied from 58 to 610. The auto- 
covariance was calculated using a CRC-102 com- 
puter. The resulting curves were averaged for each 
type of terrain. A simple theoretical approximation 
of these curves for high confidence level portions was 
found to be 


p(r) ~e7'7/B (2.2) 
where 


B=characteristic constant 
r=distance between points 


Seven cases for the experimental autocovariance are 
shown in figure 1. 

In figure 1 the autocorrelation curves for samples 
1 and 3 seem to slope off faster than the rest. In case 
of sample 1 this is caused by the lag distance for a 
given number of lags being greater than for samples 
2 through 5 and 7. In view of this, if all these 
curves were plotted on a semilog paper (as in fig. 1), 
for a given lag distance rather than for a fixed number 
of lags, these would plot approximately as straight 
lines over limited distances. These two curves can 
also be approximated as p(r)=e7'"'* cos cr, where c 
was found to be 0.2188 and 0.7931 deg/ft for samples 
1 and 3, respectively. For near vertical incidence, 
the distances involved are rather short, and therefore 
cos cr is approximately unity. It is felt that the 
apparent cosine factor in the autocovariance appears 
because of the relatively short length of the sample 
used in this case. This seems to indicate that the 
autocovariance function p(r) often varies as e~!"!/2. 
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It is known [Cohen, 1948; Isakovitch, 1952] that 
the slope of p(r) at r equal to zero should be zero 
unless the ground elevation function has infinite 
slopes. In this instance, p(r) has a slope of —1/B 
which is in most cases believed to be of the order of 
10-* (see table 1). Considering the minute size (in 
the general range of 10-* to 10-* m) of the earth parti- 
cles, this may be quite a reasonable approximation to 
the exact description of the terrain roughness. 

It was also found that the ground elevations 
generally are normally distributed. For a certain 
type of terrain, the ground elevation data along a 
random line is most probably a random function 
belonging to a large, approximately stationary en- 
semble. It is further suspected that the expressions 
obtained for the scales of roughness covered in this 
study when reduced by a scale factor, might well 
describe extremely minute variations of ground 
roughness as seen by very high frequency waves. 
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3. Scattering Coefficient 


The expressions for vector waves, in radar return, 
although easy to set up are very difficult to evaluate. 
For this reason the Kirchhoff-Huygens’ principle has 
been applied to scalar waves. This leads to difficult 
integration problems and other such complexities, so 
certain simplifying assumptions were made. Similar 
assumptions were made by others [Davies, 1954; 
Moore, 1957b]. These are as follows: 

(1) No portion of the ground is shielded from the 
incident radiation. 

(2) The ground is considered to be a perfect 
conductor. 

(3) The magnitudes of the surface currents are of 
the same order as those of a plane reflector, but the 
phase varies in a random manner, depending on the 
height of a particular point. 

(4) The reradiation from a particular small area 
on the ground is isotropic. 

(5) The antenna gain G is essentially uniform for 
—Oo<0<, and is zero outside of this range. 

(6) The results of this scalar wave approach are 
an approximation to those for the vector waves, as 
most field vectors would be nearly parallel to the 
surface involved 

The variables ¢ and ¢’ were replaced with y and 
vy’, where y’ is the sum of y and a as shown in figure 2. 
This change of variable was made following Davies’ 
[1954] work, in order to make the integration a little 
less complex. 
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Ficure 2. Geometry of the problem, showing interrelation 
between various variables of integration. 





The radar equation as applied to pulse radar 
[Moore and Williams, 1957a], gives the average re- 
ceived power P, from many scatterers as 


apes 
x2 Pr (1-25) G?(0)o0(8,A, . . -) 
P= ae dR. (3.1) 


32m R 





where 


P;=power transmitted 
R=range 
oo=scattering coefficient per unit area (assumed 
independent of ¢) 
t=time 
c=velocity of propagation (velocity of light) 
G=antenna gain (assumed independent of ¢) 
A=wavelength 
6, ¢=angles as shown on figure 2. 


The field E at 0’ is obtained by applying Huygens’ 


principle to the modified geometry of figure 2, and 
can be written as 


E=| <p (I.n)'/? cos 6e~ #*F GA, (3.2) 
where 


I,=P7G/(4rR?) 


n=intrinsic impedance of free space 


/ | 
dA=area element on the ground (=nak a) 


One arrives at the following expression for E by 
applying (3.2) to a rough surface, 


z=! 4 (% coté Pea\ 
NJ-erJ-0 R \4aR? 
4 


exp [j2k(R—6 cos 6)]dydR (3.3) 





where 
k=wave number (=) 


7=pulse width 
26,—=beam width of the antenna 
5=6(Ro+s,y+a)=elevation of ground height 
above mean ground level at a point located 
at range Ry+s and modified azimuth angle 
+a. 
hinds to the projection of a point at range R, 
on the mean ground level (see fig. 3). 


Equation (3.3) takes into account the variations in 
phase and neglects the other effects of the changes in 
range. This is based on the approximation that the 
percentage changes in range are not appreciable for 
near-vertical incidence. The received power is 


EE* Gx 
n 47 





P=; Re (3.4) 
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A typical point, at a random height 5) above mean 
ground level. 


FIGuReE 3. 


where E* stands for complex conjugate of E and Ze 
is the symbol for ‘real part of’ the expression follow- 
ing it. The difference in ranges of two points on 
ground located at (R,0,~) and (R’,6’,~’) can be ex- 
pressed as follows with the help of figure 3, 


R—R’=[R,—4(Ro, W) cos 6]—[Rj—4(Ro, ¥’) cos 6’) 
=~—s+[6(Ro+s, ¥+a) cos 0’—85(Ro, ¥) cos 6] 
=—s+(6—6)) cos 0 (3.5) 

where @~0’ (assumed) 

5=6(Ro+s, +a) 
5o=5(Ro, P) 
Ro2zRots 
The substitution of (3.3) in (3.4) results in 


{ PG? cos 6 cos 6’ 


P.= J Rip? e~IRR-RY TAA’. 


30,2 Fee 








(3.6) 


where the primed quantities refer to a point at (R’, 
6’, v’). 

The average received power can now be expressed 
in terms of R, s, y, and @ by substituting (3.5) in 
(3.6) as 


al ‘ime 1 agp PG? cot 6 cot 6’ 
os 32?) _, 20ye) —(Ro—h) Ry-— RR’ 


exp [j2k{s+(6:—4) cos 0}|dRdsdyda (3.7) 








where primed quantities refer to values at a point 
different from that for unprimed quantities. The 
limits — 26) to 24 on the variable y are used to cover 
the entire 360° of the azimuth angle. 





It is apparent from the statistical studies of contour 
maps that the probability density function for ter- 
rain elevation above the mean (say 6) can be written 
as 


p= e~87/(20), (3.8) 


1 
(2x (2ro?)'/2 


The conditional (normal) density function for the 
elevations 6 given 4) at two points on the ground, a 
distance 7 apart, is given by 


1 _ (= inp)? 
P(6;b0,") =p aq — Bae Beto. 





(3.9) 


The expression exp [j2k{s+(é)—6) cos @}] in (3.7) 
is first averaged over 6 using (3.8) and the result is 


exp [j2ks—a*k? cos? 6+72ké9(1— >) cos 6], 
(3.10) 


where 
a’?=207(1—p’). 


Now (3.10) is averaged over 4) using (3.9), and it 
becomes, 
exp [j2ks—4k’o?(1—p) cos? 6]. (3.11 


One can deduce from Davies’ [1954] approximation 
r? = R*o’?+-s* cosec’ 6 (see fig. 2) that 


s? conee" Q 


r~Rat (3.12) 


It can be now shown that the integration of (3.7) 
after substituting (3.11), (3.12), and (2.2) results in 
an expression, which when rearranged, becomes 








P=pol “) Se ( @ ) 
v 32n? R3 B) \sin 6 
« 2 2\ 2 n+1 
on ttre? cost o x __ (Aka?) *(cos? 6)"*"__| ap (3.13) 





rigs rtf 2k sin? 0+ | 


Here it has again been assumed that R’ and 6’ are 
approximately equal to R and @ respectively, for 
near-vertical incidence. From the comparison of 
(3.1) with (3.13) it is clear that the scattering coef- 
ficient, 09, is given by 








_ 4y2nB? 6 
ae sin 6 
— 48% cos? 9 > (4k?c?)"(cos? 6)"+! 
. nai (n—1)I (2k? B? sin? 6+-n?}*/2 (3.14) 
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The value of 1/B is far less than k for nearly smooth 
surfaces, and hence the scattering coefficient for such 
surfaces can be approximated from (3.14) as 


2 
oo (for nearly smooth surfaces) ae (@cot*@) (3.15) 
for 640°. 


4. Experimental Verification 


For nearly smooth surfaces, the surface character- 
istic constant 1/B~0, and (3.15) gives the scattering 
coefficient. This result compares very closely with 
published results [Nielson, 1960] for new ice as shown 
in table 2. 

For rough (not nearly smooth) surfaces, (3.14) 
describes the relationship of the scattering coefficient 
o» and other variables such as the angle of incidence 
6, wavelength A, standard deviation o« and surface 
covariance constant B, etc. Two curves of the scat- 
tering coefficient o) versus @ for each of the three 
values of \/B, 0.1, 0.5, and 1.0 for o/A equal to 0.05, 
and 0.1 are shown in figure 4. It may be noticed 
that as the surface becomes rougher, or as \/B 
increases for a specified , the scattering coefficient 
curve becomes flatter, showing the relative impor- 
tance of the contribution of the power return from 
the surface at angles other than those near zero. 


As expected, when the surface becomes: smoother or 


1/B decreases, the received power seems to come 
primarily from near-zero angles. These curves are 
quite similar to those recently published [Campbell, 
1959; Dye, 1959; Edison, 1960]. The experimental 
data [Nielson, 1960] on desert and new ice also seems 
to follow the pattern of these theoretical curves de- 
scribed above. 











TABLE 2. Comparison of theoretical versus experimental scat- 
tering coefficient 
(Normalized) 

6° | ao Theoretical | oo Experimental 

= | é 
30 1.000 1.000 
40 0. 291 0. 308 
50 . 088 . 089 
60 . 022 . 021 
70 . 004 


. 004 | 





Similar results ? obtained by an acoustic simulator 
at the University of New Mexico also verify these 
theoretical conclusions. It is interesting to note 
that an experimental expression (o.=0e"' for 
6 in radians) for the scattering coefficient [Hughes, 
1960] of the Moon, verified by photographic astro- 
nomical [Briggs, 1960] calculations, is a very close 
fit to the graph of eq (3.14) for ¢/A=0.1 and \/B=1.0 
in the range of incidence angles of 3 to 14 deg. The 
authors believe that other values of B and o may 
well be appropriate, as only a limited set has been 
tried to date. 


2A. R. Edison and R. K. Moore, Preliminary report on an acoustic simulator 
for investigation of backscatter of E.M. waves, unpublished report of the Univer- 
sity of New Mexico, Albuquerque, N. Mex. (1961). 
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5. Conclusions 


The scattering coefficient (o)) for nearly smooth 
surfaces is inversely proportional to the wavelength, 
but varies directly with (c*), (@ cot* @) and 1/B, where 
o, 6, B are standard deviation, angle of incidence, 
and the terrain characteristic constant respectively. 
For rough surfaces it has a negative exponential 
a cos? 0 

oN 
times a constant. The surface characteristic con- 
stants B and o can be calculated from the radar 
return data. Although approximate, the theo- 
retical results agree well with the experimental data; 
and therefore, suggest the usefulness of the approach. 
The application of these results may be extended to 
the moon-echo data, with proper corrections for 
Faraday and liberation effects, etc. This investi- 
gation has established that for near-vertical inci- 
dence, the normalized autocovariance for the terrain 
elevation is more often of the exponential form exp 
(—|r|/B) rather than the Gaussian form, exp (—r?/B). 
The former may well be more appropriate for finer 
terrain irregularities than those considered in this 


factor, where the exponent is made up of 
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paper. It may also be representative of the normal- 
ized autocovariance of the moon surface. An exact 
theoretical, but usable, expression for the scattering 
coefficient has been very reasonably approximated 
by the results of this study. 
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Estimates of the mutual interference expected to occur between the ground terminals 
of space communications systems and surface point-to-point systems are presented in a 


fashion suitable for engineering applications. 


These estimates are obtained from recently 


developed methods for predicting the transmission loss over tropospheric paths in terms of 
parameters such as geographic separation, elevation angle of the antenna, antenna patterns 


and frequency. 
under suitable conditions. 


1. Introduction 


Many of the problems of mutual interference 
between proposed satellite communications systems 
and surface communications systems can be investi- 
gated in terms of presently available theories. The 
prediction method used here is that developed by 
Rice, Longley, and Norton [1959], for tropospheric 
scatter propagation. The method is designed to 
predict the median value of hourly median basic 
transmission loss, Ly, and to give a distribution of 
the medians about this value; the median basic 
transmission loss is defined in terms of the ratio of 
the power, p,, radiated from the transmitting 
antenna to the available power, p,, at the receiving 
antenna, when isotropic antennas are used at both 
eads of the path. The method is applicable for 
any path configuration. The median transmission 
loss, L, is obtained from this by subtracting the path 
antenna gain, G,, [Hartman and Wilkerson, 1959] 
from the basic transmission loss. It should be noted 
that these methods give excellent agreement with 
data measured over a frequency range from 100 
Mc/s to 10,000 Me/s. Most of these data were 
taken over conventional scatter paths with the 
antennas directed at the horizon in the great circle 
plane, although some preliminary data obtained 
using elevated beams are shown later in the paper. 

Tolerable interference will depend not only on the 
relative levels of the desired and undesired signals, 
but also on the type of modulation used for both 
of these signals and in most cases must be deter- 
mined experimentally, in some cases subjectively 
(for example, subjective rating of television pic- 
tures). Further, neither desired nor interfering 
signals will be steady signals, but will be represented 
by a distribution of values. Because of this, satis- 
factory operation should usually be specified in 
terms of an_ hourly-median-desired-signal to an 
hourly-median-undesired-signal ratio required to 


1 This is a résumé of work sponsored by the Joint Technical Advisory Commit- 
tee of the Electronic Industries Association and the Institute of Radio Engineers. 





It is concluded that these systems can share the same frequency assignment 


provide a given grade of service or better for some 
percentage of the hours. A discussion of this prob- 
lem is given by Norton [1959]; and Miller, Ports, 
and Savage [1960]. However, in view of the many 
unknowns the actual problem has been somewhat 
simplified by relating the hourly-median-interfering- 
signal power to the noise power of the receiver under 
consideration. Thus the question answered here is 
what separation distances and antenna elevation 
angles are required so that the interfering hourly- 
median-signal power is equal to or less than the 
receiver noise power for a given percent of the hours. 


2. Path Parameters 


The following values of the path parameters are 
used here to illustrate the methods of a more general 
treatment [JTAC, 1961]. 

It is assumed that the satellite ground terminal 
has an antenna with the general characteristics of a 
60-ft parabolic dish, and that the surface system 
antennas are at most 10-ft dishes. The satellite 
terminal antenna beam will be directed above the 
horizon at an angle y as shown in figure 1, and the 








Ficure 1, 


Geometry. 
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surface system antenna wil] be directed at the hori- 
zon. Both antennas are to be located 30 ft above a 
smooth spherical earth. Because antenna heights, 
terrain effects, and diffraction become increasingly 
important for shorter path lengths, a distance of 100 
miles is chosen as a minimum for which calculations 
of this type will give representative values of trans- 
mission loss. At any distance a careful selection of 
the terminal site is necessary, and in some cases it 
will be possible to use terrain to give protection at 
distances shorter than those shown in this paper. 
If the specific path parameters are known, the trans- 
mission loss can be predicted accurately for any dis- 
tance using the methods of the references given 
above [Hartman and Wilkerson, 1959; Rice, Longley, 
and Norton, 1959]. Both antenna axes are centered 
in the great circle plane. Although this situation 
usually results in the maximum interference, it may 
be impossible to avoid it if the satellite systems share 
frequencies with existing systems. 


FREQUENCY | Gc 
DISTANCE 100 MILES 
10 ft DISH-60 ft DISH, SIDE ONLY 


10 ft DISH-60 ft DISH 
LOBE ONLY 


—- 





Figure 2 shows values of hourly median transmis- 
sion loss which will be exceeded by 99.9 percent of 
the hourly medians for a 10-ft dish at one end, and 
the main beam, a side lobe with the same half-power 
beamwidth as the main beam and 20 db down from 
the main beam, and an isotropic portion of the 
antenna pattern for a 60-ft dish at the other end. 
These are plotted versus the elevation angle, y, for 
the main beam of the 60-ft dish. Examples of this 
type of calculation for median values are shown in 
figure 3. Also shown in this figure are the results of 
some preliminary measurements. These data points 
represent 5-minute medians measured over a 165- 
mile test path from Boulder to Haswell, Colo. The 
pertinent parameters for this path are as follows: 
Frequency, 409.9 Me/s, the angular distances with 
both antennas directed at their respective horizons 
is 33 milliradians, the fixed transmitting antenna is a 
14-ft parabolic dish, and the receiving antenna being 
elevated is a 60-ft dish. 


FREQUENCY 2 Gc 
DISTANCE 100 MILES 


1oft DISH-6oft 





DISTANCE 150 MILES 


of DISH DISH, SIDE LOBE ONLY ~ 


ft DISH, MAIN ONLY 


FREQUENCY 6.5 GC 
DISTANCE 100 MILES 


L(99.9%),db 


10 ft ft DISH, 
SIDE LOBE ONL’ 


DISTANCE 150 MILES 


10 ft DISH-60 ft 











10 


DISTANCE 150 MILES 
10 ft DISH-60 ft SIDE LOBE ONLY 


loft DISH-60ft 





FREQUENCY 10 Gc 
DISTANCE 100 MILES 


DISH 
SIDE LOBE ONLY 


ft DISH-6Oft 
MAIN BEAM ONLY 





DISTANCE 150 MILES 





¥, deg 
Fiaure 2. Transmission loss exceeded 99.9 percent of the time. 
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Figure 3. Boulder-Haswell Test Path, transmission loss 


exceeded 50 percent of the time. 


The following procedure is used to determine the 
required separation. distance for which the undesired 
signal is less than the receiver noise power at the 
receiver intermediate frequency for 99.9 percent of 
the hours. 

Compute the value of transmission loss required to 
make the interfering signal equal to the noise power. 
This value is 


Lreq=P.—L.—-M—P, 


where P, is the power, in decibels, of the interfering 
transmitter, I, represents coupling losses in the 
system and may include transmission line losses, 
cross-polarization losses, etc., M is a term to allow 
for the situation in which the transmitted energy is 
spread over a frequency band different from that 
which will be accepted by the receiver. It will 
depend on the type of modulation as well as the 
bandwidths, but an estimate is made here by letting 
M=10 log b,/b, where b, and 5, are the transmitter 
and receiver bandwidths, respectively. P, is the 
available noise power in decibels at the receiver 
terminals and is given by 


P,=NF-+10 log kT)b=10 log kT.6 


where NF is the effective noise figure [Norton, 1959; 
Hogg, 1960] of the receiving system, & is Boltzmann’s 
constant, 6 is the noise bandwidth of the receiver in 
cycles per second, J) is a reference temperature, 
288.48 °K, and T, is the effective noise temperature 
of the receiver and antenna combination. 7’, should 
not be confused with the effective antenna tempera- 





ture. This value of Zyeq is then compared with the 
curves for the frequency under consideration to 
determine if Lreg< L’ where L’ is the smallest of the 
transmission loss values using the main beam or the 
side lobe or the isotropic part of the antenna pattern. 


3. Examples 


1. Consider first the case of interference from the 
earth terminal transmitter to a point-to-point relay 
receiver. The receiver noise power is determined as 
above assuming a bandwidth of 20 Me/s and a noise 
figure of 10 db. Then 


P,=NF-+10 log b—204=—121 dbw. 


Assume a transmitter power of 1 kw, coupling losses 
of 4 db, including line and polarization losses, and a 
transmitted bandwidth of 20 Me/s. Then 


Lng P~Lg~M—P,=147 db. 


This transmission loss value is now used with figure 2 
to determine whether the above condition, Ly.q< L’ 
can be met. It is clear that for a separation distance 
of 100 miles and for all four frequencies the trans- 
mission loss from the isotropic portion of the trans- 
mitting antenna is less than the 147 db required loss. 
For a separation of 150 miles the energy transmitted 
by way of the main beam becomes significant and 
the required conditions can be met at all frequencies 
if the main beam of the 60-ft antenna is elevated 
above the horizon by approximately 5°. 

This analysis assumes that the antenna of the 
surface point-to-point system is directed toward the 
earth terminal. Rotating this antenna so that the 
main beam does not point directly toward the earth 
terminal would increase the transmission loss. When 
the pattern of the 10-ft dish is at the isotropic level 
in the direction of the earth terminal the “‘isotropic- 
isotropic” curves of figure 2 will apply rather than 
the “‘10-ft dish-isotropic” curves. The transmission 
loss via main and side lobes will of course increase by 
approximately the same amount. Therefore, the 
requirement could be met at 100 miles and at all 
frequencies by orienting the 10-ft dish to reduce the 
interfering signal by 10 db at 1 Ge and by lesser 
amounts at the other frequencies. 

Other earth terminal to satellite systems might 
require more or less transmitter power. Suggested 
systems [FCC hearing] range from approximately 3 
to 36 dbw. The effect of changing this parameter 
in the above example would give a range of Lyreg 
from 120 to 153 db at the 99.9 percent level. 

If instead of a 60-ft dish, a 120-ft dish or a 30-ft 
dish were used for the transmitter, the only notice- 
able effect would be in the calculations for L(p) for 
the main beam. In the former case L(p) for the 
main beam would be decreased by less than 6 db and 
in the latter case, L(p) would be increased by approxi- 
mately 6db. This can alter the conclusions only by 
increasing (or decreasing) the required elevation 
angle of the transmitter by at most 2°. 
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2. Consider next the case of interference from a 
point-to-point microwave relay transmitter to the 
earth-terminal receiver of a satellite system. It is 
assumed that a wide deviation FM system is used 
in the satellite to earth link, so that the earth ter- 
minal receiver has an RF bandwidth of 100 Me/s. 
With the use of a low noise antenna and maser 
amplifier, we assume that the effective noise tem- 
perature of the receiver is reduced to 30 °K. 


P,,=10 log T.b—228.6=—134 dbw. 


Compute the transmission loss, Lyeq, assuming a 
transmitter power of 1 w, coupling losses of 4 db 
and a transmitted bandwidth of 20 Mc/s. Then 


Lnea=P ;—Lp—M—P,=137. 


This value of Lyeq is compared with figure 2 to 
determine whether the requirement can be met. It 
is seen that under these assumptions the point-to- 
point antenna could not be directed toward the 
earth-terminal at a distance of 100 miles. As in the 
previous example the curves indicate that the re- 
quirements could be met by not allowing the 10-ft 
dish to be directed toward the earth terminal. 

The effects of varying the transmitter antenna 
size and the transmitted power were noted in ex- 
ample 1. Similar statements can be made about 
varying the receiver antenna size and transmitted 
power in this example. Values of the power for 
typical point-to-point microwave relays range from 
—3 dbw to +7 dbw. 


4. Conclusions 


Our theoretical analysis indicates that space 
communications systems and surface systems of the 
conventional microwave relay type can share the 
same frequencies if care is used in locating the possible 
interfering sources. As seen in the examples, 








separation distances of from 100 miles to 150 miles 
usually suffice, and under ideal conditions, distances 
of less than 100 miles could give adequate protection. 
Estimates have been made for other systems such 
as high-powered radar and these indicate that 
harmful interference should not be experienced if 
the radar and earth terminals are separated by 
distances of 500 miles or more [JTAC, 1961]. 

The data shown in figure 3 represent samples 
during times when aircraft were not present in the 
propagation path, and although the data agree with 
the predicted values for the stated conditions, the 
estimates may not be accurate for some paths where 
aircraft are present. 


The measurements reported here were made under 
the supervision of A. F. Barghausen. 
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During 1959 and 1960 long-term transmission loss measurements were performed over 
a 223 kilometer path in Eastern Colorado using frequencies of 100 and 751 megacycles per 
second. This path intersects Pikes Peak which forms a knife-edge type obstacle visible from 
both terminals. The transmission loss measurements have been analyzed in terms of diurnal 
and seasonal variations in hourly medians and in instantaneous levels. As expected, results 
show that the long-term fading range is substantially less than expected for tropospheric 
scatter paths of comparable length. Transmission loss levels were in general agreement with 
predicted knife-edge diffraction propagation when allowance is made for rounding of the 
knife edge. A technique for estimating long-term fading ranges is presented and the results 
are in good agreement with observations. Short-term variations in some cases resemble the 
space-wave fadeouts commonly observed on within-the-horizon paths, although other phe- 
nomena may contribute to the fading. Since the foreground terrain was rough, there was 
no evidence of direct and ground-reflected lobe structure. 

In most cases comparatively high correlation exists between signals received simultane- 
ously on two antennas with 8.3 and 14 meters vertical separation. These separations were 
chosen as being representative for practical space diversity systems designed for eliminating 
the effects of fading arising from direct and ground-reflected phase interference phenomena. 
The comparatively high correlation observed suggests that space diversity will be relatively 
less successful in mountain obstacle paths with rough terrain near the terminals than on 
tropospheric scatter paths or on line-of-sight paths over smooth terrain. 

The enhancement of field strength associated with propagation over mountain ridges 
may cause concern in applications where mountains are being counted on to shield unwanted 
radio waves. Some radio astronomy installations have been located in mountain valleys for 
this reason, and it is possible that obstacle-gain effects may aggravate rather than alleviate 


interference. 
1. Introduction 


It is generally known that mountain ridges can act 
as diffracting knife-edges in point-to-point communi- 
cation circuits on VHF and UHF frequencies. The 
mathematical approach has been developed using the 
Fresnel-Kirchhoff diffraction theory [Schelleng, Bur- 
rows, and Ferrell, 1933]. Substantial reductions in 
transmission loss as compared to scatter circuits of 
similar length, and relatively fading-free signals were 
reported by several workers, both in this country and 
abroad [Dickson, Egli, Herbstreit, and Wickizer, 
1953; Kirby, Dougherty, and McQuate, 1955; Kono 
and Hirai, 1955; Crysdale, Day, Cook, Psutka, and 
Robillard, 1957; Nishikori, Kurihara, Fukushima, 
and Ikeda, 1957; Hirai, Fujii and Saito, 1958]. More 
recently, additional work in Japan [Kakita, Iwai, and 
Teiri, 1959], and so far unpublished data from Alas- 
kan knife-edge diffraction paths! have shown that 
some paths exhibit considerable fading at carrier fre- 
quencies in the 900 Me/s range. Much of the time 
the fading encountered resembles the ‘‘space-wave 
fadeout”’ phenomena observed commonly on within- 
the-horizon paths at comparable frequencies [Bean, 
1954]. It consists of relatively deep and slow fades, 
and the occurrence of these fades may show marked 





! Private communication from Western Electric Co. 
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diurnal trends and correlation with meteorological 
phenomena as well as with simultaneous ducting on 
paths just beyond the horizon. 

To obtain a better understanding of the long-term 
signal and fading characteristics, NBS established a 
test path in Colorado to study knife-edge type dif- 
fraction phenomena. This path extends in a roughly 
north-south direction along the Front Range of the 
Colorado Rockies, and uses Pikes Peak as the dif- 
fracting knife-edge. An outline map of the path with 
the locations of the terminals at Beulah and Table 
Mesa, Colo., is shown on figure 1. Pertinent path 
and equipment data are given in table 1, below. The 
path profile is shown on figure 2, based on an equiva- 
lent earth radius of 9,000 km. This is 1.41 times the 
actual earth’s radius and corresponds to a surface 
refractivity of 316 N-units. 

Operation on 751 Mc/s commenced in December 
1959, and terminated in September 1960. Data 
were recorded continuously during 5-day periods at 
the rate of one each month with more frequent opera- 
tion in June and August. Operation on 100 Me/s 
started in August 1960, on the same basis. During 
two 5-day periods in August 1960, an additional 
site was operated on top of Pikes Peak, using hori- 
zontal half-wave dipole receiving antennas on both 
frequencies. These antennas were mounted about 
6 m above ground. The length of this auxiliary 
line-of-sight path was 77 km. 
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FigureE,1. Outline map for Colorado obstacle gain path. 


TaBLeE 1. Path and equipment data 
Pikes Peak obstacle gain pcth 








EEE OAL | 223 km. 
Angular distance. _- Ch ghikaletin east _ 64.4 milliradians. 





Terminal elevations above mean sea 
level: 

Trenmmitter...........- 

REGEN a ee 


wiatiieies (Beulah) 1905 m. 
(Table Mesa) 1666 m. 


Height of obstruction above mean sea 
level: 
Eee ae ae ee 4292 m. 


Height of obstruction above mean 
terminal elevation ---_....-.......-. 2507 m. 


Antenna heights above ground: 


Transmitter_____- pote ... 7.3m for 751 Me/s, 12.4 m for 100 Me/s. 


I ek on cstuak oe 8.2, 16.5 and 22.3 m for 751 Me/s, 13.7 m 
for 100 Me/s. 
Polarization _ --_-_-- . Horizontal. 
Antennas: 751 Me/s 100 Mc/s 
Transmitting _.. 4.3m Dish Yagi 
Receiving. -- . 3m Dish Yagi 


Modulation ___- Ee AEE 


ag EE PEAK,ELEVATION 4292.7m 


_ continuous wave. 
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Figure 2. Terrain profile of Colorado obstacle gain propagation 
path. 
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The receiver outputs were recorded on strip charts 
and multichannel magnetic tape. Transmitter 
power output values were also recorded on strip 
charts, so that a continuous record of the trans- 
mitter power is available for both frequencies. The 
receivers were calibrated by signal generators 
checked against laboratory standards. 


2. Characteristics of the Received Signal 
and Variation of Hourly Median Values 


2.1. Signal Characteristics 


Figure 3 shows sample recordings of both fre- 
quencies at the two sites taken on the morning of 
August 10, 1960. The 751 Mc/s record at the Table 
Mesa receiving site shows more rapid small signal 
variations superimposed on a much slower fading 
pattern, which includes a fade in excess of 15 db 
lasting several minutes. There appears to be sub- 
stantial correlation between the signal recorded 
simultaneously on the two antennas on 751 mce/s. 
The slow fading pattern is absent on the Table 
Mesa 100 Me/s record; only the more rapid varia- 
tions are present, and they are superimposed on a 
virtually constant signal. It is believed that these 
rapid signal variations of small amplitude are due to 
a component of scattered field superimposed on the 
more slowly varying diffracted field. For the line- 
of-sight path with the terminal on the top of Pikes 
Peak this scattered component is not noticeable, as 
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Figure 3. Simultaneous chart recordings for Colorado obstacle 
gain and line-of-sight propagation paths, August 10, 1960. 
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expected. However, the slower variations are still 
present on 751 Me/s, and they do not appear cor- 
related with the slow fading observed at the Table 
Mesa site. 

The slow signal variations which are similar to 
the “‘space-wave fadeouts’’ commonly observed on 
line-of-sight paths will be discussed in more detail 
later on. 


2.2. Variations of Hourly Medians 


A typical sample of diurnal variations of hourly 
medians is shown on figure 4 for the period August 
8 to 12, 1960, at the Table Mesa receiving site. Each 
open circle denotes an hourly median of basic trans- 
mission loss for the hour given by the abscissa. The 
distribution of the points and the median curve 
shows no pronounced diurnal trend. This is in 
contrast to (smooth earth) diffraction or scatter 
paths of similar length where usually a substantial 
diurnal trend of hourly median values is observed. 

Similarly, the absence of a pronounced seasonal 
trend is shown by figure 5. Here, the median of all 
hourly medians recorded during each 5-day period 
is plotted versus time. The maximum variation is 
about 5 db. The overall medians obtained for the 
two antennas appear well correlated, and there seems 
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Figbas 4. Diurnal variations of hourly median transmission 
6 loss levels for August 8-12, 1960. 
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Ficure 5. Seasonal trend of overall weekly medians, 751 Me/s. 


to be no significant difference between the results 
obtained at the two heights used for the corner 
reflector antenna. During the three recording 
periods in June when the corner reflector antenna 
was at the 22.3 m elevation, the overall medians 
follow the trend given by the rest of the data, when 
the 16.5 m height was used. 

The 100 Me/s data at Table Mesa and the Pikes 
Peak data were not included in this study, because 
they are only available for August and September 
not enough to study seasonal trends. 

It should be recognized that the time periods on 
figure 5 are consecutive as far as recorded data are 
concerned. Weeks with no recordings have been 
omitted. 

Diurnal and seasonal variations may also be 
studied by obtaining distributions of hourly medians 
for eight time blocks as previously defined [Norton, 
Rice, and Vogler, 1955]. For this purpose, the 
hours of the day are divided into four blocks for 
summer as well as winter months, as follows: 





TABLE 2. Time blocks 





No. Months | Hours 
Bo occkiwcewaccucntcusceuona| See Oe Eo. watts Soto 0600 to 1300 
SS ee DARED eS LYS IRE Cor) November to April..................-. 1300 to 1800 
, ere ene 2 ee eee a eS 1800 to 2400 
CARS Rees ee eee ey pene Oe 3 eee 0600 to 1300 
aed da aiahctantie ature Fg ee 1300 to 1800 
as Peiareeens i ere 
j eae een ee ee pore ee pe Sree | rae ee Oetoner:.. .. 6.255 acssnesas 0000 to 0600 
Divviisanswneenancmeeeunl | Nevessber te Apefi.._..............-.. 0000 to 0600 





Cumulative distributions of the hourly basic trans- 
mission loss medians are calculated for each indi- 
vidual time block, or for combinations of time 
blocks, and plotted on log-normal graph paper. 
Thus it is possible to obtain an estimate of the range 
of observed hourly medians which is the long-term 
fading range. This will be compared with theoretical 
estimates, or with similar data obtained from other 
paths. 

The cumulative distributions of all measured 
hourly median transmission loss values are shown 
on figure 6 separately for all time blocks, and for 
both 751 Me/s antennas, with the corner reflector 
at the 16.5 m elevation. All curves are bunched 


439 





WINTER SUMMER 












CORNER REFLECTOR 
___AT 165 METERS 









+ 


3. METER DISH 
AT 82 METERS 


210 T 





3. METER 
AT 82 METERS * 


BASIC TRANSMISSION LOSS, db 








220 


oo! 510 60507 9093 99 999 AI ' 50 30507 909 99 939 


PERCENT OF HOURS 


Freure 6. Distributions of hourly median basic transmission 
loss values for Colorado obstacle gain path on 751 Me/s. 


and nearly parallel, indicating very little diurnal or 
seasonal variation. Table 3, below, shows a com- 
parison of the overall median values and the long- 
term fading range. Here the overall median is the 
median of all hourly medians within a time block, 
and the long-term fading range is defined as the 
difference between hourly median values exceeded 
1 percent and 99 percent of all hours. 

Even though the cumulative distributions of hourly 
medians plotted by time blocks on figure 6 appear 
quite similar, table 3 indicates some seasonal effect 
in the long-term fading range. The range is larger 
during summer, and is also larger for the corner re- 
flector at the higher elevation than for the dish. 
There is very little variation in the overall median 
values between time blocks, and the difference be- 
tween the dish and the corner reflector appears to be 
slightly higher in summer. 


TABLE 3. 


Comparison of winter and summer results 





Overall median Long-term fading 





decibels range, decibels 
Difference (1% to 99%) 
Time block No. | ee ee 5 ee I 
| medians | 
3m dish Corner | db | 3mdish | Corner 
reflector | reflector 
Winter: 
8 (0000 to 0600) _ 198. 4 193.2 | 5.2 8.9 14.8 
1 (0600 to 1300) - 197.6 192.6 | 5.0 8.9 10.5 
2 (1300 to 1800). 198. 3 192.8 5.5 12.4 14.0 
3 (1800 to 2400) _ 197.9 192.7 5.2 | 11.5 11.6 
Summer: 
7 (0000 to 0600) - 198. 1 192.7 5.4 | 13.5 16.2 
4 (0600 to 1300) - 197.6 191.4 6.2 | 13.8 14.4 
5 (1300 to 1800) _ 197.3 189.9 7.4 | 11.4 10.7 
6 (1800 to 2400) _ 198. 5 192.3 6.2 13.0 12.9 





2.3. Comparative Studies of Median Distributiors 


Although tropospheric scatter is not primarily 
involved in these experiments, it is of interest to com- 
pare first the measured distributions with theoreti- 
cally expected values and with data from tropospheric 
scatter paths baving similar path parameters. For 
this purpose, two hypothetical paths and the meas- 
ured results obtained from two actual tropospheric 
scatter paths were used. The two hypothetical 
paths were considered to have the same antenna 
height as the actual path. One was assumed to have 
an angular distance which a tropospheric scatter path 
of the same length would have over a smooth earth, 
and the other would have a path length correspond- 
ing to a smooth earth path having the same angular 





distance as the obstacle gain path. The actually 
measured paths were selected to correspond closely 
to these assumptions in path geometry and fre- 
9 Yad Table 4, below, shows the path constants 
used : 


TABLE 4. Path constants for comparative study 














Path | Angular 
distance | distance 
km Milliradians 
Ee CER | 223 24.0 
BN III ois cone nnctan cen ncesurusen amas | 566 | 64.4 
Actual paths: | | 
Pikes Peak obstacle gain (751 Mc/s)_-------------- 223 64.4 
Cheyenne Mt. to Garden City (1046 Me/s)--------| 364 | 26.7 
59. 8 





Cheyenne Mt. to Anthony (1046 Me/s)_.-..------- | 633 | 





Overall time block medians and the long-term 
fading ranges for the hypothetical paths were deter- 
mined from a consideration of the expected long- 
term variation of hourly medians. These distri- 
butions are obtained from the empirical V(p, 6) 
curves, which were determined from a large amount 
of experimental data representing a number of differ- 
ent geographical and climatic regions within the 
United States [Rice, Longley, and Norton, 1959]. 
The time block medians and long-term fading ranges 
given there are measured relative to a calculated 
reference value of basic transmission loss, which is 
the median of all hourly median values in time 
block 2 (winter afternoon). The curves are assumed 
to be valid for all frequencies in the 70 to 4000 Me/s 
range. 

A comparison of the long-term fading range 
determined in this manner for the hypothetical paths 
with the ones measured on the Pikes Peak obstacle 
gain path is shown as table 5 for all time blocks 
separately, for all hours during summer and winter, 
and for all hours of the year. 


TaBLeE 5. Comparison of hypothetical tropospheric propagation 
paths data with measured data from Pikes Peak obstacle 
gain path 





Long term fading range, db (1% to 99%) 








Time block | Hypothetical | Hypothetical | Pikes Peak 
| path A, 223 | path B, 566 | obstacle ony 
| km,@=24 | m, 6=64.4 | path, 223 km, 
| milliradians | milliradians 6=64.4 millira- 
| | dians (8m dish) 
i: Ee SSCS nee 42.1 21.1 | 8.9 
ed Rae aa s PRS 37.6 13.9 | 8.9 
"EES See ee 32.9 14.0 12.4 
a EN SS ES eee 39. 2 16.3 | 11.5 
{ 
fe ae | 42.2 21.4 | 13.5 
Gaia ee 37.7 28.1 13.8 
es ere Se oe i oo ewan 25.1 18. 2 | 11.4 
eS ee a en 36.5 19.0 | 13.0 
All hours, winter... -.---.---- } 38.8 19.3 | 11,7 
All hours, summer-___-------- 41.9 26.0 } 13.4 
7 26.4 | 14.3 


eR ee ee ee | 39. 





It is seen that the actually observed long-term 
fading ranges were considerably smaller than the ones 
predicted on the basis of scatter propagation, al- 
though the fading predicted for a scatter path 


440 








i 
; 
9 





having the same angular distance is a somewhat 
better estimate. 

In order to provide a more valid basis for esti- 
mating long-term fading ranges and the distribution 
of hourly medians for knife-edge diffraction paths, 
use is made of a method suggested by K. A. Norton. 
The results of this method as applied to the Pikes 
Peak obstacle gain path will be compared with the 
actual measurements in what follows. 

For the purpose of studying variations in hourly 
median field strength or transmission loss values a 
knife-edge diffraction path is considered to consist 
of two Jine-of-sight paths in tandem. The diffract- 
ing knife-edge constitutes a common terminal for 
the two paths. If the hourly median signal for 
each of the two line-of-sight paths is subject to 
random variations, V, and V2, then the random 
variation V expected on the mountain obstacle 
path would be V=V,+V2, where V, V;, and V2, 
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Figure 7. Predicted and measured cumulative distributions of 
basic transmission loss hourly medians, December 1959 to 
September 1960. 


44] 





are random variables expressed in decibels. The re- 
sulting cumulative distribution of hourly medians for 
the entire path may be determined by the convolu- 
tion of the distributions of V; and V, [Davenport 
and Root, 1958]. 

The cumulative distributions of hourly medians 
for the two individual line-of-sight paths may be 
calewated by empirical formulas as functions of the 
path distance, the effective antenna heights, the 
carrier frequency, the horizon elevation angles, and 
season, summer or winter. The methods used are 
unpublished empirical methods communicated to 
the authors by P. L. Rice. This material will be 
described in detail in forthcoming NBS Technical 
Notes and in a U.S. Air Force Technical Order 
dealing with design procedures for tropospheric 
communication circuits. Results are shown on 
figure 7, where the empirical distributions are 
compared with the actually measured distributions. 
Comparison is made for the data on 751 Me/s 
obtained using the 3-m dish separately for all hours 
during the summer and during the winter months, 
and for the 100 Me/s data for the summer months. 

Calculation of the empirical distributions of 
hourly medians shown on figure 7 includes determi- 
nation of the overall median reference value by 
Fresnel diffraction and ray optics methods which 
will be discussed in more detail in section 3 below. 
What is primarily important here is the comparison 
of the slopes of the actually measured and the 
empirically predicted distribution curve shown on 
figure 7, rather than their absolute position on the 
basic transmission loss scale. 

The comparison of measured and predicted fading 
ranges from figure 7 is summarized in table 6 below. 
TABLE 6. Comparison of actually measured fading ranges 

— empirical predictions for Pikes Peak obstacle gain 

pat 











Long-term fading range, 
in decibels (1% to 99%) 











Actually Empirical 
measured prediction 
751 Me/s: ; 
Son ies oon cn cn cadets 11.7 11.5 
a ee 13.4 16.2 
100 Me/s: 
Aa) eee, eae Lo 7.0 9.0 





Table 6 shows that the empirical method just 
described provides much better estimates of the 
long-term fading range than the assumption of a 
trodospheric scatter path of the same linear or 
angular distance. This conclusion applies to this 
specific knife-edge diffraction path, and should be 
tested against data obtained from other knife- 
edge paths in order to arrive at an estimate of the 
uncertainty inherent in the prediction process. 
However, agreement of the predicted 1 to 99 percent 
long-term fading ranges is well within the sampling 
error of the empirical data on which this prediction 
is based and this suggests that the proposed method 
of predicting long-term fading ranges is a reliable 
one. 
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Fiaure 8. Comparison of predicted and measured diurnal and 
seasonal variations of the long-term fading range. 


The data in tables 5 and 6 have been plotted on 
figure 8 to obtain a better picture of the comparison 
between predicted and measured long-term fading 
ranges and their diurnal trends. As the empirical 
prediction method just mentioned does not dis- 
tinguish between individual time blocks, it is 
represented by a straight line for all summer and 
winter time blocks, respectively. In spite of the 
large discrepancy between the measured fading 
ranges and the ones predicted on the basis of scatter 
propagation, it is of interest to note that the pre- 
dicted diurnal trend is evident in the measured 
data for the summer time blocks, although to a 
considerably lesser degree. This is not the case for 
the winter time blocks due to the relatively large 
measured Time Block 2 long-term fading range. 

An additional comparison of the measurement 
results was made with similar data obtained from 
two actually measured tropospheric propagation 
paths which represent a geographical location similar 
to the Pikes Peak obstacle gain path. In addition 
thereto, these paths correspond most closely to the 
hypothetical paths described above in path distance 
and angular distance, and at the same time, repre- 
sent a similar carrier frequency. The hourly median 
distributions for these two paths were taken from 
unpublished reports. A description of the paths 
and the various experiments conducted have been 
published by NBS [Barsis, Herbstreit, and Hornberg, 
1955]. 
| Figure 9 compares the distribution of hourly 
medians measured during the summer months over 
the Pikes Peak obstacie gain path, with the two 
paths between Cheyenne Mountain and points in 
the eastern plains. It is seen that the long-term 


fading range is smallest for the obstacle gain path; 
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Figure 9. Cumulative distributions of hourly median basic 
transmission loss values for Colorado-Kansas propagation paths. 
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it is quite large for the Garden City path which 
represents the same linear patb distance, and slightly 
larger for the Anthony path which represents the 
same angular distance, but a substantially greater 
linear path distance. Table 7 shows the comparison 
of overall median values and ranges taken from figure 
9. 


TaBLE 7. Comparison of long tropospheric propagation path 











results 
| Basic | 
Fre- transm. | Long-term 
Path quency Time period No. of loss fading 
Me/s hours overall range, db 
| median | (1% to 99%) 
db | 
eT SO SE, Ae en en eee } —— 
| | | 
Pikes Peak 751 | Summer 1960..._| 613 | 198.7 | 13.4 
obstacle gain. | | | 
Cheyenne Mt. 1046 | Summer 1952... 1646 208. & | 30.5 
Garden City. 
Cheyenne Mt. 1046 | Summer 1952 184 225. 5 | 21.5 
Anthony. (July to 
August). | 





The overall medians have been added for com- 
pleteness of the information. It is quite clear that 
transmission loss tends to be smaller for a knife-edge 
diffraction paiii as compared to tropospheric scatter 
paths of similar length. Perhaps even more sig- 
nificant is the substantial reduction in the long-term 
fading range. Although there are fewer hours of 
measurements available for the Cheyenne Mountain- 
Anthony path, table 6 shows that a iy come 
scatter path of comparable linear distance has a 









































substantially greater long-term fading range than 
an obstacle gain path, whereas a scatter path of 
similar angular distance has only a moderately 
greater long-term fading range. This agrees well 
with the results obtained by the application of 
empirical methods for determining fading range— 
the method based on scatter propagation paths 
yielding substantially greater fading ranges than 
the method based on line-of-sight paths. 


2.4. Frequency Dependence of the Distributions of 
Hourly Medians 


Figure 10 shows a comparison of all measured 
data on 100 Me/s and 751 Me/s during the two 
weekly periods when the receiver site on Pikes Peak 
was operated. These are overall summer-time dis- 
tributions; they have not been separated into time 
blocks because of the comparatively few hours 
available. As expected, the long-term fading range 
is smaller for the line-of-sight path to the top of the 
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FIGuRE 10. Measured distributions of hourly median basic 
transmission loss values for August 8-12 and 22-26, 1960. 
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peak than for the obstacle gain path. Overall 
median and ranges for this period are given in table 8 
below. 


of hourly medians 





TaBLE 8. Data from distributions 


Long-term 
| fading range, db 
(1% to 99%) 


Basic 
transmission 
loss overall 
Median, db 


Path and frequency | 









Beulah-Pikes Peak, 100 Mc/s_.-_-.-------- 128.2 5.7 
Beulah-Pikes Peak, 751 Me/ a il 132.0 6.3 
Beulah-Table Mesa, 100 Me/s-_---_------- 165.3 7.0 
Beulah-Table Mesa, 751 Me/s (dish)__--.-- 197.0 13.0 
Beulah-Table Mesa, 751 Me/s (corner) - - -- 191.6 | 14.7 








This tabulation shows a substantial frequency 
dependence of the fading range for the obstacle gain 
path—more than doubled when comparing 100 Me/s 
with either one of the 751 Me/s antennas. There is 
not nearly that much frequency dependence on the 
line-of-sight path. It has already been shown that 
the V(p, @) distributions are not appropriate for use 
in predicting the performance of knife-edge diffrac- 
tion paths, and that much better agreement is 
obtained by the use of a different empirical method 
taking into account carrier frequency and horizon 
elevation angles. 


3. Height Gain Studies 


Previous measurements made of obstacle diffrac- 
tion over Pikes Peak [Kirby, Dougherty, and 
McQuate, 1955] showed marked agreement between 
the predicted and observed variation of transmission 
loss with height. These measurements were made 
using a shorter path and with a relatively smooth 
terrain area in the immediate foreground of the 
antenna whose height was varied. 

The 751 Me/s equipment installed on Table Mesa 
permitted the collection of similar data for antenna 
heights up to 100 ft using a telescopic tower with 
dual receivers and recorders. However, the Table 
Mesa site differs from the receiving site previously 
used in that the terrain is very rough except for a 
small flat area immediately in front of the receiver 
site, which is only large enough for approximately 
one-third of the first Fresnel-zone ellipse. It has 
been previously shown that the cumulative distribu- 
tions of hourly basic transmission loss medians for 
the two antenna heights on 751 Me/s have approxi- 
mately the same dispersion and the medians are 
separated by approximately 5 db in winter, and up 
to 7.4 db in summer. Using a small stored-program 
electronic computer the height distribution of field 
strength or basic transmission loss was calculated 
by methods of geometrical optics using the vector 
sum of four rays [Dickson, Egli, Herbstreit, and 
Wickizer, 1953; Crysdale, Day, Cook, Psutka, and 
Robillard, 1957]. Each of the rays has a magnitude 
and phase value determined by reflections from the 
ground and diffraction by the knife-edge type 
obstacle. The reflecting planes were assumed to be 
parallel to a smooth earth’s surface, and to pass 
through the bases of the antenna towers. Reflection 








coefficients of 0.5 and 0.9 were assumed, with a phase 
shift of 180° occurring at the reflecting planes. The 
attenuation and phase-shift at the obstacle were 
calculated using the asymptotic expansions for 
Fresnel integrals [Schelleng, Burrows, and Ferrell, 
1933]. This was justified due to the magnitude of 
the parameters involved. For all calculations a 
refractivity value N,=290 was assumed, which 
corresponds to average winter afternoon for the 
geographical region where the measurements were 
made. Changes in the N, value within the range 
240 to 320 N-units did not materially affect the 
relative height dependence of basic transmission loss. 

A height-gain run was performed in the morning 
hours of March 25, 1960. For each position of the 
corner reflector antenna (mounted on the telescopic 
tower) the basic transmission loss was determined 
for both receiving antennas, so that the values 
measured could be corrected for changes recorded 
on the fixed antennas. The results are plotted 
versus receiving antenna height on the subsequent 
figures in order to permit comparison with calculated 
values using various parameters. 

Figure 11 shows the comparison of the measured 
height dependence with calculated values, using reflec- 
tion coefficients of 0.9 and 0.5. There is no agree- 
ment in the magnitude of the transmission loss or 
in the number of lobes. The measured data show 
only one lobe with a maximum at approximately 
21 m above ground, whereas the theoretical calcula- 
tions predict 4 to 5 lobes within the height interval 
measured. A change in the reflection coefficient 
only changes the magnitude of the maxima and 
minima. Similarly, any change in the assumed 
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Figure 11. Comparison of calculated and measured height-gain 
curves for Colorado obstacle gain path. 





elevation of the transmitting antenna has little or 
no effect upon the position and number of lobes at 
the receiver, but can modify the transmission loss 
approximately the same amount at all heights. 
Height-gain experiments using the transmitting 
antenna were not performed. 

Studies were also made assuming the cross section 
along the profile of the refracting knife-edge to be 
rounded. The terrain profile cross section through 
the top of Pikes Peak can be approximated by a 
circular are with a radius of roughly 11.5 km. Fol- 
lowing, in general, the rounded knife-edge procedure 
[Wait and Conda, 1959], the resultant diffracted 
field as a function of receiving antenna height was 
calculated again as the vector sum of four rays, 
each of which is affected by the rounded obstacle 
and by reflections from the terrain in the foreground 
of the antennas. These results are shown on figure 
12, together with the measured height-gain run 
and the theoretical curve for the ideal knife-edge. 
Both theoretical curves are based on a reflection 
coefficient of 0.9. It is seen that only the absolute 
values of transmission loss are affected by “round- 
ing’’ the knife-edge; the number and position of the 
lobes with antenna height remains the same. Addi- 
tional calculations show that the change in trans- 
mission loss depends on the radius of curvature 
assumed; as this is somewhat arbitrary for this 
terrain configuration, the assumption of 11.5 km 
depends on how much of the profile the arc is fitted to. 

Figure 11 shows that the long-term median values 
of transmission loss observed on the dish at 8.2 
m and the corner reflector at 16.5 m have approxi- 
mately the same difference as is indicated by the 
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Figure 12, Comparison of calculated and measured height-gain 
curves for Colorado obstacle gain path. 
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Ficure 13. Transverse profile through Pikes Peak as seen from transmitter site. 


measured height gain curve. For the two different 
heights used for the higher antenna (corner reflector 
at 16.5 and 22.3 m) there is little difference in the 
values of transmission loss observed. This tends to 
confirm the similarity of long-term transmission loss 
values observed at the two heights as shown in 
figure 5. 

It is quite clear from the above discussion that 
the basic transmission loss dependence on antenna 
height cannot be represented in this case by the 
simple four-ray knife-edge model, even if the knife- 
edge is assumed to be rounded. The reason becomes 
evident if the transverse profile through the obstacle 
is studied, which is a terrain profile drawn through 
the obstacle at right angles to the propagation path. 
Figure 13 is a photograph of the visual horizon as 
seen from the transmitting antenna with various 
visible ridge lines emphasized. A grid is super- 
imposed which shows the angles subtended from 
the transmitting site. It is seen that the trans- 
mitting antenna beam, with a 7.7° half-power beam 
width sees a substantial portion of the ridge line 
on either side of the principal obstacle. Similar 
conditions exist if the transverse profile is viewed 
from the receiving terminal. Thus field contribu- 
tions may exist from the ridges and other terrain 





features between the terminals and Pikes Peak 


illuminated by the antenna beams. It has so far 
not been possible to make a quantitative study of 
these contributions in order to use them for field 
computations. The effect of a transverse triangular 
mountain profile was studied by Matsuo in [1950]. 
Numerical results for a specific example were given 
but no comparison was made with actual measured 
results. 


4. Vertical Space Diversity and Fading 
Characteristics 


A detailed study of the short-term fading char- 
acteristics observed on this path is contained in 
another paper [Barsis and Johnson, 1961]. A num- 
ber of possible mechanisms appear responsible for 
the fading. These include those associated with 
prolonged space-wave fadeouts, which resemble 
phenomena observed on within-the-horizon paths, 
as well as those associated with scattering com- 
ponents. Other mechanisms appear to be due to 
the terrain configuration and may be the cause of 
an intermediate type of fading, characterized by 
fadeouts on the order of 5 db which last less than a 
few minutes. 

In view of the apparent positive correlation 
between the 751 Me/s signals recorded simulta- 
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neously on the two receiving antennas, studies were 
performed to determine the correlation coefficient 
between instantaneous envelope values for various 
sample periods. For this purpose, the transmission 
loss information recorded on magnetic tape was 
digitized and fed into a computer program. The 
tape contained adequate timing signals to insure 
proper evaluation without time lag. Table 9, 
below, contains the time and duration of each 
sample, the separation of the receiving antennas, 
and the resulting correlation coefficient. For pur- 
pose of digitizing, the magnetic tape was sampled 
once per second. The correlation coefficients listed 
apply to basic transmission loss in decibels; a check 
revealed that results were not materially different 
if receiver input voltages (proportional to field 
strength) were correlated. 


TABLE 9. Results of envelope correlation studies for Pikes 
Peak obstacle gain path 








Antenna |Correlation 


Date and time of sample (1960) | Length I \ 
| spacing coefficient 








Minutes Meters | 
14. 










June 13, 1745 to 1800 15 4.1 | 0. 491 
June 13, 2120 to 2125 5 14.1 a 
June 13, 2200 to 2220 20 14.1 -— 
June 14, 1010 to 102 20 14.1 _ 
June 15, 0738 to 0748. -.------------------- 10 14.1 + 24: 
June 15, 1200 to 1210_--------------------- | 10 | 14.1 | - 033 
June 15, 1800 to 1805----.----------------- | 5 | 14.1 | - 537 
June 16, 0350 to 0400-- _.------------------ | 10 | 14.1 | . 656 
September 15, 1700 to 1710---------------- 10 | 8.3 4 
September 15, 1920 to 1930. - -------------- 10 8.3 - 677 
September 15, 2126 to 2146_--------------- 20 8.3 —_ 
September 16, 0255 to — Fo meine mintei as 4 -s = 
Se 2 ;) OL ares 4 ‘ 

September 16, 1015 to 10 } ~ | 2 | = 


September 16, 1120 to 1130---------------- | 





The samples were picked to represent different 
types of the signal as it appeared on the visual 
recording charts for the two antennas. They are 
not to be considered taken at random in a statistical 
sense; therefore they serve only to illustrate the 
phenomenon, and cannot be used to establish a time 
or spacing dependence, or similar results. They 
show, however, that moderately high correlation 
exists between simultaneous envelope values of the 
signal received on vertically spaced antennas, and 
that the correlation coefficient seems to be higher at 
nighttime and for closer spaced antennas, and that 
these results make the effectiveness of vertical space 
diversity questionable for a path of this type. A 
more complete analysis using many more samples 
evenly distributed over the recording period would 
be necessary to draw final conclusions on the diurnal 
variation and other characteristics of the correlation 
coefficient. : 

A visual representation of the instantaneous en- 
velope correlation may be obtained by applying the 
tape-recorded signals to an z-y plotter. Figure 14 
shows two such samples, and a good qualitative idea 
of the magnitude of the correlation coefficient may 
be obtained thereby. 

Figure 15 shows examples of the correlation co- 
efficient between hourly median values of basic 
transmission loss. In general, hourly medians tend 
to be more highly correlated than instantaneous 
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Ficure 14, Instantaneous envelope correlation for Colorado 
obstacle gain path on 751 Me/s. 


values. The two graphs shown represent the two 
antenna separations used. The correlation coeffi- 
cient is substantially higher for the closer spaced 
antennas (8.3 m). Jt cannot be determined, how- 
ever, if this is entirely due to the spacing, as the two 
samples were taken at different times. The graphs 
indicate the antenna spacing, the number of hours, 
the standard deviations of hourly medians for each 
antenna, and the correlation coefficient. Hourly 
medians determined with the 3 m dish are plotted 
along the y-axis, whereas the corner reflector data 
are plotted along the z-axis. 

Finally, a new analysis of fade durations and signal 
distributions is shown on figure 16. Four periods of 
record ranging from 6 hr to more than 200 hr each 
were analyzed to determine the distribution of fades 
as well as high signal of various levels and duration 
periods. Each curve shown corresponds to one of 
these levels, and the points plotted are the percert- 
ages of the total time the signal remained either 
above or below the indicated level for the duration 
given by the ordinate. The median is the dividing 
line; for each graph curves to the left of the median 
represent fades, and the ones to the right represent 
periods of signal higher than the median. The 
asymptotic extensions of the curves toward a zero 
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Figure 15. Correlation of hourly medians for vertically spaced 
antennas at 751 Me/s. 


fade or signal duration become the distribution of 
instantaneous signal levels for each period in the 
limit. As an example, in the upper right hand 
graph (period August 8 to 12), the curve labelled 
202” indicates the distribution of durations when 
the basic transmission loss exceeded 202 db: 1 per- 
cent of all “fades” to 202 db was of at least 8 min 
duration. 
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Fiaure 16. 751 Mc/s fading characteristics observed on 3 meter 
dish. 





This type of representation permits an analysis of 
fades or short-term variations in terms of basic trans- 
mission loss levels, which may be easily related to 
any arbitrary reference, like the overall median of 
the period considered. This type of presentation 
can be used to determine, at least statistically, the 
characteristics of individual fades. For example it 
might be of interest to know how many fades occurred 
whose duration was between 1 and 2 min. The 
difference between the percentage of the total time 
the signal fell below the specified level for durations 
of 1 and 2 min respectively can be read from this 
typeofpresentation. From this information the num- 
ber of fades within this range of duration can easily 
be determined for the entire period of the analysis. 


5. Conclusions 


The analysis of 751 Me/s data from a long obstacle 
gain path showed principally that the distribution 
of hourly medians of field strength or basic trans- 
mission loss may be approximated by the convolution 
of cumulative distributions for two line-of-sight 
paths in tandem which have the obstacle as a common 
terminal. As expected, these long-term fading char- 
acteristics are quite different from the ones observed 
for tropospheric scatter propagation over comparable 
distances. No significant diurnal variation in propa- 
gation characteristics was observed on the obstacle- 
gain path, and the indication is that seasonal varia- 
tions are small. Due to the limited width of the 
obstacle, the knife-edge approximation by a semi- 
infinite plane, or half-cylinder may not be entirely 
appropriate, although a somewhat arbitrary assump- 
tion of curvature leads to theoretical values of 
transmission loss equivalent to observed values. 

Correlation studies of the signal received on 
vertically spaced antennas suggest that space 
diversity may not be very effective on circuits 
utilizing obstacle gain. Transmission loss recordings 
show slow and deep fades occurring simultaneously 
on antennas vertically separated by 14.1 and 8.8 m. 
Correlation was somewhat lower for the greater 
spacing. 

It is hoped that additional diversity experiments 
can be performed using this path with the antennas 
spaced horizontally in line with and normal to the 
path profile. 

The outstanding advantages of obstacle gain paths 
of this type are that the signal level is significantly 
higher than expected for scatter paths of comparable 
length, and the amplitude of rapid Raleigh-fading 
components is substantially reduced. However, the 
received field strength is lower than values calculated 
using the idealized knife-edge theory. This is 
probably due to the profile through the obstacle, 
which represents a rounded knife-edge, and to 
reflections and contributions from other terrain 
features. 

These and similar measurements show that rela- 
tively high signal strengths are consistently observed 
behind mountain ridges. Thus it is unlikely that 
mountain ranges can always be relied upon to shield 
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potential interference. Special consideration to these 
phenomena should be given in locating radio astron- 
omy or space communication terminals. 


The authors thank all the personnel of the Propa- 
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Messrs. H. T. Dougherty and L. J. Maloney. 
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An experimental study is made of inverted L-, T-, and related transmission-line antennas 


and of antennas with two-, three-, and four-element top loads. 


The apparent measured 


impedance is corrected for end effects; the approximate theoretical impedances are calcu- 


lated for very small heights for which measured values cannot be obtained readily. 


The 


transmission-line antennas studied are shunt-driven with one open and one closed end or 


with two closed ends. 


1. Introduction 
The inverted L-antenna and T-antenna, both 
well known in radio broadcasting, have recently 
found application in telemetry for missiles. More- 
over, a variety of structures that are useful as missile 
antennas may be analyzed by an appropriate super- 
position of different inverted L- or T-antennas with 
‘images. Theories of the inverted L-antenna and 
_T-antenna with images have been developed by 
| Harrison [1958b; 1958d], and by King and Harrison 
[1959a; 1959b]; they have been extended to related 
_ circuits by Harrison [1958a; 1958c]. Special forms 
of these structures with low silhouette so that they 
/may be useful as externally carried missile antennas 
| have been studied by King, Harrison, and Denton 
| [1958] in terms of an equivalent shunt-driven trans- 
| mission line that is terminated at each end in a reactor 
of arbitrary value. 

In order to supplement the extensive theoretical 
work, an experimental study has been made of the 
impedances of antennas of these types including 
especially structures with low silhouette. Antennas 
with multi-element top loads have also been in- 
vestigated. The series of measurements was begun 
with a half dipole over a conducting plane. Antennas 
with low silhouette were approached by successively 
folding over larger and larger sections parallel to the 
ground plane, while the overall length was kept 
/ constant. 
| The driving-point impedances of the antennas 
/were measured with a slotted coaxial line as shown 
‘in figure 1 for the inverted L-antenna; the ground 
screen was of aluminum approximately 10 by 12 
feet in size. A 300 Me/s signal with modulation at 
1000 c/s was used. The detecting system included 
‘a movable voltage probe and a superheterodyne 
receiver. 


| 
t 
| 
| 
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Figure 1. Inverted L-antenna. 


2. Theory of Measurement 


The theory of measurement summarized here 
applies to all of the experimental work described in 
this paper. 


2.1. Distribution Curve Method for Determining 
Apparent Terminal Impedance 


The part of the voltage distribution curve near 
a minimum was plotted and used to locate the 
minimum accurately. The complex terminal func- 
tion O.¢=PsatjJPsa of the load Z,, is defined [King, 
1955, p. 102, eq. 1] as follows: 

6.a=coth!(Z5q/Z,) (1) 
where Z, is the characteristic impedance. The 
phase function %,, is determined experimentally 
from the relation 

®,,=2md/A+-2/2 (2) 
where d is the distance between the location of the 
minimum with a short circuit terminating the 
coaxial line (Z,,=0) and the location of the minimum 











with the actual load Z,, terminating the line. The 
attenuation function p,, is obtained from 
Dp: . 
sin*(6W,/2)—" sin?(BW,/2) 
sinh?p,,= me (3) 
P3 


3} 
1 


where W, and W; are the widths of the distribution 
curve minima at the power levels pj and p3 [King, 
1952]. With a square-law detecting system, the 
power levels p? and p} could be read directly on a 
Ballantine voltmeter. Since p,, is small in all cases, 
the approximation sinh? p,,=p,._ is valid. From the 
experimentally determined values of p., and ®,q, 
the apparent admittance Y,,=G@.a+jB,, may be 
calculated with the following formulas [King, 1955, 
p. 102, eq. 12]: 


sate sinh 2p.4 (4) 
Isa Cosh 2psat+cos 2,, 





ee sin 2®,, 
“cosh 2p,.+ cos 2®,, 





where the normalized admittance is 
Ysa= Y hc Gos +30 s0- 


The apparent impedance of the load is Z,,=1/Y¥ sq. 
Z, is the characteristic impedance of the coaxial 
line; its numerical value is Z,=60.06 ohms. 


2.2. End Correction for the Coaxial Line Driving the 
Antenna Over a Ground Plane 

When an antenna is driven over a ground plane 
by a coaxial line, it is necessary to apply an end 
correction [King, 1955, p. 430] in order to obtain 
the theoretical impedance Z, from the measured 
apparent impedance Z,, terminating the line. This 
correction consists of a capacitance C;, in parallel 
with Z,=1/Y> so that 


Y.= YotjoCr,. (6) 
A curve of Cr,/eoaz as a function of a./a, is given by 
King [1955, p. 435, fig. 22.4]; a. and a, are the 
radii of the outer and inner conductors of the 


coaxial line. For a;/a;=2.762, C7,/eg9a2=—0.25 or 
ne= —0.24210-" farad. It follows from (6) that 


Yo=YoR.= (GotJjbo) R= Ysa—J Bri 
=Gatj(Bsa—wCr,). (7) 
The normalized theoretical impedance is 
2o=Potjto=Zo/R.=1/Yo. (8) 
When @(A+/) is near 1/2, wC7, is negligible compared 


with B,,; when 6(h+/) is near z, the contribution by 
wy, is significant. The theoretical terminal func- 


tions po and ® in @%=po+j® are also computed 
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Since 
Since the end effect 
is purely susceptive and 6) is small, the terminal 
function is given by 


from eq (1) with the subscript 0 in place of sa. 
Go= Gq it follows that po=psa- 


®=27+tan-!bo=r+ bp. (9) 


2.3. End Correction and Theory for Antennas of 
Small Height 


Owing to the finite size of the conductors, it was 
not practical in the experiment to use antennas with 
heights A (fig. 1) less than 2.5 em (@h=2ah/A= 0.161). 
Therefore, the experimentally determined curves of 
po and # were extended with theoretical values. 
These are readily determined since for small heights, 
the inverted L-antenna with its image is essentially 
a balanced section of a two-wire transmission line 
for which the radiation resistance referred to maxi- 
mum current is given by 

R?=30 B2(2h)2=120 Brh?. (10) 
When the maximum current is at the driving point, 
B(h+l)=7/2, the input resistance is R;,=H*. In 
the case of a T-antenna with image, the resistance is 
one-half of R;, for the inverted L-antenna with its 
image since the two equal resistances are in parallel. 
In general, for the antenna with n-element top load 
with image, the resistance is 1/n of R;, for the 
inverted L-antenna with its image since n equal 
resistances are in parallel. When the overall elec- 
trical length of the inverted L-antenna is 7, the 
radiation resistance R* referred to maximum current 
must be transformed through an electrical length 
of line Bw=7/2 to obtain the resistance at the 
driving point. The relevant formula is [King, 1955, 
p. 134, eq 9a, b] Ri,=R2/R* where Ri is the charac- 
teristic resistance of the two-wire line formed by 
the inverted L-antenna and its image. If h denotes 
the height of the antenna over the ground plane 
(one-half the spacing of the two-wire line) and a, is 
the radius of the conductor, FR} is given by 
Pi=120 cosh (h/a;). (11) 
The theoretical value of the terminal function py may 
be calculated from the general formula [King, 1955, 
p. 102, eq 9]. When r2>1, r3>22 this reduces to 
po=1/2 tanh! 2/rp=1/ro. (12) 
For antennas with overall electrical length (height+ 
one element of top load) equal to approximately 
m/2, ro=R*/R, and pp=Ri/R*. When the overall 
electrical length is 7, the maximum current is a 
quarter wavelength from the driving point and 


ro= (R2/R*) (1/R) = Ri/ Ke. (13) 
It follows that for small values of Bh, 
po= R*/R,. (14) 
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The terminal function # for small 6h and when 
8(h+1)=7/2 is obtained from that of a transmission 
line with input reactance X;, in series with the induc- 
tive reactance X,;=wl of the short length of con- 
ductor of length A and radius a, perpendicular to 
the ground plane. The inductance of this straight 
conductor is [King, 1945] 


L=(yoh/m)[(In (4h/a,) + (a,/2h) —1]. (15) 


The input reactance X;, includes a contribution from 
the end effect at the open end. This is approxi- 
mately equivalent to a lumped capacitance Cy, that 
is given by the following formula [King, 1955, p. 367, 
eq 13] 


méo(2h—a) 


Cre 3p (2h/ai)?] 


(16) 


The general formula for transforming a reactance 
Xr-=—1/wC7, from the end of a line of length / to 
its driving point is obtained from King [1955, p. 134, 
eq 9a, b] in the special case when the termination 
is purely reactive. The sum of the normalized 
reactances 2=2;,+2, (normalized with respect to 
Ri, the characteristic impedance of the two-wire 
line) gives the total normalized reactance terminating 
the coaxial line. ) may be obtained for a purely 
reactive termination from (1) with appropriate 
change in subscript. With 2?<1, 


®) = 7/2+-2. (17) 


The formula gives ®) for small 8h where the overall 
electrical length of the antenna is B(A+/) = 7/2. 

When the antenna length is B(h+/)=7, the 
current is small at the driving point and the reactance 
terminating the coaxial line may be assumed to be 
the input reactance of a section of open-end line of 
length (h+/). The equivalent lumped capacitance 
of the open end, Cy,, may be transformed to the 
driving point, so that 2=27,=X7,/R;. In this 
case, since 22>1, 


= r—1/x=7+0l p,. (18) 


In the preceding discussion, approximate correc- 
tions have been made for reactive end effects at the 
load end of the coaxial line and at the open ends of 
the open-wire lines forming the antennas. It is to 
be noted that an accurate analysis of end effects 
when the generator is in the terminal zone is not 
available. In the case at hand the diameter of the 
coaxial line is of the same order of magnitude as the 
smallest heights A of the antennas so that the electro- 
magnetic field in the vicinity of the driven end of 
the entennas and the load end of the coaxial line is 
complicated. An accurate representation in terms 
of an equivalent lumped load for the coaxial line and 
equivalent point generator for the antenna is not 
possible and a close quantitative agreement between 
theoretical results based on such a representation 
and measured values cannot be expected in general. 





3. Experimental Results for Inverted L-, T-, 
and Multi-Element Top-Loaded Antennas 


The measured values of the apparent terminal 
functions O4=Psat+jPsq are shown in figure 2 for 
the inverted L-antenna with 6(h+/)=3.2. The cor- 
responding corrected functions 6)=po+j®) with 
their theoretically derived extensions to small values 
of Bh are also shown. The apparent terminal im- 
pedances Z,, and the corrected impedances Z,— Ry 
+jX, are shown in figure 3 for the inverted 
L-antenna with B(h+/)=3.2. Note that Z,, is 
the impedance actually observed as terminating 
the coaxial line; it includes end effects. Zp, is the 
ideal impedance of the antenna in the absence of 
end effects. Figure 4 shows the impedance of the 
inverted L-antenna, T-antenna, and antennas with 
two-, three-, and four-element top loads with 6(h+1) 
=1.6. 


4. Antennas of Transmission-Line Types 


It has been shown by King, Harrison, and Denton 
[1958] that antennas of transmission-line types may 
be treated by transmission line theory corrected for 
radiation if the height h of the antenna is sufficiently 
small. A general formula for the input resistance 
at the driving point is given by King, Harrison, and 
Denton [1958]. This formula does not take account 
of the effect on the driving-point resistance of the 
conductor of length 2A in series with the generator. 
Nor is this effect included if the inductive reactance 
X,=oL of the section of length 2A is added to the 
input reactances X1;, and X>;, of the parallel com- 




















0.11 T T T T T T T 
bs e 
0.10 (=== Blh+e)=m 4 
0.09;- 
0.08F- 43.26 
Approx. 
theoretical 
0.07-- 43.24 
p 
0.06F 3.22 
® 
0.05—- 43.20 
(measured) 
0.04F, 43.18 
! Approx. 4 
ed theoretical —_ 
o002Fr 43.14 
0.0IF i 4312 
‘ 
‘7 
- L 1 Ll l 
l lL 
0 4 r) i2 16 20 24 268 3.2 


Fiaure 2. Terminal functions for inverted L-antenna, B(h+1) 


Ge 


451 








B(h+?)=m7 












Xsq (measured) 


U 


‘ 
Pa diea theoretical X 


1 L l 4 4 1 1 
fe) 4 8 12 16 20 24 2.8 
Bh 











3.2 


| since the impedence is simply # 


binations of the two transmission-line sections in 
the form 


XiinXoin 


Xiint Xan — 


rr 
X= Art+X aed in 


Zin=Rint+j Xin Where 
Rj, is assumed to be given by King, Harrison, and 
Denton [1958]. An error due to this effect may be 
expected whenever X,, becomes large and the current 
in the generator is not near a maximum. The 
inclusion of a series inductance X;, in (19) 
meaningful only when the current in the generator 
is near a Maximum owing to the fact that there is 
actually no point generator but rather a distribution 
of driving field at the end of the coaxial line that 
extends over the full length A. For purposes of 
comparison X;, is calculated both with and with- 
out X;,, and curves of both X;, and Xj, are shown. 
The shunt driven transmission line with one end 
open and the other terminated in a straight con- 
ductor has three variables s;, s. and h. Figure 5 
shows the measured impedances when @s,=7/2 and 


is itself 





Bh=0.161 and 0.225. 


In 


this case the left-hand 


section is terminated in a straight conductor and the 


corresponding 


terminal function is given by King 
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Impedances of inverted L-antenna, 8(h+1)=3.2 
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Figure 4. Impedances of inverted L-, T-, and mutti-element 
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top-loaded antennas, B(h+1)=1.6. 





[1955, pp. 123-123]. The right-hand section is 
adjusted in length to be near a quarter wavelength 


long. 


Let Bs.=2/2— 6’ where 


6’ is a small length. 


The right end is open but it is not an ideal open end. 
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Impedances of transmission-line antennas, BS2=7/2. 
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The terminal or end effect is equivalent to a lumped 
capacitance C;, given approximately by eq (16). 
Tbe corresponding normalized susceptance with 
Ri=120 In (2h/a;) is bp. = RiwCz,. The corresponding 
terminal function may be obtained from King 
(1955, p. 119, eq 8]. The theoretical resistance 
curves are shown in figure 5. It is seen that there 
is good agreement between the theoretical and 
experimental values. 

The inductance of a straight section of wire of 
radius a, and length 2h is given by King [1945]. 
The input reactances of the sections of lossless line 
are 


Xi in= — RP, cot (88; +p) ; Xoin= — RR, cot (882+). 
(20) 


For the two cases Bh=0.161 and 0.225, Xz, Xin and 
X2in have been calculated and X;, determined. 
For the antenna over the image plane, the input 
reactance is X;,/2. The theoretical values are 

22.5 ohms for B6’=0.0582| ¢.. 27 a 1¢ 

[ar ohms for 86’=0 J for Bh=0.161 
Awle= 
35.5 ohms for 86’=0.0735| ¢ 97 9 99% 
48.2 ohms for 86’=0 J for Bh=0.225. 





(21) 
The experimentally observed values are 


22.3 ohms for Bh=0.161 
Xin 2= 
37.0 ohms for Bh=0.225. (22) 


The value for Bh=0.161 in (22) corresponds closely 
to the theoretical value for 86’=0.0582 in (21); in 
this case Bs.=1.511 or s,=0.242 X. Similarly the 
experimental value for Bh=0.225 agrees well with 
the theoretical value for s,=0.245 \. The input 
reactance as a function of Bs; is evidently very sensi- 
tive to the particular setting of s.. A very small 
change in s, alters X,, significantly. The theoretical 
expressions involve formulas for end effects that 
provide corrections that are comparable with the 
line spacing but not comparable with the radius of 
the conductors. Evidently changes in Bs, of the 
order of magnitude of Ba,=0.025 or small multiples 
thereof are precisely the order of magnitude of the 
difference between the observed and calculated 
values. 

Figure 6 shows the impedance of the same structure 
for Bh=0.161 and 6s,=0.44 when gs, is varied. 
Resonance occurs near 8(s;+s.)=7/2. The theoret- 
ical resistance P;, and reactance X;, are also shown. 
Evidently the length 8(s,;+s2) is a critical quantity 
in the manner discussed in the preceding special case 
with Bs,=7/2. 

The measured and theoretical impedances of the 
sbunt-driven transmission line antenna with both 
ends terminated in straight conductors (m-antenna) 
are shown in figure 7 when Bh=0.225 and Bs.=7/2. 
It is again evident that the setting of Gs, is critical. 
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Impedances of transmission-line antennas, Bs:= 
0.44, Bh=0.161. 
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A more detailed version of this paper appeared as 
Technical Report No. 305, Cruft Laboratory, 
Harvard University. 
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There is experimental evidence that VLF signals propagating from west to east suffer 
less attenuation than from east to west. Earlier work treating the case of nonreciprocal 
a op eos along the magnetic equator is extended in latitude. The nonreciprocity shown 

y ||R|| for highly oblique propagation along the magnetic equator persists when the re- 
flection point moves towards a magnetic pole, but at the pole itself ||R|| is reciprocal. To 
a first approximation 1R1 is reciprocal at all magnetic latitudes. The conversion co- 
efficients ||Ri, 1 R||, are greater for east-to-west propagation than for propagation in the 
opposite direction, except at a magnetic pole where they are equal. 








1. Introduction 







There is now a considerable amount of experimental evidence [Crombie, 1958; Taylor 
1960] showing that VLF propagation is nonreciprocal, since VLF signals received from a west- 
erly quarter appear to suffer less attenuation than those received from the east. A qualitative 
explanation of such behavior in terms of interactions between the longitudinal movements of 
ionospheric electrons in the plane of incidence, and the transverse component of the earth’s 
magnetic field was given by Crombie [1958]. 

Subsequently, Barber and Crombie [1959] calculated the reflection coefficient of a sharply 
bounded ionosphere and showed that this explanation was correct for the special case of waves 
polarized in the plane of incidence, propagating along the magnetic equator. This special 
case simplified the calculations since the magnetic field is horizontal and transverse to the 
plane of incidence. 

It has also been shown [Crombie, 1960] that the nonreciprocity which could be expected 
from the lack of reciprocity exhibited by the reflection coefficient does occur when mode theory 
is used. Wait [1960, 1961] and Wait and Spies [1960], using some coefficients derived from 
the numerical data of Johler [1961], have shown that mode propagation depends on the direc- 
tion of propagation. More recently, Dobrott and Ishimaru [1961] have also concluded that 
VLF propagation along the magnetic equator may be nonreciprocal. 

It is the purpose of this paper to derive formulas for the reflection coefficient of a sharply 
bounded ionosphere in the rather more general case where plane waves are incident from the 
west, or from the east at any magnetic latitude. 


















2. Introductory Theory 






The starting point of this work is a paper by Yabroff [1957a and b] in which, following 
an outline by Bremmer [1949], he derives perfectly general relationships between the fields in 
a sharply bounded ionosphere, and obtains the appropriate boundary equations. Yabroff 
considered the coordinate system shown in figure 1. The vertical direction is x and the earth’s 
magnetic field H, is contained in the zz plane. The angle of incidence of the waves is §;, and 
the magnetic field of the earth is inclined to the vertical at an angle 8,. The incident, reflected 
and transmitted waves are described in the coordinate system 2’,y’,2’ in which for a wave 
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+B; y’ . 


Fiaure 1. Coordinate system. 


z The positive values of 8; and a; are measured clockwise 
+; from the positive z and z axes respectively. 











propagating in the a;,8; direction, 2’ is in the a;,8;+90° direction; y’ is directed in the a;+90°, 
90° direction and 2’ is in the direction @;,8;. The lower boundary of the supposedly uniform 
ionosphere is the ry plane. If z is magnetic north, then y is magnetic east. 

The incident plane wave F;; can be written as 


F; =Re{ Fr exp| j2 ~ (ct—z cos B;—y sin B; sin a;—2z sin B; COs a; |} 


in which w=wave angular frequency, cis the velocity of light and tis time. Since the boundary 
conditions must be satisfied it follows that the wave in the ionosphere F,,, arising from this 
incident wave, must be of the form 


, -W ° ° ° 
F,,=Re { Fm eXP [ : (ct—Dar—y sin B; sin a;—z sin B; cos x) | } 
D may be regarded as the cosine of the (complex) angle of refraction in the ionosphere. 


Maxwell’s equations and the equations of motion of an electron in a magnetic field can be 
written as 








VXE oH : 
XL — Ot 
VXxH=—NeV+e oF > (1) 
ia ‘i ot 
m oe mvV—p (VXH.) 
dt ~ iit ia ii 


in which, —e, m are the charge and mass of an electron, 

V is the electron velocity, 

v is the collision frequency, 

€9, Mo are the permittivity and permeability of free space. 
From these three equations the three components of H and V may be eliminated giving three 
equations relating the components of E. In matrix form these equations are: 





8 s—hz ghr fre % 4 ; 
aman? et arD— 3p: a,D ary h?) : E.. 
yh 9 ¥ s gh 7 
arD+- Re 1 -a-D a at “ E, |=0 (2) 
hrh hi, i —hr . 
| aD+— = ae ara, +3 — hp 1—D?—a7— ra “| } He | 
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where 


s=(1—jZ)/X h=Y/X, hy=h cos B, hr=h sin B, 


X= Ne?/meqw?= (normalized plasma frequency)’. 


Y=y,eH,./wm=normalized gyro frequency. 


Z=v/w=normalized collision frequency. 


a,=sin B; COS a;, dp=sin B; sin a;, a=sin Bj. 


In order that Z be finite the determinant of (2) must be zero. Thus 


where 


b,=8(s? 


bD*+bsD°+b.D?+b,D+b=0 “ 


—h?)—s?+hi 


b;=2azhrhr 
b.=2s[(1—a*)h?— (s—1) (s—a?s—1)]—A?-(1 — az) — (2—a*) hi. 
= —2(1 —a*)ayhyhr 


bo= (s—1) {[(1—a?)s—1]?— (1 —a?)7hi- cos? a;} 


— (1—a?’)[(1—a?)s—1] (A? sin? a; +h2). (3a) 


The four complex values of D obtained from this quartic represent the two upgoing and two 
downgoing characteristic waves which can propagate independently in the medium. Values of 
D with positive real and negative imaginary parts represent upgoing waves. 


2.1. Boundary Equations 


Since the total tangential fields Z,, L,, H,, and H, must be continuous across the boundary, 
the following boundary equations are obtained: 


From E,: 


Ey, cos a;+E,,; cos B; sin a,+EF,,, cosa,—E,,, cos 8; sin a,x=P\E n+ PE» (4a) 
From E£,: 
—E,,sina;+E,,; cos B; cos a;— E,,, sin a;—E,,, cos B; co8aj=E yt» (4b) 


From H,: 


Ey cos B; sin a,—E; cos a;—Ey, cos B; sin a,—Ey, cos a,=SiE a +S,E.. (4c) 
From H,: 
Ecos B; cos a; + E,,; sin a,—E,,, cos B; cosa; +E, sina=TE,+T.E 2 (4d) 


In these equations 


P=EyE-M“{ ns 


Q=EJE=M-{ nd 


595362—61—_—_-4 


e st—hz gh gh , Ach 
1—a— ans | | axa = alte 7D+ #5] [ aud 2 s( aor | } 


(5) 








' hah ; ; jh jh 
D+ ame | E —aj—-D— 5° |+| aur | [« D— FF, | 


(6) 
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where M is given by 


?—h sae. 3 ih ih q 
M=[1—a— 8 [1 -at—D 55 LH aed- 2%, |[ on + 5 | (7) 








and 
R=nH,/E,=ar—a,P (Sa) 
S=nH,/E,=a,Q—D (8b) 
T=nH,/E,=DP—arzQ. (Se) 


In which 
20> ~ Mo/€o- 


2.2. Reflection Coefficient R 


This can be written in matrix form as 


Ry Ry 
=F * ° 9 
. Ky ts | ©) 
The first subscript denotes whether the electric field is parallel (||) or pendicular (1) to the 


plane of incidence, and the second subscript refers in the same way to the electric field of the 
reflected wave. The components of F? are defined as follows: 


R= E/E; (10a). LR = E,,,/Ey. (10¢) 
WR =E£,,,/E., (10b). tR =E,,,/Eyi. (10d) 


The components of F are obtained by solution of the boundary eqs (4) for the appropriate 
pairs of electric fields. 


3. Application to Propagation Along a Line of Magnetic Latitude 


The work of Yabroff given above is perfectly general and expressions for the reflection 
coefficient could be derived for any direction of propagation, and for any orientation of the mag- 
netic field. The labor involved would be excessive, however, and the resulting expressions 
cumbersome. The more practical approach to this general case is to use a computer for the 
solution of the boundary eqs (4). This has been done by Johler and Walters [1960] for propa- 
gation in the magnetic meridian. Johler [1961] has also calculated numerically, reflection coef- 
ficients for propagation at angles of 0°, 60°, 126°, 180°, 240°, 300° to the magnetic meridian, 
omitting the directions 90° and 270°, with which this paper is concerned. 

However, in the case of propagation along a line of magnetic latitude the expressions 
obtained for the reflection coefficient are relatively compact as will be shown below, although even 
in this case the algebra is very tedious. 

For magnetic east-to-west or west-to-east propagation a;=90° or 270°. Thus, a,=0. 
Making this substitution in eqs (4) gives for the boundary equations: 


Ey: E/ ,C—E,',C=P,E,+P.E 2 (lla) 
Hz: E,! +E ,,=T:Eat+T.En (11b) 
Ez: —E,’ ,—E,',=Ent+En (11¢) 
Ay: E,’ ,C—E,’,C=S,E,+S.E 2 (11d) 

where C=cos B;. 
(8b) and (8c) reduce to T=DP—a,Q (12) 
S=—D (13) 


458 


F 
4 


RUSE ETE 


——————— 











while (5), (6), and (7) reduce to 


P=P'M-'=M-'4 j[ 1-a- py |+*2( arD+- =i) ans (14) 


Santee a = yd dhs hr Ss : ht = 
Q=Q'M-'=M {=i (aD— 5; 3)" 2 (1—D— <p) rae (15) 


M=|1-a—— =F |[)- —[Pr— zat ee “yi (16) 


The quantities M, P, Q, S, T, carry the same subscripts as the D which they contain. 


3.1. Solution of Boundary Equations 


The boundary eqs (11) can be solved in the following way for the components of R. \R, 
is taken as an example. From (11a) and (11b) 


; , T,\C—P,+(Ep/Ea)(T2,C—P») 5 
ey ee i BEST Be Se Se seed. Pid. 
R=E,’ JE,’ = T "OFP,A( z -/E.a)\ TC rep) (17) 





From (11¢e) and (11d), on putting E,’,=0 it follows that 
F»/En=—(C+D,)/(C+D,). (18) 


Substitution of (18) in (17) gives 


R,=A-"[(7,C—P,) (C+D) — (T,C—P?2)(C+D,)], (19) 
A=(T,C+P;)(C+D,)—(T,C+P2)(C+D,). (20) 

Similarly 
R,=A-[2C(D,—D,)}, (21) 
1R\=A7"[2C(T,P,—TP))], (22) 
LR, =A" (T,C+P,) (C—D,) —(CT.+P,2)(C—D,)]. (23) 


relatively compact expressions for Ry and , Ry 


Ee ai {3 | (1-2) ) CDD.— spam. pageant —Dp bf =5 


Substitution for 7,, P,, JT, and P2 in these equations gives, after much work, the following 


~1|(0.+) + j(205+") e) ert | (a ) (D+ D.D,) +" Fg) als 
ih =ay" [{: | Dac D.—D.)—D,D,{ 1 iff ton '+D,+D; = mA 
—1]+D)+| 5201+") 42? || (a2 p-1) OH DD) +e a |}. @ 


in which 


a.=a sep 1D» |= {i | (1-2) ead. (C+D.+D>) | 


4“r8P (0(C+D+ D2) + DD; i} ee 7 —1]D+D,) 


A? DD, ‘s 
+59 | (ata) OH DD) + | 6) 
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and Dj=1—a*—1/s, the value of D? when there is no magnetic field present. 


Equation (22) reduces to 


LRy= a7" g4,| (1-2) D+ Ic. 27) 


Equation (21) reduces to 
: s—h? ; 
yRa=Ay"- 20MM, ("), (28) 
Ap, 
When MM, is expanded, it is found to contain terms involving Dj7+D} and DjD}. These are 


respectively equal to —b2/bp and b,/b) in eq (3). Substitution of these in (28) gives, after con- 
siderable reduction, 


eee i3| (1-2) D3 + Jc (29) 


which is identical to , Rj. 


4. Discussion of Reflection Coefficient Formulas 
These formulas have been checked under various limiting conditions to ensure that they 
conform with other published results. 
4.1. No Magnetic Field 
If the magnetic field is zero then Ry and ,R, reduce to 


(D2+<a? )C '"—J), 


a= : ~ (DRR+a?)¢ +Dy 
—Dy 
1hi= aaa 


while ,2,=,R)\=0. These limiting forms were given by Yabroff [1959b]. 


4.2. Horizontal Magnetic Field 


If the magnetic field is horizontal, h,=0, hr=h, and (19) reduces to 


; 1 “mp , Ararl2 
[s2.( 1—1) Arar | 0_(jpp4"2r2) 
[am(1— peg er ao ss) 

s s 
From figure 1, hy is negative and for west-to-east propagation @7 is also negative. In Barber 
and Crombie’s paper [1959] their term equivalent to hyay was taken as negative for west-to- 
east propagation, because of the different coordinate systems used by them. Bearing this 


in mind it is possible to show that (30) is identical with their earlier result. 
If in (25) hz =0 then (25) reduces again to 


R C—D, C-Dy 
ina i= C+D C+D, 


(30) 





R= 


confirming the result of Barber and Crombie [1959]. It is easily seen from (27) and (29) that 
when h,=0 both , Ri) and \R, become zero. 
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4.3. Vertical Incidence, Vertical Magnetic Field 


Here hy, a7, a become zero and C=1. 
Substitution in (24), (25), (27), (29), give 


D,D,—1 
(Di+1) (D2+1 


D,—D, 
(D.+1)(D2+1) 








y= 1 R= y Rk) =)Ri=J 


which are results previously obtained by Bremmer [1949]. 


4.4. Reciprocity 


The expressions (24), (25), (27), and (29) all contain Azar. This is positive for west-to- 
east propagation and negative for east-to-west propagation. Thus, provided that A,r is not 
zero, all the reflection coefficient terms may exhibit nonreciprocity. 


5. Numerical Work 


In order to evaluate the expressions for the reflection coefficients it is necessary first of 
all to find D, and D, from (3). For propagation along a line of magnetic latitude a,=0 and 
(3) becomes a quadratic in D?. Thus 


b,D‘+6,D?+b,=0, 
and 
D¥=[—b,+ yb3}—4b b4] (2b4)-1. 
Substituting a,=0 into eqs (3a) gives 
b.=2s [(1—a?)h? —s(s—1) D3] —hz — (2—a*) hi 
b}—4b b,=h*—2 (1—a?) h*hi.+-4s (s—1) Dghz + (1—a?)*hi, 
bs=s (s*?—h?) — (s*?— hj). 


The values of D,, D. having positive real and negative imaginary parts are chosen since these 
1) 5 So fo) “ : ad 
represent the upgoing “ordinary” and “extraordinary” waves in the ionosphere. 
Some special cases are of interest in obtaining the results of section 4. 


For h=0; D?=D3}=D?=1—a?—1/s. 


2 9 2 s—1 
For h,=0; hr=h and i= 0; Py ae 
. 1 
For hr=0; hp=h and a—U, a alana 


Some approximate calculations have been made for daytime conditions to show the effect of 
the angle of dip of the earth’s magnetic field on the four reflection coefficient terms using the 
formulas (24), (25), (27), and (29) developed above. The calculations are based on the fol- 


lowing parameters: 
w= 10° corresponding to an electron density ~310/ce, 
v=10"/sec, w=10°/sec, a; ~78°80’, C=0.2. 
The earth’s magnetic field has been taken as 0.5 gauss and the inclination of the field to the 
vertical B,=90°, 60°, 30°, and 0°. The results of these calculations are shown in table 1, 
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TABLE 1 
























































| 
| Teal 1P1 (Fa =1)\* 
West to East to West to East to West to East to 

| east | west east west east west 
90°._..--.-----|  0.72/175° | 0. 38/175° || 0. 76/171° | 0. 76/171° || 0.0 0.0 | Horizontal field. 
60°_...........| .69/176° |. 38/176° . 76/162° |. 76/162° |} .042/136° | .081/192° |_......-.--.-.------ 
peas --|  .67/179° | —_. 42/188° Fg be de ee a 
P=. ----| .59/180° |. 59/180° ||. 66/150° . 66/150° || .083/153° |. 083/153° | Vertical field. 





*For propagation in the Northern Hemisphere the phase angle of Ry and Fi should be reduced by 180°. 











It is clearly seen from this table that Ry, )R, and ,R) all show nonreciprocal properties, 
whereas ,R, does not. In the case of Ry the difference between east-to-west and west-to-east 
propagation is greatest at the magnetic equator where the magnetic field is horizontal, and 
gradually decreases as the reflection point moves nearer a magnetic pole. Nevertheless, the 
difference is very appreciable even when the angle of dip is 60° (8,=30°). At the pole, where 
the magnetic field is vertical the two reflection coefficients are equal, since there is then no 
transverse field. 

At the magnetic equator )R, and , Ry are zero. When the magnetic field is inclined ,R, and 
1A become finite, but again their values for east-to-west and west-to-east propagation differ. 
This difference again disappears at the magnetic pole. 

Despite the fact that eq (25) for , R, contains a7hz, it seems from table 1 that, for the values 
of the parameters used, ,R, is essentially reciprocal. Thus in (25) the variations due to the 
sign of a7hy must, to a first approximation cancel. This presumably will not necessarily occur 
for other values of the parameters. 

The above numerical calculations have been made for highly oblique incidence, since long 
distance VLF propagation is of interest at present. However, the formulas are applicable at 
any angle of incidence. Thus the predictions of the formulas could be tested by making obser- 
vations of Ry and/or ,R, at points situated at equal distances to the west and east of a VLF 
transmitter. Alternatively, observation of the ground-interference pattern to the west and 
east of the transmitter could be made with airborne equipment. 


6. Conclusions 


This paper extends the earlier work of Barber and Crombie [1959] to the case of west-to- 
east and east-to-west propagation at any magnetic latitude. The nonreciprocity shown by 
2, for highly oblique propagation along the magnetic equator persists when the reflection 
point moves towards a magnetic pole, but at the pole itself Ry is reciprocal. To a first approx- 
imation at least, ,R, is reciprocal at all magnetic latitudes. However, the conversion coef- 
ficients ,R,, ,R\, which are shown to be equal, are greater for east-to-west propagation than for 
propagation in the opposite direction, except at a magnetic pole where they are equal. 


The work described in this paper was done while the writer was on a short visit to the 
Boulder Laboratories of the National Bureau of Standards. The writer wishes to express his 
appreciation for this opportunity, particularly to Dr. J. R. Wait through whom it was arranged. 
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When noise radiation of roughly one or a few kilocycles per second is emitted in the 
higher atmosphere, part of it (an extraordinary wave) is propagated downward into the 
space between earth and ionosphere. Reflection at the earth and ionosphere leads then to 
a standing wave in this space, whose intensity for a given incident power flux varies very 
much with frequency. Maximum field strength of the standing wave is derived for the 
resonance frequencies of the space. The incident wave fronts are assumed to be horizontal. 
Only clearly defined wave fronts being of a sufficiently wide extension and showing no 
noticeable irregularities will lead to a definite resonance. The space between earth and 
ionosphere is comparable to an air gap between two parallel plane reflectors. The lower 
ionosphere is however only a partial reflector, allowing radiation to enter the resonance 
space and causing at the same time some loss of energy out of the resonance space (leakage). 
Stratification of the lower ionosphere has some influence on the resonance phenomena. 
With a D layer, an additional resonance at a lower frequency is obtained. 

It seems too early to decide whether in any observed noises the resonance spectrum of 
the space between earth and ionosphere becomes apparent, but it is expected that noise 
spectra observed on the ground are modified by the resonances. In case of monochromatic 
emissions, the received intensity depends on the position in the resonance spectrum. The 
resonance spectrum should be received in case of emission of a white-noise spectrum, provided 


Resonance of the Space Between Earth and Ionosphere 





the wave fronts are appropriate. 
1. Introduction 
1.1. Cavity Between Two Plates 


In a space between two infinite parallel plates 
resonance is observed when half a wavelength or 
multiples of a half-wavelength fit into the space. 
The electromagnetic field in this case of resonance is 
that of a standing wave whose wave planes are 
parallel to the plates. If one of the plates is made 
partially transparent, a model of the earth and 
ionosphere is obtained. The semitransparent plate 
like the ionosphere lets waves enter from outside—or 
from the higher atmosphere—and_ still permits 
resonance oscillations. The phenomena encountered 
in the case of two semitransparent plates (wire 
gratings) or one conductor together with a semi- 
transparent plate (grating or dielectric plate) have 
been described in the literature [Casey and Lewis, 
1952; Wait, 1954; von Trentini, 1956]. 

The field intensity in the resonance space is 
strongly dependent on the frequency of the incident 
waves, showing peaks at the resonance frequencies. 
The partial transparency of one of the reflectors 
entails imperfect reflection and contributes thus to 
a limitation of the resonance peaks. 


1.2. Fundamentals 


The ground is a good reflector for frequencies of 
the order of 1 ke/s. With a conductivity of 107% 
mho/m, for example, the amount of the complex 
refractive index of the ground for 1 ke/s is |n|=134, 
a value sufficient for good (almost perfect) reflection 
in vertical incidence, which is considered exclusively. 
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The ionosphere is only a partial reflector. The 
extraordinary wave at very (or extremely) low fre- 
quencies finds in the ionosphere (except its lowest 
part) conditions of free propagation, as is known 
from the audiofrequency whistlers. The refractive 
index of the extraordinary wave under these condi- 
tions is in general high (1). When the increase 
to these high values is fairly rapid in terms of wave- 
lengths, we may expect that part of the extraordinary 
wave passes through the height interval of this rapid 
increase and part is reflected there. Thus we have 
for extraordinary waves the combination of an almost 
perfect reflector (the ground) with a semitransparent 
reflector that allows waves to enter the resonance 
space. 

The descending wave fronts and consequently also 
the wave fronts in the resonance space will be as- 
sumed as horizontal. An eventual slight inclination 
of the wave fronts in the ionosphere would, because 
of the refraction in leaving the ionosphere, cause a 
considerable inclination of the wave fronts in the 
resonance space with the result of an essential 
modification of the resonance spectrum. Random 
distortion of the wave fronts would correspond to a 
mixture of propagation directions and consequently 
lead to a smearing-out of the resonance spectrum. 
The lateral extension of the wave fronts must be 
wide enough to have the simple plane-wave theory 
applicable. 

There is some relationship between the resonance 
effect under consideration and the “tweeks’’, se- 
quences of pulses caused by multiple reflection of 
sferics in propagation to a distant location [Potter, 
1951]. The lightning flash from which a sferic 








originates is nearly a point source. The radiation 
thus is omnidirectional and the tweek spectrum 
may not be interpreted as a resonance spectrum 
corresponding to horizontal (or otherwise simply 
defined) wave fronts. The asymptotic limit fre- 
quency of the tweek spectrum, however, is supposed 
to relate to horizontal wave fronts and should conse- 
quently be identical with our fundamental resonance 
frequency, at least in case of a well-defined reflection 
level in the ionosphere. In other respects, the 
behavior of tweeks and that of our resonance 
phenomenon are not expected to be the same. | 

The energy of the resonance oscillations is limited 
by losses. In our case the losses are due to partial 
transparency of the lower ionosphere (leakage), 
absorption on the path in the ionosphere below the 
reflection level, and imperfect ground reflection. 
Only the partial transparency is included in the 
theoretical treatment to follow. At the rather low 
frequencies the partial transparency is thought of as 
being the main cause of loss under most conditions. 
It was noted before that the ground is an almost 
perfect reflector. An estimate of the absorption 
in ionospheric propagation is made in appendix 1. 
The absorption coefficient found for normal daytime 
conditions may not prevent the resonance oscilla- 
tions, but we have to be aware that absorption tends 
to depress the quality factor of the resonances and 
to smooth the resonance curves. This may prefer- 
ably happen in daytime and at the higher resonance 
frequencies. 

The refractive index of the ionosphere for our 
frequencies is in most heights fairly well represented 
by the “quasi-longitudinal” approximation (eq (14) 
of appendix 1, slightly simplified) 


v=1F a (1) 
ffu| cos a| 
The two signs refer to the ordinary (—) and the 
extraordinary (+) wave.’ fis the plasma frequency, 
f, the gyrofrequency of the electrons in the ter- 
restrial magnetic field, and @ the angle between the 
wave normal and the direction of the terrestrial 
magnetic field. Equation (1), sometimes further 
simplified by neglecting unity, has frequently been 
used to describe whistler propagation [Storey, 1953], 
(more references given by Poeverlein [1959a]). The 
equation is based on the assumptions: /o>f, f</u, 
a not near 2/2, and negligible collision frequency. 
In the lowest ionosphere there is, of course, a height 
interval in which fy is comparable to f. At our 
frequencies this height interval is found in a rather 
low altitude and is obviously narrow in terms of 
wavelengths (or decay distances). It has therefore 
no remarkable influence on the propagation and eq 
(1) may be used throughout. At higher frequencies, 
at which f) is comparable to f in a significant height 
interval, the refractive index formula has to be more 
rigorous and coupling or transition between ordinary 
1 This choice of signs complies with Storey’s [1953] classical paper on whistlers 


and with the conventional definition of ordinary and extraordinary waves on 
basis of the Appleton-Hartree formula. The opposite notation, however, is 


occasionally found. 








and extraordinary waves has to be taken into con- 
sideration. This is known to be the case at fre- 
quencies of the order of 10 or 20 ke/s. 

Contributions of positive ions to the refractive 
index become noticeable near the gyrofrequency of 
the ions (approximately 50 c/s for 0+). All deriva- 
tions from eq (1) are therefore useful only for 
frequencies well above the ion gyrofrequency. 


1.3. Objective 


This theoretical study is concerned with the modi- 
fication of noise signals originating in the higher 
atmosphere by resonance of the space between earth 
and ionosphere. Frequency response curves of this 
space will be derived, which show for the frequencies 
under consideration the amplification or reduction 
of the signal intensity on the ground compared with 
the intensity of the incident wave. 

First some simplified models of the ionosphere are 
assumed, exhibiting one or two reflecting boundaries 
of homogeneous media. Later, samples of nighttime 
and daytime stratifications of the lower ionosphere 
are considered. With a simple reflecting boundary, 
resonance is obviously obtained when multiples of 
half-wavelengths fit into the resonance space. 

The study applies to wide horizontal wave fronts 
only. A noise emission coherent over a wide area is 
therefore required. Such an emission may arise from 
large-scale geophysical processes (e.g., incident cor- 
puscular radiation) rather than from local events. 

White noise or short pulses of a sufficient coherence 
over a wide area would lead to an immediate observa- 
tion of the derived frequency response curves. Mono- 
chromatic noise undergoes an enhancement or 
reduction of intensity depending on the position of 
the frequency in the frequency spectrum (response 
curve). Generally the noise spectrum observable on 
the ground is the result of the noise emission, eventual 
propagation effects on the path through the iono- 
sphere, and the response of the resonance space. 

Processes of noise generation in the higher at- 
mosphere have been treated to some extent in the 
literature [Gallet and Helliwell, 1959] (some more 
references are given by Gustafsson, Egeland, and 
Aarons [1960]). A survey of any amount of observa- 
tional material as for indications of our resonance 
effect has not yet been made. 

This study does not show much of a relationship 
with familiar theories of ELF or VLF propagation. 
This is mainly due to the consideration of waves 
that come down from high altitudes and travel ver- 
tically. Other investigations are concerned with 
propagation from a source to a receiver, both between 
earth and ionosphere. Some remarks on absorption 
and on a wave theory of propagation that might be 
of more general use appear in the appendixes. 


2. Simplified Models 
2.1. Earth and Sharply Bounded Ionosphere 


_ The earth will be considered as a perfect reflector 
in all cases that will be discussed. The simplest 
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model of the ionosphere is a homogeneous medium 
with a sharp boundary at the bottom. The refractive 
index of the ionosphere according to eq (1) has a high 
positive value for the extraordinary wave and in 
general an imaginary value for the ordinary wave. 
The extraordinary wave will be dealt with now. 

The simple model of the resonance space and the 
ionosphere is shown in figure 1. The electric field in 
the resonance space is drawn in for various resonance 
cases. The field may originate from some noise 
emission in the ionosphere. A descending wave 
generated in the ionosphere will reach the bottom 
of the ionosphere and lead to a standing wave in the 
resonance space below and to a reflected wave trav- 
eling back within the ionosphere. The standing wave 
in the resonance space may be interpreted as the 
superposition of a descending and an ascending wave 
traveling in a vertical direction. 

Ascending and descending waves will be denoted 
by F and G respectively. Strictly speaking, the 
symbols will denote the horizontal electric field 
strength components of these waves. In the aniso- 
tropic ionosphere an additional vertical (longitudinal) 
field strength component may exist. 

Figure 2 is a representation of the horizontal 
electric field strengths F and G in the complex plane. 
The vectors indicate amplitudes and phases of 
ascending waves (F) and descending waves (G@). 
The polarization of the waves (disregarding vertical 
E components in the ionosphere) is assumed constant 
over the entire propagation path. In case of a 
reflecting artificial dielectric, for example, the polar- 
ization may be linear in a fixed direction. In case 
of the extraordinary wave which is determined by 
the ionospheric propagation characteristics, the po- 
larization is circular in the quasi-longitudinal 
approximation that is being used now [Ratcliffe, 
1959]. The sense of the rotation is clockwise for an 
observer looking in the direction of the terrestrial 
magnetic field. 

The construction of figure 2 starts at the earth’s 
surface. The field strengths F, and G, there are 
opposite to each other in consequence of the assumed 
perfect reflection. In traveling upward the two 
vectors rotate in opposite senses. The rotation in 
the air space amounts to kh and is, with a time factor 
exp (iwt), clockwise for the ascending wave (fF). In 
traversing the bottom of the ionosphere there is a 
sudden change from F, and G, to F; and G3. 
The total horizontal field strength components 


Enor= F+G 


, ao V2 n(F—@) 
Ho 

must be continuous at the boundary. ,,, there is 
imaginary, Ayo, real (see fig. 2). In passing the 
boundary, therefore, the end points of the vectors 
in figure 2 are shifted symmetrically along a hori- 
zontal line with a reduction of the real compo- 
nents by a factor equal to the refractive index 
ratio (n of the ionosphere in our case). It may be 
recalled that of the extraordinary wave is assumed 
to be real and greater than unity. 


(2) 
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FiagurE 1. Earth and sharply bounded ionosphere. 


The broken lines represent the distribution of the electric field strength in 
various resonance cases. The numbers denote the levels for which field strength 
vectors are shown in figure 2. 


It is superfluous to consider the longitudinal elec- 
tric field strength components of the waves at the 
boundary. Continuity of a resulting normal D com- 
ponent at any boundary is guaranteed by requiring 
continuity of yor. 

From our construction it is easily seen that the 
reduction of the field vectors F and G in going up 
through the boundary is maximum if the vectors at 
the boundary are real or 


khang, Qe, Se, . . . . (3) 


This is the resonance condition. In all cases given 
by eq (3), a relatively strong field in the space 
between earth and ionosphere (7; and G;) arises 
from a weak wave incident from above (G3). 

Because G; is the field of the primary wave, one 
might like to start the computation or construction 
of the field strengths with G;. This would however 
necessitate the consideration of multiple reflections 
between earth and ionosphere, which did not become 
apparent in the method chosen. 

The procedure of figure 2 is also usable for more 
complicated models of the ionosphere. A multiple 
layer with several boundaries requires an iterative 
application of the procedure (section 2.2 deals with 
a double layer). A step approximation allows us 
to transfer the method to a continuously varying 
medium (section 3 and appendix 2). 

Figure 2, of course, is nothing but a graphical 
representation of the analytical expressions describ- 
ing the waves in the present case. From the figure 
or from the analytical formulation we derive the 
ratio between the field strengths arriving at the 
ground (G,) and incident on the boundary of the 
ionosphere from above (G3) 


f IMAGINARY 











Ficure 2. Electric field strengths in the complex plane (F 
electric field strength of the ascending wave, G electric field 
strength of the descending wave). 
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The ratio of the Poynting vectors may be con- 
sidered as more significant than the ratio of the 
field strengths. The horizontal field strength com- 
ponent G of the descending wave in a medium of 
refractive index n corresponds to a vertical com- 
ponent of the Poynting vector proportional to 
nG (cf. eqs (2)). An eventual horizontal component of 
the Poynting vector is not of interest in this study. 
The “ratio of the Poynting vectors” in the following 
always refers to the vertical components of the 
Poynting vectors which are proportional to nG’. 
The ratio of the Poynting vectors corresponding to 
the field strength ratio of eq (4) thus is 





1 |G)? _ os : (5) 
n |G3| n® sin?kh+cos* kh 


Equations (4) and (5) show that resonance is obtained 
under the condition of eq (3), with an enhancement 
of the Poynting vector 


Li 
n 


G, 
Gs 


2 


— (6) 








Marimum reduction of the Poynting vector is ob- 
tained for 

"ae. Wee . 

) 


kh=5 575%) >> (7 


with a reduction ratio 


16 
n |G 


ae! 
- 


(8) 





For comparison the results of the Fresnel reflection 
theory for waves traversing a boundary may be 
noted: If the waves come from the medium with 
high-refractive index, a fraction 2n/(n+1) of the 
electric field strength or a fraction 4n/(n+1)? of the 
Poynting vector appears in the second medium. 
Equation (6) shows that the resonance effect in our 
case seems to improve the transparency of the 
boundary (in terms of Poynting vectors) by a factor 
of approximately n?/4. The reduction in the 
minimum case amounts to approximately 1/4 (eq 8). 

The enhancement of the Poynting vector in the 
resonance case does of course not mean a violation 
of the law of conservation of energy. With perfectly 
reflecting earth, there is no net power flux in the space 
between earth and ionosphere. The Poynting vector 
considered in this space is the Poynting vector of 
one of the two progressive waves of which the 
standing wave is supposed to be composed. 


2.2. Earth and Double-Layer Ionosphere 


The single-layer model (fig. 1) may be a fairly 
good approximation of the lower ionosphere at night. 
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Ficure 3. Earth and double-layer ionosphere. 





The daytime ionosphere in lower heights is better 
approached by two adjacent dielectrics, representing 
the D and the £ layers.?, This model is shown in 
figure 3. The heights are chosen in accordance with 
normal daytime conditions. The refractive indices 
for the extraordinary wave at a frequency around 
500 c/s can be assumed as 


Np=8 
il }- (9) 
Nep=50 

The vector diagram, figure 2, must now be com- 
plemented by a rotation of the vectors corresponding 
to the phase shift in traveling up through the D 
layer and one more horizontal shift of the vector end 
points corresponding to the transition from the D to 
the E layer. In passing through the boundary 
between the D and £ layers, the horizontal com- 
ponents of the field strength vectors in the vector 
diagram are reduced by the refractive index ratio 
Ne/Np. 

In the present case two kinds of resonance may be 
Matiowtiead, the one resulting from maximum 
reduction of the electric field strengths F and @ 
(fig. 2) at the bottom of the D layer, the other 
from maximum reduction at the boundary between 
D and E layers. The first kind of resonance requires 
that the field vectors F and @ in reaching the D 
layer are real, or in other words, that an integral 
number of half-wavelengths fits into the space 
between ground and D layer. With a height of 70 
km for the bottom of the D layer (fig. 3) this 
yields the resonance frequencies 


f=2.14XK10°xm ast 


10 
m=1,2,3, . _ 


The minimums between the resonances are found at 
26 eee 
M=3, 3) 3; - . a : 
The second kind of resonance, corresponding to 
maximum reduction at the boundary between the 
D and the E layer requires real field vectors at this 
boundary. With the model of figure 3 and the 
refractive indices of eqs (9) resonance frequencies of 
509, 1386, . c/s and minimums at 151, 945, . 
c/s were derived. 
2 A brief account of a resonance case in the two-layer model was given in an 
earlier paper [Poeverlein, 1959b]. 
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A simple approximation for the second kind of 
resonance is obtained by assuming np>1. With 
large np the horizontal components of F and G (fig. 
2) in passing the bottom of the D layer are reduced 
to almost zero. This means that F and G@ imme- 
diately above the bottom of the D layer become 
imaginary. A phase shift of 3, $7, $7, . 
within the D layer makes then the vectors F and @ 
real at the top of the D layer. The corresponding 
resonance frequencies under our conditions are 469, 
1406, 2340,. . . e/s 


3. Realistic Cases 


3.1. E Layer, Extraordinary Wave 


A nighttime F layer profile given by Nertney 
({[Poeverlein, 1959a] fig. 18) may now be used for the 
study of resonance of the extraordinary wave. The 
profile and an approximative step function are 
shown in figure 4. 

The refractive index in the ionosphere (eq 1) 
depends essentially on the plasma frequency fp. 
Values of fy following from 

f§=80.7X10°XN (11) 
are noted in figure 4 (fy in c/s, the electron concen- 
tration N per cm*). In the refractive index formula, 


fa |cos a| =1.316X 10° c/s (12) 


may be assumed (corresponding to some temperate 
latitude) and hence 


f3/( fz |cos a|) =61.3N. (13) 


The step function of figure 4 similarly to a combi- 
nation of several dielectrics (sections 2.1 and 2.2) 
allows the iterative application of the graphical 
method (fig. 2). The step-by-step construction of 
the field vectors was begun at the perfectly reflecting 
ground and conducted up to 104 km height. It 
appeared that not much contribution to the reflec- 
tion of the waves comes from heights beyond 104 
km, though some inaccuracy may arise from stopping 
at this somewhat low level and from the large width 
of the individual steps. 
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Fieure 4. E layer profile. 
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Fiaure 5. Ratio of the Poynting vectors on the ground and in 


the ionosphere for the profile of figure 4. 


The graphical construction was used to derive the 
intensity variation of the descending wave. The 
ratio of the Poynting vectors on the eround and in 
104 km height, (1/mios) X | Go/ Gros] * , is plotted versus 
frequency in figure 5. The positions of the maxi- 
mums and minimums are not much different from 
those for a sharply bounded ionosphere (section 2.1), 
in 90 or 100 km height. The gradient of the electron 
concentration apparently is too steep to have much 
influence. A shallower gradient doubtless has some 
effect and must lead to diminished resonance fre- 
quencies as does a D layer below the £ reflection 
level (next section). The greatest enhancement of 
the Poynting vector (fig. 5) is seen to amount to 5 
(at 1.63 ke/s), the strongest |reduction to 0.11 (at 
0.60 ke/s). 
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Ficure 6. Profile of D layer and E layer. 





3.2. D Layer and E Layer, Extraordinary Wave 


In daytime, with presence of the D layer, the 
profile of the ionosphere may according to Sigeeg 
be that of figure 6 ([Poeverlein, 1959a], fig. 17, 
curve (b)). The step approximation of : Be profil 
leads to the Poynting vector ratio, (1/Mr00) X | | Go/Groo | ?, 
represented in figure 7. The maximum enhance- 
ment is 1.02, the strongest reduction 0.06. 

The sequence of the resonance peaks seems to be 
irregular. In fact, the stratification of the iono- 
sphere resembles the simplified model of section 2.2 
for which two series of resonances were found. The 
first resonance peak in the present case (at 440 ¢/s) 
corresponds to the first peak of the simplified model 
(at 509 or 469 c/s) which was explained by the be- 
havior of the field at the bottom of the E layer and 
the phase shift in passing through the D layer. In 
the further course of the resonance curve (fig. 7) it 
seems impossible to separate the two series of reso- 
nance peaks of the simplified model. 

There is no pretension to accuracy of figure 7. 
At the higher frequencies, for which the subdivision 
of the step function (fig. 6) is not fine enough, some 
estimates were used in the values plotted. 


3.3. Ordinary Wave 


For the ordinary wave, which is obtained with the 
minus sign in eq (1), the refractive index is imaginary 
in most parts of the ionosphere. The graphical 
method of figure 2 is still usable for imaginary n; 
however, the ‘electric field strengths F and G become 
purely imaginary and have unequal amounts. The 
two waves (F and @) are attenuated, each in pro- 
gression in its direction of propagation. 
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Ratio of the Poynting vectors on the ground and in 
the ionosphere for the profile of figure 6. 


Ficure 7. 





The resonance case now is the case of disappear- 
ance of G in higher altitudes, strictly speaking in the 


altitude in which the signal is generated. In this 
case the enhancement of the Poynting vector is 
infinite, because no initial Poynting vector is present. 
A finite enhancement would result from considering 
electron collisions in the refractive index formula. 
The present study in which collision losses are dis- 
regarded allows therefore only the computation of 
the resonance frequencies. The first resonance fre- 
quency in the case without D layer (fig. 4) is found 
to be 1.52 ke/s, nearly the same as for a sharply 
bounded ionosphere. 

Whereas resonance of extraordinary waves may 
be excited by waves coming down from any altitude 
within the ionosphere, the ordinary-wave resonance 
requires emission of radiation not too far beyond the 
reflection level, in order to have the ordinary waves 
only slightly attenuated in traveling down. The 
radiation source may be somewhat above the reflec- 
tion level of the ordinary waves or in the part of the 
pry below the reflection level (lower tail of 
E, or D layer). The polarizations of ordinary and 
setesitidind-+r waves show opposite senses of rotation. 


4. Comparison With Observations 


As noted in section 1.3, only an emission of white 
noise or short pulses in the ionosphere, fulfilling the 
coherence requirement, would lead to an observa- 
tion of the derived resonance functions (figs. 5 and 
7). In ease of a given noise frequency the resonance 
function determines the enhancement or diminution 
of the signal on the ground compared with the noise 
intensity descending from higher altitudes. 

The field observable on the ground should be a hori- 
zontal H. On a perfectly reflecting ground there 
should be no horizontal E. Imperfectness of the 
ground reflection, however, makes a horizontal 
possible. Circular polarization as expected from 
extraordinary (or ordinary) waves would correspond 
io two horizontal F (or H) components with 90-deg 
phase shift. 

Noise emissions exhibiting two or sometimes three 
peaks in the spectrum were observed in the auroral 
zone in connection with geomagnetic disturbances 
[Aarons, Gustafsson, and Egeland, 1960; Gustafsson, 
Egeland, and Aarons, 1960]. The observed spectrum 
resembles at least roughly the resonance functions. 
The noise peaks in most cases were near 750 c/s and 
between 2 and 3.5 ke/s. The peak near 750 ¢/s is 
rather stable in frequency. In one event of noise 
emission, for example, it is reported to have varied 
as little as +150 c/s. At other occasions there were, 
however, larger deviations. The observers suggest 
gyroresonance of protons as cause of the 750-c/s peak. 

The second peak in the observed spectrum is more 
variable. It may be that this peak is caused by 
resonance enhancement of practically white noise 
and the first peak by gyroresonance radiation, pos- 
sibly in near coincidence with our first resonance 
peak under the prevailing conditions. A clue to the 
nature of the noises could be obtained from observa- 
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tion of their polarization. In proton gyroradiation 
ordinary waves are emitted. Our standing-wave 
resonance, on the other hand, may preferrably be 
excited by extraordinary waves, as the remarks in 
section 3.3 indicate. Simultaneous noise observations 
on the ground and in the ionosphere might also allow 
a decision about the character of the noises. 

For a discussion of VLF and ELF noises with 
respect to characteristics and possible origin one may 
refer to Gustafsson, Egeland, and Aarons [1960]. 


5. Additional Remarks 


From a pulse-type signal emitted in the ionosphere, 
a series of pulses will be received on the ground. The 
first pulse reaches the ground immediately after 
leaving the ionosphere and traveling down through 
the air space. The others traveled back and forth 
between earth and ionosphere once or several times. 

As previously, only noise signals that have exten- 
sive horizontal wave fronts are considered. The 
spectrum of the received series of pulses is then neces- 
sarily identical with the resonance spectrum obtained 
from white noise emissions in the ionosphere, be- 
cause a pulse (provided it is short) has the same spec- 
trum as white noise. In the simplified model of a 
sharply bounded homogeneous ionosphere with a 
refractive index independent of frequency, there is 
no distortion of the pulse in the ionospheric reflec- 
tion, only a reduction of the field strength by a fac- 
tor (n—1)/(n+1). The geometric series of pulses 
thus obtained has a Fourier spectrum corresponding 
to eq (4) of the above resonance theory. This proves 
the equivalence of the two considerations in this 
special case. 

The question arises whether the whistlers observed 
in the audiofrequency band are also affected by 
resonance of the space between earth and ionosphere. 
Whistlers in coming down from high altitudes may 
have extensive wave fronts, but presumably not in a 
horizontal position. A possible variation of the 
resonance phenomena resulting from this inclination 
is not so much of interest now. The modification of 
the whistler by the resonance effect means nothing 
but a multiple repetition of the whistler signal, just 
the same as in the case of pulses. It will not be 
discussed now whether this repetition may or may 
not be recognizable in recordings. An important 
feature of the resonance spectrum is the first mini- 
mum, which appears in figure 5 at 600 c/s and will 
be slightly shifted towards higher frequencies in case 
of inclined wave fronts. This minimum may con- 
tribute to a low-frequency cutoff of the observable 
whistler spectrum. 

The resonance effect this paper deals with refers 
to standing waves between two parallel reflectors. 
The curvature of the actual reflectors (earth and 
ionosphere) may necessitate some correction of the 
theory at the lowest frequencies.* The lateral ex- 
tension of the field must be wide in terms of wave- 
lengths, but does not have to include the entire 
spherical shell. 


3 For the general theory of VLF propagation between curved earth and iono- 
sphere, refer to Wait [1960]. 
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A resonance effect caused by the spherical shape 
of the earth is the resonance of the zero-order mode 
propagated between earth and ionosphere which results 
from winding-up of the waveguide to a spherical 
shell [Schumann, 1957; Konig, 1959; Wait, 1960). 
The fundamental frequency of this resonance effect 
corresponds to a phase variation of 27 along a path 
around the earth and amounts according to Schu- 
mann [1957] and Kénig [1959] to approximately 
9 ¢/s. 

The fundamental frequency of our resonance effect 
in case of no leakage into the ionosphere is identical 
with the cutoff frequency of the first-order mode be- 
tween earth and ionosphere, because this cutoff re- 
lates to vertically traveling waves. The higher-order 
resonance frequencies similarly represent the cutoff 
frequencies of the higher-order modes in this case. 


Appendix 1. Effect of Collisions 

The quasi-longitudinal approximation of the re- 
fractive index derived from the Appleton-Hartree 
formula with inclusion of the collision frequency v 
of the electrons is * [Ratcliffe, 1959] 


5 (14) 


a 
w(w— iv + wy! COS a} ) 





nr=1 


This approximation is usable at sufficiently high 
frequencies and at very low frequencies, provided 
the angle a between wave normal and _ terrestrial 
magnetic field is not near 90 deg. In both cases 
the approximative formula can be further simplified. 
The high-frequency approximation of n with the 
assumption of w>w, w>wy, and v<w is 


(15) 


At sufficiently low frequencies we may assume 
Wy > w,w2 Pwwy, v<wy and obtain then for the ex- 
traordinary wave 


=( = ) 0-53 —_— ) (16) 
wy |COS a} 2 Wy |COS a) 


The absorption in ionospheric propagation is 
w ; . , 

— Imag (n) - ds), The absorption 
* 


coefficient thus is at high frequencies 








given by exp 


w 1 wv a 
— |Imag (n)|=5— - (li 
c! g ( Fee ) 
and at sufficiently low frequencies 
w o: a? v ? 
— |Imag (n)|== ———_, = 18 
c | g (n)| 2 (wy |cos a )?/? ¢ (18) 


The absorption coefficient for high frequencies, eq 
(17), shows the familiar frequency dependence of 


4 Angular frequencies are written in place of the frequency quantities of eq (1). 





the “nondeviative” absorption (with neglect of wy 
compared with w). The absorption coefficient for 
very low frequencies, however, decreases with de- 
creasing frequency. This has the consequence that 
at our rather low frequencies the absorption in the 
lower ionosphere is only small. 

A numerical evaluation based on eq (14) and 
daytime data of the D layer (fig. 6 and » as noted) 
gave the following amounts of the absorption 
coefficient of the extraordinary wave at 1 ke/s (for 
fy\cos a|=1.2-10® e/s): 


0.021 per km in 72 km with v=10’ per sec, 


0.019 per km in 80 km with »=2-10° per sec. 


These values should lead to a considerable absorp- 
tion on paths of 50 km length, but the absorption 
coefficient decreases rapidly above 80 km and the 
decisive vertical path lengths are shorter. 

A closer investigation of the absorption and its 
influence on the resonance curve are still desirable, 
but the above estimate suggests that the absorption 
is small and may not deteriorate badly the resonance 
curves, particularly not at the lower frequencies. 
At night the absorption coefficients are much smaller 
than the ones derived for daytime and in daytime 
the lowest resonance frequency is below 1 ke/s, as 
figure 7 shows. 

It may be pointed out that the absorption co- 
efficients for high and very low frequencies are 
different not only in their variation with frequency 
but also in their dependence on w. The absorption 
coefficient, of course, is not in an immediate connec- 
tion with the reflection coefficient of the lower 
ionosphere at very low frequencies. 


Appendix 2. Wave Theory Based on the 
Vector Method 


The construction of the complex field strength 
vectors (fig. 2) can be used for a rigorous wave-theo- 
retical treatment of propagation in a medium that 
varies continuously with one coordinate. It was 
seen that the method of constructing field vectors 
is applicable to a step approximation of a continu- 
ously varying medium. The start of the graphical 
procedure was at a perfect reflector, in our case the 
ground. The solution consequently is a standing- 
wave solution. 

If propagation of plane waves in a continuously 
varying medium is to be studied, we replace the 
medium by an approximative step function and add 
a perfect reflector somewhere outside, thus creating 
conditions for the vector method and a standing-wave 
solution. A traveling-wave solution or any more gen- 
eral solution follows from superposition of two in- 
dependent standing-wave solutions. These two 
standing-wave solutions may either be obtained by 
assuming the perfect reflector at two different loca- 
tions or by assuming two alternative reflectors, one 
that does and one that does not reverse the electric 
field strength in reflection of the incident wave. The 
only approximation introduced is the step function. 





By muking the steps differentials the method becomes 


rigorous. Nevertheless, the method has certain 
limitations. Before showi ing them, a brief analytical 
formulation will be given. 

The medium will be assumed to vary with the 
coordinate z. For waves propagated in the direction 
of the gradient, the field strengths may then be 
written 


E,=C(z) sin [} n(z)kodz+¢(2)] 


H=i/2n(20(2) cos [[n(2hete+0(2) 
(19) 


The time variation factor exp (iwt) has been omitted. 
The amplitude function C(z) and the phase function 
¢(z) aresought. EF, and H, represent transverse field 
strengths only. An additional longitudinal electric 
field component EF, may however exist, due to the 
anisotropy of the medium. The polarization of the 
wave in the z,y plane must be constant in order to 
have the continuity considerations of the next para- 
graph valid. The polarization may be linear (with 
E, and H, only) or circular (for example with 
E=<8,, H,=—iH,). 

"The continuity of E, and H, at boundaries requires 
that at the steps of the step function 


pornisd (20) 
A(nCc)=0 


Herein s and ¢ are abbreviations for the sine and 
cosine appearing in eqs (19): 


s=sin [| nkede-+0) 
(21) 
c=cos [| nkude +¢] 


The integral term in the argument of these functions 
varies continuously. Only ¢ is supposed to jump at 
the steps of the step function. The differences of 
s and ¢ needed in eqs (20) are therefore 


As=cAd 
iene 8 2) 


Yquations (20) together with eqs (22) yield 


AC An , 
—~—=——¢ 
‘By n 
(23) 
Ad=— 8¢ 


Making the steps differentials leads to the rigorous 
wave-theoretical differential equations for C and 





qo __dn 
Cn (24) 
ste 





n 
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It has to be noticed that the quantities s and ¢ 
in eqs (24) are functions of ¢ (eqs (21)). The last 
of eqs (24) allows the computation of the phase 
function @ without considering the amplitude funce- 
tion, whereas the first equation involves both the 
amplitude function and the phase function. Equa- 
tions (24) are also readily obtained from Maxwell’s 
equations. 

A set of equations closely related to eqs (24) has 
been derived by van Cittert [1939] [Kofink, 1947]. 
Van Cittert considered oblique incidence on a one- 
dimensionally stratified medium. His equations 
therefore are more general in this respect. The 
medium he assumes as isotropic. The equations are 
necessarily different for the two cases of (linear) 
polarization. 

Slow variation of the refractive index suggests av- 
eraging the factors c? and se over many periods of 
cand s. Equations (24) then become approximately 


dC ldn \ 


do=0 
whence we obtain 
1 
Ca —= 
yn . (26) 
¢@= const 


This is the behavior of amplitude and phase in the 
geometric-optical approximation. 

The above study was limited to normal incidence 
on a one-dimensionally stratified medium. <A general- 
ization for oblique incidence (on a one-dimensionally 
stratified medium) may not be difficult, as van 
Cittert’s formulation shows. The more severe limita- 
tion arises from the demand of constant polarization. 
In case of oblique incidence the polarization in the 
planes of the stratification has to stay constant. 
This is easily achieved in isotropic media. In 
anisotropic media constant polarization is obtained 
in normal incidence if the polarization can be as- 
sumed as circular, independent of the location in the 
layer. In more general cases of propagation in 
anisotropic media (oblique incidence or elliptic 
polarization) the polarization in the planes of the 
stratification varies. The continuity at the steps of 
our step function can then be preserved only by 
simultaneous consideration of ordinary and extra- 
ordinary waves with variable amplitude factors. 
This indicates the well-known phenomenon of cou- 
pling between the two types of waves. 





The step method is not applicable near a zero of the 
refractive index. In approaching the zero, the re- 
quired step length decreases towards zero. However, 
a rigorous expression for the refractive index that 
includes the collision frequency has no zero at any 
location. 

Complex refractive index leaves the analytical 
formulations valid, but makes the graphical pro- 
cedure of figure 2 more intricate. 





The author thanks James R. Wait for a stimulating 
discussion and Thomas J. Birmingham for help in 
preparing the manuscript. 
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Propagation attenuation rates for frequencies below 1 ke/s in the ELF region (30 ¢/s to 
3,000 c/s) were calculated from the spectra of atmospherics observed at widely-spaced sta- 
tions. Data are presented for east to west propagation under sunset approaching the eastern 
station. Under these conditions, the attenuation rates are about 1 db/1,000 km at 75 c/s and 
increase with increasing frequency, attaining about 3 db/1,000 km at 200 c/s. The attenuation 
rates observed seem to be consistent with a two-layered ionosphere model with its lower 


region 90 km above the earth. 


1. Introduction 


It is now accepted that radio waves in the ELF 
(extremely low frequency) region, extending from 
c/s to 3,000 c/s, propagate mainly in the waveguide 
bounded by the earth and ionosphere in the principal 
or zero-order mode. In recent years, a number of 
workers have investigated the propagation of ELF 
waves using atmospherics as signals, in lieu of more 
suitable man-made signals [Holzer and Deal 1956, 
Chapman and Macario 1956, Liebermann 1957, 
Pierce 1960, Smith 1960, and Watt 1960]. 

In most of the work carried out to date, observa- 
tions were made at a single station of a large number 
of atmospherics having their origins in lightning 
discharges at different ranges from the receiver. The 
locations of individual lightning discharges were de- 
termined using direction finding systems or, in some 
cases, estimates of the storm locations were arrived 
at from meteorological information; thus, Chapman 
and Macario [1956] recorded waveforms at a single 
station while the locations of the lightning flashes 
were obtained from an atmospherics direction finding 
network. Propagation attenuation rates were then 
extracted from differences between the spectra of 
waveforms received at different ranges and from 
different sources. In this technique, some assump- 
tions regarding the constancy of the radiation spectra 
of the sources must be made. 

It would be of interest to determine propagation 
attenuation rates from waveforms of the same at- 
mospherics observed at different ranges; especially 
with the observing stations lying approximately 
along a great circle path passing through the origin 
of the lightning discharges. In this favorable cir- 
cumstance, fewer assumptions regarding the radia- 
tion spectra of the sources need be made. 

Such observations were made recently through the 
cooperation of the USAF who supplied atmospheric 
direction finding records from which the locations 





and times of occurrence of individual lightning dis- 
charges were determined. 

In this paper some preliminary results of atmos- 
pheric waveforms in the ELF band are presented. 
Data are presented here resulting from fortunate 
alinements of storms and two recording stations 
approximately along great circle paths. The re- 
sultant attenuation rates are compared with theoreti- 
cal and other experimentally-determined values. 

The simultaneous recordings of atmospheric wave- 
forms of the vertical electric field were carried out at 
Boulder, Colorado; Fairbanks, Alaska; and Maui, 
Hawaii in the spring of 1960. The equipment at 
each of the stations consisted of a 25 ft vertical 
antenna with an isolating cathode-follower con- 
nected directly to its base. The input impedance of 
the cathode-follower was approximately 400 meg. 
A coaxial transmission line connected the output of 
the cathode-follower with the amplifiers located in 
the recording building approximately 200 feet away 
from the antenna. The remainder of the equipment 
consisted of an amplifier, filters, oscilloscope and 
strip-film camera. The recording equipment and its 
rs vam were the same at each station. 

The overall amplitude versus frequency response 
of the recording apparatus, including a capacitor of 
100 puf to simulate the antenna, is flat within 1 db 
between 30 c/s and 400 c/s and within 3 db between 
15 c/s and 1 ke/s (see fig. 1). The decrease in ampli- 
tude response at the higher frequencies, to eliminate 
interference from transmitters operating in the VLF 
region, was accomplished through the use of RC 
filters. The decrease in response at frequencies below 
45 c/s is a result of the time-constant of the interstage 
coupling capacitors and a reduction in voltage trans- 
ferred from the capacitive antenna to the cathode- 
follower. 

All of the receiving sites were located two miles or 
more from commercial power lines to minimize the 
reception of power line interference. Electrical 
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Figure 1. Amplitude versus frequency response of recording 


equipment. 


power was supplied to the recording equipment by 
diesel-electric generators. After taking these pre- 
cautions, it was found that the predominant inter- 
ference was received from the distant power lines. 
At Maui and Boulder, the 60 c/s interference (before 
filtering) amounted to about 5 mv/m and 10 mv/m, 
respectively. By comparison, the amplitude of 
received signals varied over a range of about 5 
mv/m to 100 mv/m. It was found advantageous to 
use a narrow-tuned filter to reject the 60 c/s power 
line interference. In plotting the spectra of the 
atmospherics, spectral components at 60 c/s were 
found by interpolating from components at fre- 
quencies either side of 60 ¢/s. 

Recordings were made simultaneously at each 
station on continuously-moving film, which was 
advanced past the cathode-ray tube at a speed of 
50 cm/sec. The film motion provided the sweep 
necessary to resolve the ELF components. 


2. Theoretical Considerations 


Let G(t,d) represent the amplitude variations 
of an atmospheric waveform in time, ¢, and at 
distance, d, such that 


| G(t,d)=0 for t<0 and t>r, (1) 
where + is the duration of the waveform. The 
Fourier spectrum of the waveform is then 


E(wd)= f "@(td) + e-**tdt= f "G(t,d)e-*"dt, 
0 0 
(2) 


where w is the angular frequency. It is assumed 
that the source is equivalent to a vertical dipole and 
that the ELF waves of interest here propagate in 
the earth ionosphere waveguide in the zero-order 
mode. The spectrum amplitude of such a pulse, 
propagating in the zero-order mode, can be repre- 
sented by an expression of the form 


E(w,d)~A(@)[a sin (d/a)}-/2e-«"4, (3) 


| where A(w) represents the amplitude spectrum of 


the source function, a is the earth’s radius in kilo- 


| meters, and a) is the attenuation in nepers per 
unit distance. 


Equation (3) is valid at ELF pro- 
vided the sources are at least 1,000 km from the 
receiving antennas [Wait, 1960]. It is also necessary, 
of course, that the observer is not near the antipode 
of the source. 

If the waveform of the same atmospheric is 
observed at two distances, d,; and d:, where d,>d,, 


| it follows that 








|E(od,)| sind Jay 
exp (als) a= Fea), Lain dja] 


Using this expression, the attenuation rate as a 
function of frequency was obtained from the spectra 


_of atmospherics observed simultaneously at distances 


d, and d,. Then the attenuation rate in decibels 
per 1,000 km of path length is given by 


sin d,/a]_ _ 10° 
°8 sin d,/a_| dy—dy’ 
(5) 





a(w)=| 20 log eee —101 


where the logarithms are to base 10. 

The amplitude spectra of atmospheric number 1 
observed at the three stations are shown in figure 2. 
The spectrum of highest amplitude represents 
observations made at Boulder, the shortest range; 
the spectrum of smallest amplitude represents 
observations made at Maui, the greatest range; while 
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| Ficure 2. Amplitude spectra of atmospheric number 1 
observed at Boulder, Fairbanks, and Maui. 
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the intermediate values represent the Fairbanks 
observations. It is evident from these spectra that 
the frequency of peak amplitude decreases with 
increasing distance. For example, the frequency of 
maximum energy of atmospheric number 1 is about 
300 c/s at 4,500 km and 100 ¢/s at 9,600 km. It will 
be noted that the spectral width of the atmospheric 
narrows with increasing distance. The range of 
frequencies over which the attenuation can be cal- 
culated is limited to the width of the narrowest 
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Figure 3. Amplitude spectra of atmospheric number 2. 
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Figure 4. Amplitude spectra of atmospheric number 8. 





spectrum. The spectra and bandpass characteristic 
of the receiver indicate that the attenuation rates 
for the Boulder to Maui path should be valid over a 
frequency range of about 50 c/s to about 250 e/s. 

In a similar manner, the spectra of atmospherics 
numbers 2, 3, and 4, observed at Boulder, Fairbanks, 
and Maui, are given in figures 3, 4, and 5. 

In figure 6, the propagation attenuation rates in 
db/1,000 km derived from atmospherics’ numbered 
1, 2, 3, and 4 are plotted against frequency for the 
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FicurE 5. Amplitude spectra of atmospheric number 4. 
10 T T T T T T T T 
BOULDER - MAUI 
8 ; | | 
NY 
= | . 
Sst Sf/* e 
; i / 
be lV: 
= ~ “4 4 
= - 
5 WF |: 
S4b Y, i SFERIC 4 - 
= é dia 
=< fi | : 
© ¢ —e 2 = 
y anes 
2 i: / J CALCULATED 4 
0 i 1 l l 1 1 L i 
0 200 400 600 800 


FREQUENCY, c/s 


FicurE 6. Propagation attenuation rates versus frequency 
derived from atmospherics 1-4 compared with calculated 
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Boulder to Maui pair. Attenuation rates derived 
from numbers 1, 2, and 3 were recorded when sunset 
was at the eastern end of the path (1835 MST) and 
have minimum values of about 1.25 db/1,000 km 
near 75 c/s. Somewhat lower rates are observed at 
75 c/s for atmospheric number 4 which was recorded 
an hour later. It is tempting to attribute the lower 
attenuation rate observed at 1935 to the increased 
amount of darkness along the path—although the 
occurrence of sunset along the path complicates the 
situation. 


3. Presentation and Discussion of Results 


During two, 2-minute recording intervals beginning 
1835 and 1935 MST on March 8, 1960, waveforms 
of five atmospherics were recorded at Boulder, 
Maui, and Fairbanks from storms in the Atlantic 
Ocean at approximately 39° N and 58° W. This 
information is tabulated below. 








Approximate range | 








Atmospheric | from source to Time of reception 
of —— at. 
Pine) oe | Boulder 
No. | Location Boulder Maui | 
| | ee 
1 | 39°N, 59°W | 4500 9600 | 1835:52 MST 
2 | 39°N, 65°W 3500 8700 | 1836:25 
3 | 37°N, 61°W 3500 8800 | 1836:07 
4 | 39°N, 58°W 4500 | 9600 | 1935:58 
| 








The propagation paths from the sources to Boulder 
were in darkness. At the time of the first observation 
(1835 MST), ground sunset was approximately 600 
km west of Boulder and approximately 1,900 km 
west of Boulder during the second observation. 
In the first case, approximately 90 percent of the 
path from Boulder to Maui was in daylight; in the 
second case, approximately 75 percent of the path 
from Boulder to Maui was in daylight. 

The paths of propagation from the lightning dis- 
charges to Boulder and Maui, respectively, are 
approximately east-to-west and transverse to the 
earth’s magnetic field. Waveform observations were 
also made at Fairbanks, Alaska and attenuation 
rates were calculated from the Boulder and Fairbanks 
data also, although this is a less favorable situation 
since the paths of propagation are not along great 
circle paths. 

As a matter of interest, calculated attenuation 
curves are also shown in figure 6. A two layer 
ionospheric model [Wait 1958] is employed; the 
height of the bottom of the reflecting layer is denoted 
h, and the effective conductivity parameter is 
w,, At height d above this, the effective con- 
ductivity is assumed to jump to a new value of w, 
The curves shown in figure 6 correspond to the two 
cases where 


(i) w,=10°, h=90 km, 
curve) 


d= (upper calculated 


line). 





(ii) w =105, A=90 km, d=20 km (lower calculated 
curve) 


w, =10° 


The first of these cases corresponds, of course, to a 
homogeneous ionosphere. The resulting attenuation 
should be proportional to (frequency)?. This would 
lead to a convex curve which is clearly not in accord 
with the data. The second of these curves has the 
required form; the ionosphere is thus behaving as a 
stratified conductor. For the lower frequencies, 
the upper well-conducting layer is becoming effective 
and the attenuation is reduced. Of course, other 
combinations of d and w, could be chosen to fit 
the facts, but it would always be necessary to have 
the effective conductivity increase with height in 
some manner. 

An average of the attenuation rates observed for 
the Boulder to Maui path is shown in figure 7 (solid 
There appears to be minimum in the atten- 
uation rate of about 1.0 db/1,000 km of path length 
near 70 c/s. The rate increases with increasing fre- 
quency attaining values of about 1.5 db/1,000 km 
at 135 c/s, 2 db at 160 c/s, etc. The increase in 
observed attenuation at frequencies below 70 c/s 
is doubtful; it is thought that the data are influ- 
enced by noise at frequencies below about 70 c/s. 
However, it is possible to construct a layer model 
which does have an attenuation increase as the 
frequency becomes very much lower. Such a model 
might consist of a three-layer ionosphere with para- 
meters w,’, w,’’, and w,’’’. If w,’”’ is taken to be 
somewhat less than w,’’ and w,’, a lower frequency 
will be reached where the attenuation is controlled by 
the outermost layer of conductivity w,’’’ [Wait 
1957]. It is also possible that heavy ions may 
influence the results in frequencies below 100 ¢/s. 

Attenuation rates derived from the Boulder 
Fairbanks and Fairbanks to Maui data are also 
plotted in figure 7. These data are subject to ad- 
ditional limitations in interpretations; however, the 
rates are about the same as those presented for the 
largely sea-water Boulder to Maui path. Jn view 
of the theory for ELF propagation [Wait, 1960], the 
effects due to ground conductivity would be 
expected to be small if present. 

The propagation attenuation rates in figure 7 are 
compared with more complete data for day and for 
night by Chapman and Macario in figure 8. The 
present data plots between the daytime and night- 
time values reported by Chapman and Macario. 
Perhaps this is to be expected since the paths were 
partial dark at the east end, and in daylight to the 
west. It is of interest to recall that the Chapman 
and Macario data were recorded at ranges up to 
4,000 km while the present data results from obser- 
vations at a minimum range of about 4,000 km and 
a maximum range of about 10,000 km. Although 
the present data are not very comprehensive, they do 
confirm attenuation rates resulting from other ex- 
periments using different techniques. In addition, 
the feasibility of deriving propagation data at ELF 
from spectral analyses of atmospherics is established. 
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Figure 7. Average of propagation attenuation rates versus 


frequency for atmospherics 1-4. 
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Previous solutions for the ELF mode series are discussed briefly. 
The excitation of the modes for vertical and hori- 


is paid to the height-gain functions. 
zontal dipole excitation is also considered. 


1. Introduction 


In a number of earlier papers [Wait, 1960a and b] 
the author discussed the theory of electromagnetic 
wave propagation in the extremely low frequency 
range (i.e., 1.0 to 3000 c/s). The model adopted is 
a smooth homogeneous earth surrounded by a con- 
centric reflecting layer which is to represent the 
ionosphere. The field was represented by a sum 
of modes for both vertical and horizontal electric 
dipole excitation. Essentially the same model has 
been used by Schumann [1954] but some of his results 
differ from those of Wait [1960a]. It is the purpose 
of this note to discuss certain extensions and limita- 
tions of the theory. 

Initially the earth is taken to be flat. This ide- 
alization simplifies the discussion. Furthermore, for 
ELF, influence of sphericity can be introduced later 
in the analysis. 


2. Formulation 


Choosing the usual cylindrical coordinate system 
(p,6,2) the earth’s surface is at z=0 and the iono- 
spheric reflecting layer is at z=h. Free space of 
electrical constant ¢€,49 occupies the space 0<z<h. 
At z=0 and z=A, it is further assumed that the 
fields satisfy the Leontovich boundary conditions. 
Again this is not an essential restriction but it does 


simplify the discussion somewhat. Thus 
can? at 2=0 (1) 
E,=Z,H, 

and 
H,=Z,H, at z=h (2) 
Eo= —Z.H, 

where 
Z g= (inqw/og)!=A,/n, 

and 


Z = (tuqw/o,)'=Aj/n, 


in terms of the (complex) conductivities o, and o; of 
the earth and the ionosphere, respectively. The 
preceding forms for Z, and Z, are valid at ELF 
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595362—61——_5 


Particular attention 


when the earth and the ionosphere are both equiv- 
alent to isotropic conductors [Wait, 1960b]. In the 
above, 721207 chms and A, and A; are dimension- 
less complex quantities. 


3. General Form of the Mode Series 


First the source is taken to be a vertical electric 
dipole (of current J and length ds) located at z=2p, 
where 0<a<h. From a previous analysis the fields, 
in the space 0<(z<h, are obtained from 


E,=WE), 
= —SE,, 
Hy=—TE,/n, (3a) 
where 
Eo=1(n/d) Ids (e~**) /p, (3b) 


is a reference field and W, S, and 7 are dimensionless 
factors which are expressible as a sum of modes. 
Tn the usual notation, the latter are given explicitly 


by 


Wa—inb et 328,S2H® (kSep)fn(20)fa(2) (A) 


n=0 


Sa et p> 5,S,H® (kSap)fn(20)9n(2) (5) 


r=—¥ eite >» bnSnL{? (kSnp)fn(Z0)fn(2) (6) 











where 
. IR (C)R(C)/07] n 
[1+ UR (ORC) lone, © 
_C—A, 2 
RAOQ=rA (8) 

C—A,; 
Ri(C)=GT a? (9) 
and 

2fn(2)=[R,(C,)] ten +[Ry(C,)te“n?, (10) 
2gn(2)=Crl Re (Cy)] 2e*n*—C,[R,(C,)]'e-*°n?, (11) 








H§ and H? are the Hankel functions of the second 
kind of order zero and one. C, are the solutions of 
the mode equation 


R(C)R,(C)e~" =e" (12) 


and 
S,=(1-C3)}, 


4. Height-Gain Functions 


While the above result was given in previous 
papers [Wait, 1960a and b]itssignificance was not fully 
discussed there. Furthermore, the heights of source 
and observer above the ground were taken to be zero. 
It is of interest to look into the nature of the solution 
under more general conditions. 

The function W is considered first. It is the ratio 
of the actual vertical electric field EF, to the reference 
field E, of the same dipole source if it, and the ob- 
server, were on the surface of a flat perfectly con- 
ducting plane. If kp>>1, corresponding to the 
“far-zone’’, it follows that 


_{2ap = os : 
Way (FF) Si 5, 812-5 F(2)ful2), (13) 
/ n=0 





where } is the free-space wavelength. The vertical 
field Z, is thus proportional to a sum of modes whose 
amplitudes are proportional to 


i exp [—Im 27S, /\]. 
At short distances, such that 
p<h, W- +(1-{ Lk (14) 


provided z=2z=0. This is the expected behavior 
for a dipole if the reflecting layer (at z=h) was neg- 
lected. Numerical values of W were presented on 
revious occasions for the intermediate range where 
p is unrestricted but z=2=0 [Wait and Carter, 
1960]. 

As can be seen from (4), the height dependence of 
the individual modes is described entirely by the 
function f,(z). When z=0, we see that 


C2 


(F(a a (15) 


For most applications to ELF propagation, |Ci|> 
|A7|, and thus f,(0) can be replaced by unity. This 
approximation is contained implicitly in earlier work. 
For present purposes it appears to be preferable to 
define a “height-gain function” in the manner 


fn(2) 


fn(0) (16) 


F,(2)= 




















Using (10) this can be rewritten in the form 


en? + RC, )e7 Hn? 











$= TTC) 
2 
=1+ take Eee) 
+terms in (kz)*, (kz)*, ete. (17) 
Since 
2m » A; tay. Ai 
Comet pa ata = 
it is seen that for the zero mode 
: ik?27A; 
Fy(z)=1+iA,ke— Och Bois (19) 


For most ELF applications /,(z) departs only slightly 
from unity. Also, it should be noted that if the 
ground was a perfect conductor the linear term in 
height is absent. Similar remarks apply to the 
higher modes. In fact, if A,=0, we have 


F,,(2)=cos kC,z 
for all n. 

It is seen from (6) that the height-gain functions 
for the H; component is identical to F,, (z). In 
the case of the E, field it is desirable to define a height 
factor G,,(z) as follows 


Gn(Z) 
Gn(@= FCO)" 


For a given mode, this quantity is proportional to 
the ratio of the horizontal electric field at height z to 





(20) 





the vertical electric field at height zero. From (11), 
it is seen that 
2 C (ei*n?— R eg *,*) 
G,, 2 _Gn( >. 5 n g 
=F) itR; 
‘. kC,,2)? 
=A,+ikCre FO nt... (21) 
For mode zero 
2 .,,k# 
G)(z) =A,—Ag—ddgoet . (22) 


Again, under most practical conditions, only the 


first term of this series need be retained. In fact, 
if z=0, we have exactly that 
G,(0) =A, (23) 


which is simply a consequence of the Leontovich 
boundary condition imposed at z=0. If the ground 
is perfectly conducting 


G,(2)=iC, sin kC,2. (24) 
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5. Excitation of the Modes 


We now consider the excitation coefficient 6, 











defined by (7). Since 
ra) pore oR; 
sphRi= RSA +Riae 
3 2h; Oy 
Rig rs sO meee 
it readily follows that 
1 
b= —— (26) 
n 1A; \ 
lepotat 
kh(C2—A,)? | kh(C2—A?) 


If |A, |<|Ail, which is the usual case, the second 
term in the denominator can be neglected. Then 
to the same approximation 


ice 1 
mee 1 _ sin 2khC, (27) 
+ ORR, 
where use has been made of R;=e?**". Actually, 


(27) is exact if A,=0 and it can be regarded as an 
excellent ‘approximation provided |A,|<|A;|. For 
practical applications at ELF, this may be further 
approximated to 


ies since khCQ)<1 


and ; 
6,1 since 2khC,, ~2nr. 


6. Further Extensions 


The extension to a vertical magnetic dipole source 
is very simple. In view of the symmetrical prop- 
erties of Maxwell’s equations we can write down 


the corresponding field expressions. Thus 
H,=W*H, (28) 
H,=—S*H) (29) 
Ey=T*Hon (30) 


where W*, S* and 7* are dimensionless factors and 


Hy=1(1/nd) Kds(e~™*) /p (31) 
where K is the magnetic current in the vertical 
magnetic dipole of length ds. The explicit expres- 
sions for W*, S*, and 7* are identic ‘al in form to 
those given by W, S, and T (i.e., (4) to (6)). Now, 
however, the ‘modes are solutions of the equation 


RYC)RY(C)e~ 2 =e- 8 (32) 





where m is an integer and 


—Ah ; 
RO= G75 (33) 
and 
— Ah 
RMC)= a (34) 


1 m- : F 
where ater and Ata: The superscript h is to 


denote the horizontally polarized nature of these 
reflection coefficients. 
Jn this case we note that 


_Cn 
C2— (a 





(fa (OP= iC, (85) 
which is a very small quantity compared with unity. 
Consequently, the excitation of ELF waves by a 
vertical magnetic dipole on the ground plane is 
very poor. However, when the dipole is raised the 
situation is greatly improved. The height-gain 
function of the modes for this case is 


em? + REC n)e - ikC 2 
1+Ri(Cn) 


(kOn2)? 
5 + 


FA(2)= 





=1+72k2z/A,— 





(36) 


Since A, is a small quantity, the second term of this 
series may be quite large even for small heights. 

The height-gain function F%(z) describes the 
behavior of both the H, and EZ, components of the 
vertical magnetic dipole. The appropriate height- 





gain function for the horizontal magnetic field 
component H, is defined by the function 

C ikC ,2_. PPh g — tkC,2 — 

G?(z)= Gn( 2 > m(€ vge ) (37) 


ia 1+? 
in complete analogy to (21) for the electric dipole 
case. 

The excitation of the waveguide modes for a 
horizontal electric dipole was treated previously 
[Wait, 1960a]. For an electric dipole of moment 
Ids situated at z=2) and oriented in the direction 
¢=0, it was shown that 


E.,=cos ¢E,S (38) 


where E, and S are given by (3b) and (5), respec- 


tively. Also 
noH.=—sin ¢E)T* (39) 
where, explicitly, 
Tt=—T ett > bm Sm (Simp) f2.(20)fnx(2). (40) 


For most — the fields derived from S (i.e., 
the vertically polarized waves) are only significant. 
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The fields associated with 7* (i.e., the horizontally 
olarized waves) are hardly excited at all when 
baat which is the usual case. 
The field of a horizontal magnetic dipole of 
moment Kds has the same structure as the fields of 


the horizontal electric dipole. Explicitly 

H,=cos ¢H,S* (41) 
and 

E,=sin ¢H)T/n». (42) 


The solution for the horizontal magnetic dipole has 
been discussed recently by Galejs [1961]. He 
obtains his results by an application of the reciprocity 
theorem as outlined by Wait [1960a]. 

The remarks in the preceding section concerning 
the height-gain functions for the vertical dipoles 
can be carried over to the horizontal dipoles without 
modification. 

The generalization of the foregoing results to 
include the effects of earth curvature is not difficult. 
Provided one stays away from the antipode of the 
source dipole it is only necessary to multiply the 


field expressions by 
( : )? 
sin 6 


where @ is the central angle subtended by the source 
and observer at the center of the earth. The justi- 
fication for this step was given previously. It might 
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be added, however, that in addition to the restriction 
in the neighborhood of the antipode it is also neces- 
sary that 


C2>h/a (43) 


and 


(= FO 1, (44) 


Tt is not difficult to verify that these conditions 
are well justified at ELF. 
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Calculations of whistler ray paths in the outer ionosphere are shown for a variety of 
electron density profile models including exponential, constant, and columnar profiles. 
The Haselgrove formulation of the ray equations was used with the magneto-ionic repre- 
sentation of the wave refractive index to develop a set of differential equations for ray tracing 


suitable for inhomogeneous, anisotropic medium. 


The variation of paths with frequency, 


latitude, initial wave-normal angle, and other variables are examined for the purpose of 
providing a preliminary basis for comparison of the theoretical with some of the experi- 


mental results. 
1. Introduction 


Little work has previously been done to describe 
the mechanics and physics of electromagnetic wave 
propagation in an inhomogeneous anisotropic, me- 
dium. Of particular interest are certain phenomena 
taking place at frequencies below the plasma and 
gyrofrequencies of the ionosphere. The subject 
has become important during the last few years 
because of the increasing interest in VLF wave 
phenomena such as “whistlers.”’ 

Whistlers propagate in what is known as the 
“extraordinary” mode in magneto-ionic theory 
[Ratcliffe, 1959]. This mode of propagation is 
also described by the terms ‘‘whistler-mode”’ or 
“magneto-ionic duct.” Details of the phenomena 
and additional references are given elsewhere [Helli- 
well and Morgan, 1953]. he basic theory of 
whistler propagation was first given by Storey 
[1953]; however, his description is valid only for a 
restricted set of conditions, many of which do 
not apply over the complete whistler path. 

The work reported here is an attempt to demon- 
strate some of the complex phenomena of whistler 
propagation. The set of computations described 
here is in no sense complete or comprehensive; it 
is meant only to demonstrate a few interesting 
details of this mode of propagation and to stimulate 
further thought and work. Very little attempt 
has been made to interpret the form and shape of 
the results shown; such interpretations must arise 
from an examination of the form of the general 
solution of the ray path equations. 


2. Discussion 


2.1. General Procedure 


Equations which are in a form useful for compu- 
tation of ray paths in the ionosphere can be obtained 
from substitution of the familiar Appleton-Hartree 
expression for the complex refractive index into 


1 Based on a report prepared by Stanford Research Institute for Stanford Uni- 
versity under Prime Contract AF 18(603)-126 (Dec. 1959). 





ray equations derived by Haselgrove [1954]. The 
resulting equations in spherical coordinates r, 6, ¢ are 
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dt rsiné\y oo °* dt ee er ’ 

where 

u=the real part of the complex phase 
refractive index, 

Pr, pe, pg=the physical components of a vector 
of length y, and directed normal to 
the phase fronts (also called the wave 
normal or refractive index vector), 

t=time of phase travel along the ray 
(that is, (fAt)/e=number of wave- 
lengths in the medium along the 
ray path), 

f=wave frequency, 

c=speed of light. 


The quantity » and its derivatives are calculated 
from the Appleton-Hartree formula as follows: 


_2X(1—X—j2) 
D 


D=2(1—X—jZ)(1—jZ) —¥? sin? y+S 


M*=(u—jK)*=1 








S=+/Y‘ sin‘ y+4Y2(1—X—jZ)? cos? y 
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where 
M=complex phase refractive index=yp—jk, 
X=normalized electron density = (Ne?) /(w?me), 
Y=normalized magnitude of the earth’s magnetic 

~ > 

Y= (ueH,)/(wm), 


Z=normalized collision frequency = (v%)/w, 


field vector, 


y=angle between ‘ and Pim cos! [(p,Y,+ peYes 
+po¥¢)/(uY)I, 

N=electron density, 

e=charge on an electron, 

m=mass of an electron, 

w=2nf=2r- wave frequency, 

€)=dielectric constant of free space, 

= permeability of free space, 

y=collision frequency in collisions per second, 
K=imaginary part of the complex phase refractive 


index. 


The two values of M corresponding to the plus and 
minus sign on S represent the two modes of iono- 
spheric propagation commonly called the “ordinary” 
and “extraordinary” modes, 

The derivatives of u are 


Ou _ Ou OY Ou (piY cosy—uY; 
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where 1=the coordinates r, 6, ¢. 


[Nore: When ¥-0, du/d0¥-0, OY/Op—-- but dp/dp;—0.] 
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OM _ = ee eo 
of =ReS5= —Im (7p{ X+"—D 


[2—x—-2jz Ok) cos? Th) 


where Im=imaginary part of. 





ae 
[Nore: Oy/07 is calculated holding p constant and so measures 
the change in the direction of the earth’s magnetic field in 
space.] 


The quantities 0X/0i and 0Z/di must be derived 
from the space variation of Y and Z which are 
assumed. An arbitrary inhomogeneous ionosphere 
can be represented in this manner. If an earth- 
centered dipole field approximation is used for the 
magnetic field, then 


ee ene 
Y=f!, (2) V¥1+3 cos? 6 


where Y,=Y at the equator at the surface of the 
earth ~(8.7-10°)/f. The vector components of Y 
in spherical coordinates are 


¢ r ta 
ia ar__3¥ 
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= —3Y cos 6 sin 6 
00 =. 1+ 3 cos? 6 
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oo iba. 
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where 
é=colatitude and 


¢=longitude. 


The derivatives of y for a dipole magnetic field are 


oy _ 
5 = 

oy _Oy dy 
6 dy 08 


He 
where y=angle between Y and the radial (see fig. 1). 


|NorE: y=tan-! (—% tan 6) for the dipole field.] 
oy | 
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[Note: When y¥-0, both numerator and denominator of 
Oy/Oy go to zero and 0y/07—1.] 
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FigurE 1. Coordinate system for two 
dimensions, r-0. 
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It is usually desired to compute the following 
additional quantities along with the ray path: 
(1) Path length, S 


dSol ey 1 1 
—_=— On/O)? = oer = . 
dt we” + (On/Ov) ray refractive index yu cosa 





(2) Time of travel, 7, also called group delay 


a7’ _1 (1 fm 
dt ¢ wot 
Qu > OM Se ie 


9 
[2-2 jZ—jZX +5 Y* cos? y 


(—X—j2)(1+X) | }) 


(3) Absorption, A, in nepers 


dA Kk. 


ra w=2nf 


2.2. Notes on the Method of Digital Computation Used 
a. Initial Conditions Needed 


Initial values for 7, 0, ¢, p;, ps, ps are needed. To 
get p,, pe, and py, the wave normal direction can be 
inserted as an initial condition. 4 is then calculated 
and the p’s obtained. For example, let 2;=(p;/p)= 
direction cosines of wave normal with respect to the 
local coordinate axes. Then, since p=p, py=pa;. 


b. Stability Condition 


The six basic ray differential equations contain 
one more degree of freedom than there are con- 
straints explicitly expressed. That is, it only takes 
5 quantities to define a vector field (three space 
coordinates and two direction cosines). The sixth 





, 6a ee ial >| 
equation is implicit in the definition of |jp|\=y. The 
round-off errors involved in calculating » and in 


I>] 
integrating dp,/dt will result in pl <u after a few 
steps of integration unless an additional constraint 
is imposed. Thus, to maintain the equality, we 
correct the p,’s at each point in the manner 


__b 
Pt correctea— Pi 
p 


This correction is made as early as possible in the 
derivative calculation, i.e., immediately after caleu- 


J ive ; >. 
lation of M. Since the magnitude of p is not used 
in the calculation of M (only the wave normal direc- 


. . . . a 
tion is used), no errors arising from |p| #y are allowed 
to accumulate by this procedure. 


c. Interval Size 


The Af to use depends upon the accuracy required, 
the frequency, the size of refractive index gradients 
encountered along the ray path, and the stability 
and accuracy of the numerical method of solution of 
the differential equations used. For the whistler 
ray paths, shown later in this report, the interval 
size was changed when the local truncation error 
became either too large or too small. The local 
truncation error was measured by noting the differ- 
ence between the predicted values for the next point 
and the values obtained by application of the first 
corrector (an Adam’s predictor-corrector method of 
solution of the D.E. equations was used). When 
this difference became less than 5-107", the interval 
size was doubled using an extrapolation formula in a 
special subroutine. When the difference became 
greater than 10~‘, the interval size was halved. The 
corrector was applied only once every 5 or 10 points 
and only a predictor was used for the in-between 
points. This was done because application of a 
corrector and subsequent doubling or halving costs 
several “points” worth of time. Thus it is undesir- 
able to do this very often. The interval size used 
for these whistler calculations was approximately 
one-half an earth radius. Calculation was started 
at the top of the F, layer. 
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d. Other Comments 


The bulk of computation time is spent in calculat- 
ing the derivatives of M, yu, etc. Only a small frac- 
tion of the total time is actually spent in the numeri- 
cal method of solution calculations. Therefore, a 
method of solution should be chosen that minimizes 
the number of times the derivatives have to be 
calculated, and that makes maximum use of the 
derivatives already calculated. Thus the Runge- 
Kutta methods are long and cumbersome for these 
equations. A predictor-corrector method with a 
minimum number of corrector cycles applied seems 
to be about the best. Perhaps even a high-order 
predictor (5th or 6th order) might be best, since this 
would increase the interval size one could use without 
increasing the calculation time proportionately. 

The main weakness of the program used for these 
calculations is the instability of the extrapolation 
formula used for doubling. This formula has to 
jump ahead two intervals. It would be better to 
carry more points along with the calculations and 
use a higher-order extrapolation formula that only 
jumps ahead one space—or better still, carry suffi- 
cient points to use no extrapolation formula at all. 

If these equations are to be used for frequencies 
high enough that the earth’s magnetic field can be 
ignored (frequencies much higher than the local 
gyrofrequency) then a separate routine should be 
written for this case, since the equations simplify to 
a much easier and faster computation. If the same 
routine is used, many of the expressions for the 
derivatives become indeterminate when Y=0. Most 
of these derivatives actually become zero, but special 
tests would have to be installed in the program to 
handle this case. 


2.3. Computations 
a. General Procedure 


The differential equations were programed for 
numerical solution on an Electrodata 205 computer. 
Only the two-dimensional case of propagation in the 
magnetic meridian was utilized in these computa- 
tions. Also, no losses were included (v=0). The 
calculation of a whistler path was begun with the 
assumed initial conditions at an altitude of 300 km. 
This altitude corresponds roughly with the maximum 
ionization level of the F, layer. 

An electron density profile was assumed for alti- 
tudes greater than this. 

The computations are organized into the following 
categories according to the electron density profile 
used and the parametric variation examined: 


(1) Exponential model—variation of frequency, 

(2) Exponential model—variation of initial latitude, 

(3) Exponential model—variation of initial wave 
normal direction, 

(4) Constant density model, and 

(5) Miscellaneous. 


The exponential model used in all of the first three 
categories is: 








N=180,000 exp [—4.183119 (r—1.0471)] 
where 


N=electron density in electrons per cubic centimeter. 
r=radial distance from center of earth in earth 
radii. 


(r=1.0471=300 km above the surface of the earth 
at the beginning of the path computation). 

These particular values of electron density and scale 
height were obtained from Maeda and Kimura [1956]. 


b. Exponential Model—Variation With Frequency 


Figure 2(a) gives the results of ray path computa- 
tions at the frequencies 5, 10, 20, 25, 30, 50, 100, 
200, 400, 1000 ke/s. Figures 2 (a) and (b) show the 
ray path for each of these frequencies. The initial 
wave normal direction for each of these paths is 
vertical (wave normal angle=0°). In order to 
exhibit the variation of certain variables in these 
computations, figures 2(c) to 2(i) have been plotted. 

Some of the interesting features of these computa- 
tions may be listed as follows: 


c. Exponential Model—Variation With Latitude 


Figure 3 gives the results of ray path computations 
initiating at north latitudes of 10, 20, 30, 35, 40, 45, 
50, and 60°. Figures 3 (a) to (i) show the ray paths 
for each of these latitudes as well as the variation of 
wave normal direction over the path. Additional 
data from these computations have been plotted in 
figures 3 (j) to (0). 


d. Exponential Model—Initial Wave Normal Variation 


Figure 4 gives the results of ray path computations 
with initial wave normal angles 10, 5, —5, and —10° 
measured positively clockwise from vertical (along 
the radial). Figures 4 (a) to (d) show the ray paths 
and wave normal direction over these paths. The 
wave frequency is 10 ke/s. Additional data from 
these computations have been plotted in figures 4 (e) 
to (g). 

The behavior of the wave normal over the path 
changes only negligibly with these changes in initial 
wave normal direction. 


e. Constant Density Model 


Figure 5 gives the results of ray path computations 
with constant electron density. The characteristics 
of a constant density model are as follows. Since 
there is no electron density gradient, the total refrac- 
tive index gradient is formed by the magnetic field. 
Also, since the refractive index gradient determines 
the change in wave normal which in turn determines 
the bending of the ray path, the bending of the ray 
path is determined only by the guiding effect of the 
magnetic field, not by the electron density. Thus 
the ray path is the same regardless of the magnitude 
of the electron density. 
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f. Miscellaneous Computations 


Comparisons with Maeda’s computations—Maeda 
(1951) used a very approximate method of numerical 
integration as well as an approximation to the refrac- 
tive index gradient which is only very roughly 
correct. In order to compare his computations with 
the more exact formulation of Haselgrove, one of 
Maeda’s published ray paths was compared with a 
similar computation using the present program in 
figure 6. It can be seen that the difference is 
appreciable. The lack of comparison is not due to the 
quasi-longitudinal approximation of the refractive 
index function which was used by Maeda since this 
approximation was used also in Haselgrove’s formula- 
tion for this particular computation. The difference 
between the Q. L. approximation and the full expres- 
sion is too small to be shown. The differences 
between the two computed paths may be due to the 
relatively crude method of integration used by Maeda 
to numerically solve the differential equations. 

Field-alined column model—In a recent theoretical 
paper it has been shown that if the ionization is 
field-alined it can act like a waveguide in trapping 
whistler energy [Smith, Helliwell, and Yabroff, 1960]. 
This would explain some of the very distinct and 
clear whistlers which have been recorded. In order 
to answer certain questions regarding the effects of 
field-alinement, the following models have been 





N=180 o00f 1+C exp | oo ~ |} 


exp [—4.183119(r—1.0471)] 


where 


b=r/sin2@ (the equation of a field line is b= 
constant), 

b)=2.094 (this is a field line corresponding to 300 
km altitude at a latitude of 45° and is the 
center of the column), 

C=the modulation factor, i.e., the relative 
increase of maximum ionization at the 
center of the column over that of the sur- 
rounding background level of ionization, 
and 

D=the standard deviation of the column. 
(The column is of gaussian electron density 
distribution with the standard deviation 
measured in terms of the local distance 
between the field lines.) This makes more 
intuitive physical sense than making the 
thickness of the column constant at all 
points on the field line, since the individual 
electrons will tend to diffuse easily along 
the lines of force but not across them. 

For these computations D=0.029 which corres- 
ponds to a standard deviation of about 35 km at an 
altitude of 300 km, and about 215 km at the top of 
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Figure 2(c). Mazimum altitude versus frequency. 


The maximum height of the path remains essentially constant up to the frequency of minimum time ‘delay [the nose 
frequency, which is seen to be at about 20 ke/s in fig. 2(d)] at which point the maximum path height begins decreasing. 
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The time delay decreases in the expected ‘“‘ Eckersley” fashion until the nose rerion is reached. 
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FIGURE 2(d). 
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Total path delay versus frequency. 
This region corresponds ;to 


the frequency reaching approximately 45 percent of the gyrofrequency at the top of the path. These calculations confirm this 


behavior. 


The time delay at higher frequencies increases until the effect of the shortening of the path length at high frequencies 


overbalances the decrease in group velocity to cause the time delay to begin to decrease again above 200 ke/s. 


the ‘‘Eckersley”’ region. 
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FIGURE 2(e). 
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Final latitude versus frequency. 


The final latitude remains fairly constant at frequencies up to the nose frequency. This suggests some latitude focusing in 


Above the nose frequency, the final latitude starts shifting toward the pole until the effect of the 


shortening path length at higher frequencies overcomes this tendency and results in a shift back toward the equator. 
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Ficure 2(f). Wave normal—field line angle along the path. 


At frequencies below the nose frequency, the wave normal swings steadily 
below the field Jine as the path progresses. As the nose frequency is approached, 
the wave normal tends to oscillate about a position approximately ten degrees 
above the field line. At higher frequencies, the wave normal swings outward 
after an initial dip toward the field line. 
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FiGuRE 2(g). Wave normal—field line angle along the path. 
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4 Figure 2(h). 
along the path. 
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In the first two computations shown in figures 7 
and 8 the initial wave normal direction is vertical. 
When the modulation is 0.2, good guiding occurs. 
When the modulation is reduced to 0.1, the ray is not 

uided because the initial wave normal direction is 
ar from the field-line direction and the ray gets 
pulled out of the column before the wave normal gets 
properly lined up with the field line. This is in 
agreement with theoretical predictions [Smith, 
Helliwell, and Yabroff, 1960]. In the third case 
shown in figure 9, the wave normal was set along the 





field line. In that case, the guiding is even better 
than for the higher modulation of case 1. The final 
case shown in figure 10 shows the modulation reduced 
to zero with the initial wave normal direction on the 
field line. The guiding of this ray is still good, but 
less so than for the previous case of 0.1 modulation. 

These calculations support the theoretical pre- 
dictions [Smith, Helliwell, and Yabroff, 1960] that 
wave components whose wave normals lie within a 
certain angle with respect to the direction of the 
magnetic field will be trapped. 
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The maximum height of the path is seen to increase very rapidly with increas- 
ing initial latitude until an altitude is reached at which the plasma frequency 
becomes equal to the wave frequency. This represents the usual reflection level. 
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Figure 3(k). Total path delay versus initial latitude. 

The total time delay is seen to be fairly constant over a wide range of initial 
latitude. This may be a bit surprising since it was previously observed that the 
path length changes very rapidly with initial latitude. Clearly the average 
group velocity must be much higher for the longer paths which are at a higher 
latitude. 
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Figure 3(l). Final latitude versus initial latitude. 


This plot of final latitude against initial value shows some focusing of rays in 
the region of 45° and 33° latitude. 
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At initial north latitudes greater than 35°, the wave normal swings from the 
outside to the inside of the field line along the path. For initial north latitudes 
less than 35°, the wave normal swings farther away from the field line. 
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The total time delay is also seen to remain roughly constant with changes in 
initial wave normal angle. This is again consistent with the idea of average 
group velocity being roughly proportional to path length. 
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Final latitude versus initial wave normal angle. 


The final latitude is seen to remain constant with initial wave normal variation 
even though maximum path height changes. 


Figure 5. 


Constant density model. 
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On the Analysis of LF Ionospheric Radio 
Propagation Phenomena* 
J. Ralph Johler 
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(Received April 27, 1961) 


Recent theoretical work which employs the classical magneto-ionic theory for a special 
model of the ionosphere applicable to transmission via the ionosphere at or close to grazing 
incidence is employed to analyze LF propagation data. The results of the analysis illustrate 
a practical model of the ionosphere by a detailed study of transmission via the first time- 
mode in particular. 


1. Introduction 


Recent contributions and extensions to the classical geometric-optical theory of propaga- 
tion [Wait and Conda, 1958; Wait, 1960a; Johler and Walters, 1960, 1961] of LF waves about 
the earth provide a valuable analysis tool for the interpretation of various experimental data. 
This paper summarizes the mathematical formulas required to predict the field with particular 
emphasis on a type analysis most suitable for large scale electronic computers. Attention is 
focused on certain available experimental data and the results of the analysis are detailed with 
emphasis on technique. However, the physical phenomena which can be deduced as a result 
of the analysis of these data are given considerable attention to emphasize the value of the 
technique. 


2. Mathematical Theory 


A geometric-optical theory [Bremmer, 1949] can be employed with the aid of a sharply 
bounded model ionosphere [Johler and Walters, 1960, 1961] emplaced at various altitudes h 
above a spherical earth to describe the ew field, H(w,d), radiated from a dipole source current 
moment, J,/. It is common engineering practice at LF to relate the amplitude of the field, 
|\E(w,d)| to the well known concept of radiated power [Ballentine, 1924], neglecting earth 
losses, 

P,=1.6(10-") w?(Ipl)?/Zo, (1) 
where Z)~1207, a constant. Such a field, E(w,d), implies a similar Hertzian dipole receiver 
(vertical polarization). The total field, Z(w,d), is the sum of 7+1 reflected ionospheric waves, 
j=0, 1,2,3..., 


E(w,d)=S2E;(w,d), (2) 
j=0 
in which the zero order (j=0) reflection, E)(w,d) is the groundwave. Each skywave reflection, 
j=1, 2,3..., using the positive time function, exp (iw), can be written, 
E;(w,d) =iwdD;'C exp (it) GiGi aF jC), (3) 
where [Johler, Kellar, and Walters, 1956], 
E,(w,d) =[E(w,d)] exp [—i(wb+¢,)], (4) 
and where, 
C= Iplb?/4axd?= (10-")/d,(Igk=1). (5) 
The local time '! for the groundwave, ty is, 
ts=t—b, (6) 
b=nd/e, (7) 


*This is a revised version of a paper presented at the Avionics Panel Meeting of the Advisory Goup for Aeronautical Research and Develop- 
ment (AGARD) NATO, Istanbul, Turkey, October 1960. 
1 The notion of local time was apparently introduced by Lorentz [1906, p. 57]. 
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FicurE 1. Relative sky wave delay for various altitude, h, of the lower boundary of the model ionosphere. 


where c is the speed of light, e~3(10°) m/sec, and m is the index of refraction of air, m~1. 
Similarly, the local skywave time, or skywave delay, t}(j=1, 2,3 . . .) is, 


tj=t—b,, (8) 
b,=9,D),/e, (j=1, 2,3... .). (9) 


The quantity, 6;—), is frequently called the relative skywave delay, figure 1 (relative to the 
groundwave). The physical length of the ray, D;, can be evaluated geometrically for a reflec- 
tion at an altitude, 4, above the surface of the earth of radius, a, 


D,=2j [(a+h) cos ¢;,;—a cos 7,], (10) 


where, figures 2, 3, ¢; is the angle of incidence of the ‘ray’ on the ionosphere and 7; is the 
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corresponding angle of incidence on the earth and the subscript ‘‘j”’ reminds the reader that the 
equation refers to the particular ionospheric reflection under consideration. The angles 
¢; and 7; are evaluated quite simply from the geometry, 


sin ¢; ;=Aj! a sin a (11) 

COs $;,;=Aj' [a(1—cos 6/27) +A], ~ (12) 

sin r;=Aj;! (a+h) sin 6/2), (13) 

cos 7;=Aj![a(cos 6/27—1)+h cos 6/27), (14) 
A;=[2a(a+h)(1—cos 6/27) +h?J}, (15) 


where @ is the angle at the center of the earth subtended by the distance, d, along the surface 
of a spherical earth, or, simply, 


d=ae. 


The factors, Gi and Gj refer to the transmitter and receiver complex antenna patterns re- 
spectively, and again the subscript, ‘j’’, is employed to designate the particular time-mode 
under consideration. The vertically polarized Hertzian dipoles considered in this paper reduce 
to, Gj=G5=sin 7;. Of course, other antenna complex pattern factors can be introduced. 

Since plane reflection coefficients are employed in this analysis, the convergence-divergence 
coefficient, a;, is necessary to convert the plane reflection coefficients to spherical reflection 
coefficients or from the viewpoint of “ray theory’ it is necessary to take account of ray focusing 
by the ionosphere and a corresponding defocusing of rays by the earth. The ‘‘classical’’ con- 
vergence-divergence is modified by a “convergence correction” [Wait, 1959], A,, in this paper, 
since the values of the first time-mode studied are at distances close to the “caustic” of rays or 
geometric-optical horizon, and hence the complete expression, a, 


aj=(1+h/a)[(27 sin 6/27)/sin 6]! {[a(1—cos 6/27) +h]/[(a+h) cos 6/2)—a]}#A,, (16) 
in which A; can be evaluated from the cylindrical Hankel function of order, n=}, of the second 


kind, 


1 
A,~ [5 kya cos® 74/3 sin? + J HY [ia cos? 7,/3 sin? r| 


Xexp { —1 [ 5x/12—ka cos? 7 ,/3 sin? ni} (17) 


The factor,|A;, approaches unity (A;~1) at shorter distances, (77/2). Near the geometric- 
optical-horizon, 7;~7/2, and beyond, the correction becomes quite appreciable. The Hankel 
functions, H\?(z) can be readily evaluated on an electronic computer by a consideration of the 
integral form, 


HP (2)=* a exp (—7z sin p+inp\de—2; | exp (—z sinh p) {exp np+exp (—np+inz) }dp, 
(18) 


where the order, n=4, and the argument, z=h:a cos* 7,/3 sin? r,;, in which the wave number, 
ki=wm/c~w/c, for the case under consideration are real numbers. Nevertheless, it is inter- 
esting to note that the method is also applicable to complex argument, z, and complex order, 
n, provided the multiple branches of the Hankel function are carefully considered. The follow- 
ing substitution is made, 


B=exp (—p), 


and a resultant expression for the Hankel function. 
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“ij Pa ~ 
Hp (= { exp (—iz sin p+inp) dp—= | Ga exp (ins) | exp |-§ (5-6) | dg, 
(19) 
or, 


°F *1 
H®(z)= | fle,e.nde— | F(z,8,n)dB, (20) 


is found to comprise integrals with finite limits. These integrals can then be evaluated in terms 
of Gaussian quadrature,” 
oe : ; a : 
H® (2)=>3 W2, [Ref(2,om,n)+% Im f(2,pmsn)]— 5 WS [Re f(2,8m,n) +7 1m f(2,Bm,n)|+eQ2) (21) 
m=1 m=1 


m=1, 2,3... M, where, e(M) can be made arbitrarily small by increasing M, 


“ss 


where, 


W? =cH,,/2, 

Wé=H,,/2, 

Pm=TLy,|/2+27/2, 

Bn=Zm/2+1/2, 
and the constants, z», of the theory of Gaussian quadrature can be evaluated as the roots of 
the Legendre polynomials defined by, 


qd’ 
ao i )¥—2™M ! Py(z)=0, (22) 
or, 
P,(z)=1, 
P; (x) =, 
4 ] 
P, (x) =50°—5? 


P;(z)=52°—5 2 
i a 


35 15 3 
P,(2) a hada dD 


_ (23) 
where polynomials of higher degree are determined by the use of the recursion formula, 
(M+1)Py4i(z) +MPy-1(x)— (2M+1)2Pu(x)=0. (24) 


The weight coefficients, H,,, are evaluated from the roots, 2m, 
H,,=2/(1—a}) [P(tm)}?. (25) 


The factor, F,, accounts for the presence of the earth at the transmitter (source) and re- 
ceiver. If it can be assumed that the “ray” reflected by the earth is not too close or beyond 
the geometric-optical horizon, the Fresnel approximation of the ground reflection coefficient 


suffices to determine F,, or, 


Fy~l1+Ri(r +R 2 (75), (26) 





‘This type analysis has been used extensively by the author in previous papers to evaluate Fourier integrals [Johler and Walters, 1959]. 
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where the superscripts, ¢ and r, refer to the transmitter and receiver respectively and the 
subscript, e, refers to vertical polarization, and, 





Rr s)= {3 cos 75/k}—[k3/ki—sin? 7,]$} /{k3 cos 7 5/k?+[k3/k}—sin? 7,}}}, (27) 
or for horizontal polarization, 
R,(7s)= {cos 7;—[k3/kj—sin? 7,}*} / {cos 7,+[/3/kj—sin? 7,}3}, - (28) 
where, 
oa ¢c 2) 4 


and ¢, is the dielectric constant of the earth (e,.~15) and o is the conductivity of the earth 
(e~0.005 for land). Since the examples considered in this paper involve distances which, 
for the first time-mode, were close to the geometric-optical horizon, it was found necessary to 
replace the Fresnel reflection coefficients with a more rigorous treatment. 

Close to the geometric-optical horizon a calculation of F;=FjF; can be accomplished by a 
numerical evaluation of the contour integral, [Wait and Conda, 1958], 


Fi? ~ (x)—3 exp [—iha6’] exp [—7(kia/2)*0’p]/W;(p)—qWi(p) dp, (30) 


exp [—i2x/3] 


6’ = (d—dy)/a, where dy is the distance * from transmitter to the geometric-optical horizon 
(r;~7/2) and, 
» fka\3 ky ki ; 
-~i(%) ks IT (31) 


in which o=g; is the ground conductivity at either transmitter Fj or receiver F; and e, is the 
corresponding dielectric constants (relative to a vacuum), and, 


“ 


W,(p)=exp [—2ni/3] yx/3 (—p)*H? [2/3(—p)}], (32) 
? ° 79 1 -3 77 (2) 2 a (9) 2 3 
W;(e)=exp [—2ni/3] yx/3 a2 (—p) 3H; 3(—e)! +pH; 3(—p)} , 
W's (0)=exp [—i42/3] yx/3 pH Ea (33) 
MK? Z)=% cue (Z) =exp [—i2/3] HP (Z) sz HY? (2). (34) 


The Hankel functions Hj” (Z) can be calculated by previously described quadrature methods 
(20-27). The contour integral can be written [Wait, 1960b], 


a ee Te —— © exp [—*(hia/2)*6"a] 5 exp [—2(kia/2)30’a’] , , 
Ft'=(x)—* exp [—ikao’) { W (a) —qWi(a) da— | * W? a’) —qWila’) da’ >» 








(35) 
a’ =a exp (—i2z/3). These integrals can then be evaluated with a Gaussian quadrature. 
Close to but beyond the geometric-optical horizon (r;~7/2) [Wait and Conda, 1958], a 
calculation can be performed as a residue series summation, 


=, exp [—i(ha)'0’77] , 
Frt~—Qivn exp [—ik,(d—dy)] p> Er ACES (36) 





3 a ~ 6.367 (108) m. In diffraction problems it is common practice to account for the vertical lapse of the permittivity of air at the earth’s sur- 
face by use of an ‘‘effective earth radius’’, a, in place of a such that £~0.75-0.85, [Bremmer, 1949, p. 145]. 
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where the complex conjugate 7, of 7* comprises the special roots of the differential equation 
of Riccati tabulated by Johler, Walters, and Lilley [1959], employing the time function, 


exp (— wit) ’ 

















dé, 2 pe ‘ 
Tp 2rd 1=0, (37) 
where 
Pye ikgad/k? ; Zo 
“ae 2 Bi ser 
(oa)'] 1 ]*— Gaays] =e]? sl 
? o— We i 
i 
and, k,=* m~o k, implies the complex conjugate of k2, (3), and the limiting roots are found ; 


from the equations, 


HY E (—2r,)# 0, ime, 


PERL: AIR Asie 


HY [5 (—27.)! Jeo, 5.=0, 
— ee 


Finally, the effective reflection coefficient of geometric-optics, C;, can be evaluated for 
vertical polarization, 





1 di 1+ Az ‘ 
Cm FeLi aa. (39) 
where, 
A,=—R, y Mae 
A,=R,T,.+ Rul nm 
A,=R.R,[— T.2T amt Vitel: (40) 


Thus, the essential nature of the propagation about the earth via the ionosphere can be described 
in terms of four reflection coefficients, 7... Tem, Tne, and Tnm. The reflection coefficient, T,,, 
refers to the vertical electric polarization of the incident plane wave and a similar vertical 
electric polarization of the reflected wave. The coefficient, T.,,, describes the generation of the 
abnormal component by the incident vertical polarization (reflected horizontal polarization 
for vertical excitation). Similarly, 7,,,, refers to the incident horizontal electric polarization 
and the corresponding reflected horizontal electric polarization. Also, the abnormal componeut 
generated by the horizontal electric polarization (reflected vertical polarization for horizontal 
excitation) is described by the coefficient, T,,... Thus, figure 3, 


T..=E,,/Eys, Tan =Ey,|Exs, Tn =Ex, Ey Tne= Ey [Ex (41) 


where the subscripts 7 and 7 refer to incident and reflected wave respectively at the ionosphere 
boundary. The calculation of these reflection coefficients has been previously described in 
detail and the results of various computations have been tabulated [Johler and Walters, 1960]. 

The reflection coefficients are completely specified by the angle of incidence, ¢;, figures 2, 3, 
electron density N(El/cm*), figure 4, collision frequency v(c/s), figure 5, magnetic field intensity 
H,, (gauss), magnetic dip J and magnetic azimuth ¢,. The effective reflection coefficients, 
C;, can then be written in terms of the ionosphere reflection coefficients, 7’, and the ground 
reflection coefficients, R, for each local ionosphere or ground reflection, 


C=C, (71)= 1 1m 

C= C,(72) we Tee, Teele, +Ra, Tom, T ney: 

C= C3(r3) =R.,Rn, Y ay Tn t+R.,Re, 1 Tom, Test Re, Re, Tec, on Deeg + Rim, Ring Tnmg 1 ne T ney 
(42) 
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> 7 4 ‘yy ry ry 
C= C474) =R eMeglie, Da, Vy P40, +f Re R.Rng y a Toms y P.0, +R. Rn Re, 1 T ney Dems Dee, 
> lg } 4 al 7 
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> > AS 
+Rn fi malt €3 y Pm Tneg Tum, Tem, + Ry PRngRng Tne Twins T nme Doms (43) 
where the subscripts 1, 2, 3, 4... refer to the first, second, third, fourth, etc., reflection 
regions of the ionosphere with corresponding separate values for electron density, N, collision 
frequency, v, and the intensity and direction of the earth’s magnetic field vector, 7,,, or, for 
the ground, refer to the separate values of conductivity, o, and dielectric constant, é:. 
For a homogeneous ionosphere and ground with a single set of values, N, v, Hp, ¢, and 


Y 


€2, the effective reflection coefficients, C;, reduce as follows: 
O- FT, 
Cy T3284 Ra TT ne 
C3=2R Rn T eT mT met RET 2+ B37 mm Tem T mes 
C= BRT Tet BRT eT el eet PRRET oT xml eda t Ret ad nl not RET%., 


(44) 
The wave number of the ionosphere propagation medium, 
Ww 
k= Noy 
® é 
ig * Ney (45) 


implies two upgoing waves excited in the model ionosphere with index of refraction, -yo, for 
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of the calculations was the altitude, h, of the assumed sharply bounded model ionosphere. 
This implies a variation of the electron density according to the N(h) profile as the variable 
his changed. Also, small adjustments in the angle of incidence, ¢;, are implied. 

yeophysical data on the ionosphere employed in the calculation are shown, tables 1, 2. 
The magnetic data, shown in table 2, were scaled from Magnetic Charts [1954, 1957]. Appar- 
ently a nominal intensity value H,,=0.5 gauss (table 2) is an adequate assumption for many 


practical cases. 





TABLE 1 
h v Na Na No 
km c/s, X10*6 El/em3 El/em’ El/em3 
» 3 175 ‘ eee 10 
60 | 50 ee poe 1000 
65 24 | 10 1205 2400 
67.5 | 16 | 56 | See eee 
7 | 10.9 150 | 2675 | 5200 
7 | 92 | . == sciecs : 
72.5 7 200 |. é : 
7 | 4.5 160 4505 8850 
77.8 2.6 140 |__- e Saee 
80 1.6 500 | 7950 15400 
gl 1.2 3500 | 10550 | _- 
82.5 0. 88 11000 | 16100 |..---- * 
85 . 465 15000 | 20500 26000 
90 .155 10000 | eee a are - 
TABLE 2.— Magnetic field, Hm, data 
Adak-Nome, d=965.5 statute miles (1,550 km) 
Reflection , I da 
number | 
Gauss deg deg 
1 Ist 0. 5187 | 68. 68 12. 27 
2 Ist . 4996 | 65. 75 12. 52 
2 2d . 5329 | 71. 52 11.72 
2 Ist . 4952 64. 82 12. 63 
3 2d . 5187 | 68. 68 12. 27 
3 3d - 5388 | 72. 46 11. 66 
Adak-Kodiak, d=1040.1 statute miles (1,670 km) 
] Ist 0. 5035 | 67.18 | 51. 08 
2 Ist . 4950 | 65.15 52. 84 
2 2d . 5170 | 69. 04 | 49. 73 
3 Ist . 4923 64. 42 | 53. 42 
3 2d . 5035 | 67.18 | 51. 08 
2 3d 5262 | 69. 64 | 49. 35 





The modal (skywave time-mode) and the total (including the groundwave) calculated 
fields (amplitude and phase) for quiescent, disturbed, and blackout conditions of the ionosphere 
are presented, figures 6 to 9. The effect of the ionosphere focusing corrections, A;, is also 
illustrated, and, indeed, was found to be appreciable, both for the total field, E(w,d), and the 
first modal field, Z;(w,d),j7=1. The diffraction correction F; was also modified from the clas- 
sical Fresnel reflection coefficient as discussed with the aid of calculated data of Wait and Conda 
[1958] which was installed in the electronic memory of a small computer.® 

The ground wave, /y(w,d), which was calculated with the classical Bremmer-van der 
Pol theory [Johler, Kellar, and Walters, 1956] and corrected for land-sea boundaries with the 
aid of the Millington [1949] method is not shown alone graphically, since the groundwave is 
not a function of the emplacement altitude, h, but can be summarized as follows: 








TABLE 3 
Eo(w,d) 
ees ___| arg Eo(w,d) 
| 
IJ= | Iol=2050 
v/m | vim Radians 
Nome..- ; 7.56(10-% 1.55(10-5) | 2. 30 
Kodiak. _ .| 4.64(10-*) 9.51(10-) | 2. 43 





5 TBM-650. 
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where the subscripts 1, 2, 3, 4... . refer to the first, second, third, fourth, etc., reflection 


regions of the ionosphere with corresponding separate values for electron density, N, collision 
frequency, v, and the intensity and direction of the earth’s magnetic field vector, /,,, or, for 
the ground, refer to the separate values of conductivity, o, and dielectric constant, é:. 





For a homogeneous ionosphere and ground with a single set of values, N, v, H;,, ¢, and 
€, the effective reflection coefficients, C;, reduce as follows: 


C\ == re 

C,=72,R,+RnTomTme 

gay fey iy AR Ty Be Dh py Sy Oe 

C,=R,R2T in Te + 3R2R eT 2, Tem Tme-+ 2B Re, Tce T nm T em Tet R37 2m T om T met FET be, 


(44) 
The wave number of the ionosphere propagation medium, 
® 
ie z 
‘ c Noy 
ie) ~ 
=" ne, (45) 


implies two upgoing waves excited in the model ionosphere with index of refraction, mo, for 


595362—61——-7 513 








of the calculations was the altitude, A, of the assumed sharply bounded model ionosphere. 
This implies a variation of the electron density according to the N(h) profile as the variable 
his changed. Also, small adjustments in the angle of incidence, ¢;, are implied. 

Geophysical data on the ionosphere employed in the calculation are shown, tables 1, 2. 
The magnetic data, shown in table 2, were scaled from Magnetic Charts [1954, 1957]. Appar- 
ently a nominal intensity value H,,=0.5 gauss (table 2) is an adequate assumption for many 


practical cases. 


























TABLE 1 
h v Ny | Na No 
km | e/s,X10*6 El/em’ El/em’ El/em3 
55 | 175 PPA Sa) EPR at Shiro. 10 
60 | 50 1000 
65 24 | 10 | 125 | 2400 
67.5 | 16 5 leer aa oe : 
70 | 10.9 | 150 | 2675 | 5200 
71 9.2 =| gt Ee an Sine 
72.5 7 aie asta I 
75 | 4.5 160 | 4505 8850 
77.5 | 2.6 | Ea: Eee ; 
80 | 1.6 500 | 7950 15400 
gl 1.2 3500 | 10550 |___- 
2.5 | 0. 88 11000 | i eae 
85 . 465 15000 | 20500 26000 
90 .155 de! ee a 
TABLE 2.— Magnetic field, H,,, data 
Adak-Nome, d=965.5 statute miles (1,550 km) 
Reflection Ha I | da 
number 
Gauss deg deg 
1 Ist 0. 5187 68. 68 12. 27 
2 Ist 4996 | 65. 75 12. 52 
2 2d 5329 | 71. 52 11.72 
2 Ist 4952 64. 82 12. 63 
3 2d . 5187 68. 68 12. 27 
3 3d . 5388 72. 46 11. 66 
Adak-Kodiak, d=1040.1 statute miles (2,670 km) 
1. Ist 0. 5035 | 67.18 | 51. 08 
2 Ist . 4950 65.15 52. 84 
z 2d . 5170 69. 04 49. 73 
an Ist | . 4923 | 64. 42 | 53. 42 
7 2d | . 5035 67.18 | 51. 08 
| en 3d | . 5262 | 69. 64 | 49. 35 
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The modal (skywave time-mode) and the total (including the groundwave) calculated 
fields (amplitude and phase) for quiescent, disturbed, and blackout conditions of the ionosphere 
are presented, figures 6 to 9. The effect of the ionosphere focusing corrections, Aj, is also 
illustrated, and, indeed, was found to be appreciable, both for the total field, E(w,d), and the 
first modal field, E;(w,d),7=1. The diffraction correction F; was also modified from the clas- 
sical Fresnel reflection coefficient as discussed with the aid of calculated data of Wait and Conda 
[1958] which was installed in the electronic memory of a small computer.® 

The ground wave, /(w,d), which was calculated with the classical Bremmer-van der 
Pol theory [Johler, Kellar, and Walters, 1956] and corrected for land-sea boundaries with the 
aid of the Millington [1949] method is not shown alone graphically, since the groundwave is 
not a function of the emplacement altitude, h, but can be summarized as follows: 














TABLE 3 
Eo(«,d) 
bicie== ace | arg Eo(w,d) 
| 
| Il =1 Tol =2050 | 
ear | , 
v/m | vjm Radians 
Nome..-----.--| | 7.56(10-®) | 1.55(10-5) | 2. 30 
Kodiak. __.-__-|  4.64(10-%) 9.51(10-%) | 2. 43 
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All computations were made with the assumption of a dipole current moment, J)/=1. 
The total radiated power, Pr=1.3 kw, for the arctic data considered was found to correspond 
to a dipole current moment, Jp/=2050 amp-m; therefore, the calculated field amplitude, [Z(w,d)], 
or [E£;(w,d)], has been increased by this factor for comparison with the data. 
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Figure 10. Comparison of quiescent and disturbed-blackout conditions of theoretical reflection coefficient, T, 
(j=1), for various fixed heights, h, of sharply bounded model ionosphere for LF transmissions from Adak to 
Kodiak, Alaska. 
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All computations were made with the assumption of a dipole current moment, Jj/=1. 
The total radiated power, Pr=1.3 kw, for the arctic data considered was found to correspond 
to a dipole current moment, Jp/=2050 amp-m; therefore, the calculated field amplitude, [E(w,d)], 
or [E;(w,d)], has been increased by this factor for comparison with the data. 
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In addition to the total field, H(w,d), and the modal field, Zj(w,d), figures 6 to 9, the 
reflection coefficients, 7, (horizontal and vertical polarization; normal and abnormal com- 
ponents), figures 10 to 11, and the transmission data of the ordinary and extraordinary upgoing 
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Attenuation rate and phase lag rate at various emplacement altitudes, h, for upgoing “‘ordinary’”’ and 


“extraordinary”? waves excited in the model ionosphere by Adak to Kodiak, Alaska transmissions during 


quiescent and disturbed-blackout conditions. 
First time-mode: H»,=0.5035, ¢2=51.08°, I=67.18°. 
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and disturbed-blackout conditions. 


First time-mode: H»=0.5187, ¢2=12.27°, [=68.68°. 
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waves excited in the ionosphere by the reflection process, figures 12, 13, were evaluated. 


The 


unabsorbed modal field, £j(w,d)/C;, and the effective reflection coefficients, C;, are separately 


illustrated, figures 14 and 15, respectively. 


The emplacement height, h, of the model ionosphere was deduced by comparison of the 


computed field strength, E(w,d), with the measured value. 


The corresponding four reflection 


coefficients, 7.2, Tem, Timm, Tme, (amplitude and phase) were then evaluated from the theory. 
The results of this process which were readily scaled from graphs such as illustrated, figures, 10, 
11, are presented in table 4, and the reflection coefficients deduced for the first time-mode only 


are presented. 


The higher order time-mode, j=2, 3, reflection coefficients are obviously quite numerous 


since a separate set exists for each ionosphere reflecting region. 


These reflection coefficients 


are determined by the analysis by reference to the emplacement altitude, h, table 5, deduced 
from the field, E(w,d). 
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Effective reflection coefficient, c;, for Adak to Nome and Kodiak, Alaska, transmissions under 








os 


























| ] 

Condition | Mode N h | |[E(@,d)| | arg Ee 
El/em3 km v/m radians 

See ee ero! Total __. 10 | 55 7. 59 (10-*) | 5.5 
| IRE Poe | Total_.. 1000 60} 6.60 (10-*) | 2. 84 
SE j=1.--- 10 | 55 | 1.18 (10-19) | 2.10 
RE ae -| Total -.. 10 | 65 | 9.77 (10-) | 5. 66 
Di iinduesacones ces Seen 10 | 65 | 2. 42 (10-9) | 2. 33 





It is interesting to note emplacement altitudes, h, less than 65 km under quiescent conditions. 
Indeed, the highest altitude for these daytime values was 69 km (winter 1954). A good nominal 
daytime emplacement altitude is 67.5 km. The Jow value (|7|<0.1) reflection coefficients are 
also quite interesting. Note the highest value, |7|~0.27. Note also, the median value at 
Nome (1954 equinox), |7..|=0.08. A low value |7|<0.1, and “black’’, | 7|~0 appearance of the 
ionosphere indicates a low level absorbing region in the lower ionosphere at grazing incidence. 
A region of low electron density, N=10-100 El/cm’, is indicated by the model employed, 
figure 4. The transmission into such a region at grazing incidence, ¢,~80°, figures 11, 12, 
indicates almost complete absorption of both the “ordinary” and the “extraordinary” rays if a 
“virtual height,” A+Ah, is assumed of the order of 70 km. The transmission curves, figures 
11, 12, demonstrate that the “extraordinary” ray is very highly attenuated such that it is 
necessary to consider only the “ordinary” ray. Assuming 1 db/km, for example, figures 11, 12, 
and a virtual height, h+ Ah, of 70 km, noting that appreciable attenuation rate begins at 65 km, 
the total attenuation can be estimated; total decibels ~2Ah/cos ¢;~57 db. This, of course, 
neglects a reflection absorption at the 70 km virtual height plus the transmission absorption 
back into the region below the ionosphere. The values, Ah, can be decreased and a correspond- 
ing decrease in the attenuation noted. But a repetition of the process leads to the interesting 
conclusion that most of the ray reflection must occur in a fraction of a kilometer, and, hence, 
the sharply bounded model can be reasonably applied under these conditions to predict the LF 
field, H(w,d). A similar examination of the second and third time-modes leads to an almost 
similar reflecting region since the angle of incidence, ¢,, figure 2, does not change rapidly until 
the distance, d/j, is sufficiently short. It is, nevertheless, conceivable that the amplitude of the 
skywave reflections, 7=4, 5, 6, . . . at shorter distances could become quite large as a result 
of reflection from a high level “virtual height,” h+Ah=70 to 85 km. Under these conditions, 
however, the sharply bounded model would not be valid and it would be necessary to develop 
theoretical techniques for a diffusely bounded model ionosphere. 

Information on the lower ionosphere N(h) profile during nighttime conditions is quite 
sparse and hence the theoretician can only speculate on its nature. However, under certain 
reasonable assumptions, a possible explanation of the increase in field strength experienced at 
night may be deduced from the daytime model ionosphere employed. 

Despite the seeming “evaporation” of the lower regions of the ionosphere at night with a 
corresponding abrupt increase of the electron density to a large finite value, it seems reasonable 
to assume that this abrupt increase does not occur within \/10 or 0.22 km at 135.6 ke/s. This 
makes it possible to utilize the present sharply bounded model ionosphere, if ii is emplaced at a 
higher level. Recent work [Ellyett and Watts, 1959] indicates a height of 80 to 85 km is not 
unreasonable. 

Since the change of the angle of incidence, ¢;, figure 2, with ionosphere height, h, is not 
great at this range, we may neglect its effect on the field strength, figures 6 to 9, and the effective 
reflection coefficients, figure 15. The only factor changed, therefore, is the unabsorbed field, 
figure 14. It therefore seems reasonable to expect the first time mode, |/;(w,d)|, to increase 
on the order of at least twice, whereas, reference to the arctic data and figure 16, “Battle Lake 
data”’ (to be discussed below), indicates an increase of perhaps 4 to 10 times. ‘This could easily 
be explained by a slight change of slope in the N(h) profile employed as evidenced by the steep 
slope of the C; curves (j=1), figure 15, in the region involved. On this basis, it is reasonable 
to expect a nighttime reflection coefficient, |7|~0.3, 0.4. Indeed, at temperate latitude, these 
values are confirmed [Ross, 1959]. 
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It is not surprising that a large phase change could be expected during the sunset-sunrise 
transition periods, figures 6 and 8. Indeed, such phase changes were described at 100 ke/s 
by Doherty [1957], while operating the Loran-C on the first time-mode of the skywave. 

The behavior of the LF propagation during disturbed and blackout conditions which are 
especially frequent in the arctic-auroral region is not quite so clear. Figures 7 and 9 illustrate 
disturbed-blackout conditions based on the geophysical measurements of electron density 
profiles, N(h), of the lower ionosphere, figure 4. The corresponding reflection coefficients, 
figures 10, 11, and the transmission of upgoing ordinary and extraordinary waves excited by 
the reflection process are also illustrated, figures 12, 13. The disturbance of the earth’s 
magnetic field during extreme blackout conditions is negligible since the total field, H,,~0.5 
gauss and the fluctuations are of the order of 0.001 gauss [Chapman, 1940]. The most 
significant change in the propagation is caused by the drastic change in the electron density 
profile, N(A), figure 4. Unfortunately, the shape of the lower levels, <60 km, of the “black- 
out”? N(h) profile is unknown. Since the analysis under quiescent conditions indicates low- 
electron density reflections, the curve, figure 4, was extrapolated, and the corresponding field, 
FE(w,d), calculated for Adak to Nome transmission. The results are tabulated in table 5. 

Note the 7=1 mode for a region of N=10 El/cm* has undergone a 26 db decrease and 
the total field is therefore almost entirely groundwave. Appreciable as the decrease is, it is 
still well within the range of observations of arctic data. Notice, figure 16 “Battle Lake 
data’’ (to be discussed below), that the range of observations in temperate latitudes is approx- 
imately 18 db. The additional fluctuation may then, reasonably, be attributed to the change 
in electron density of the extreme lower ionosphere. If this change were to be so abrupt as 
to reach N= 1000 E}/cm in, say, 4/10 km, the field strength would show a considerable increase, 
instead of decrease. It is reasonable to expect that a more realistic value of somewhat less 
than N=1000 El/em? would also produce an increase in the received field strength. The 
shape of the N(A) profile, thus, becomes quite important. 
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pulse amplitude measurements with theoretical first 
time-mode in the diffraction zone beyond the geo- 
metrical-optical horizon. 
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The calculation of the amplitude and phase of the total field, E(w,d) as a function of 
distance is illustrated in figures 16 and 17. The corresponding reflection coefficient 7, 
(vertical polarization, normal component) is illustrated in figure 18. 

Amplitude measurements *® during the day which illustrate the diffraction of the first 
time-mode (j=1) near and beyond the geometric-optical horizon are illustrated in figure 16. 
The agreement with the theoretical curve j7=1 appears to be quite close except at the range 
of 1,800 miles. It is interesting to note the apparent decrease in the range of amplitude 
observations in the daytime, indicating greater amplitude stability with range in the 
diffraction zone. 
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4. Conclusion 


The analysis of these data demonstrates the application of propagation theory together 
with a suitable model for the ionosphere to the prediction of LF phenomena. The precision 
of the model for the case of transmission via the ionosphere at or close to grazing incidence is 
considered to be adequate upon considering the state of knowledge of the many parameters 
which influence such propagation. Nevertheless, the absence of an abundance of specific 
information does not deter the author from certain general conclusions which can be made as 
a result of this research. 

At oblique incidence on the ionosphere, a condition which is certainly met under most 
practical situations, the reflection of LF waves during the daytime occurs in a region of very 
low electron density and high absorption, with corresponding low values of reflection coefficients, 
|T|<0.1. The dominant skywave time-mode is, in general, the first; however, the second, 
third and higher order time-modes may offer appreciable contribution to the total field as the 
angle of incidence decreases as a result of reflection from higher levels and shorter transit 
distance in the ionized medium. 

Periods of ionospheric disturbance (blackout) can produce a decrease in the field as a result 
of a reduction in the amplitude of the ionospheric time-modes. However, it is equally possible 
for the skywave mode to be incremented. This phenomenon is dependent upon the slope of 


6 From the author’s persona! notes on the 1953-1954 testing of the Loran-C (Cytac) System. 
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the electron density-altitude profile, N(h) h. However, if ordinary communication is an 
engineering objective, LF seem to be quite reliable since the characteristic high absorption of 
HF waves does not appear to exist. 


Since the writing of this paper the extensive theoretical work of Barron [1961] and the 
extensive experimental and engineering work of Belrose et al., [1959] have come to the attention 
of the author. It is interesting to note quite good agreement, figure 18, between the reflection 
coefficients deduced from the sharply bounded model presented in this paper and those pre- 
sented as a result of Barron’s full wave computations. Barron’s data does, however, indicate, 
as might be expected, higher attenuation at the higher frequencies (>150 ke/s). 

It should be noted in conclusion that the adequacy of the sharply bounded model is subject 
to severe limitations, especially at the higher frequencies and the results of this paper suggest 
also the introduction into this geometric optical analysis of more rigorous or flexible models 
for the reflection coefficient which would treat a more diffuse ionosphere lower boundary. 
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Mean electron density variations of the quiet ionosphere 3— 
May 1959, J. W. Wright, L. R. Wescott, and D. J. Brown, 
NBS Tech. Note 40-3 (PB151399-3) (1960) $1.59. 

The CRPL has initiated a program for large-scale computa- 
tion of electron density profiles from ionospheric vertical 
soundings. Scaling is performed at field stations permitting 
computation of hourly profiles at the Central Laboratory. 
These profiles are combined to form hourly mean quiet pro- 
files for each station and month. The results of this program 
for the month of May are illustrated graphically. This report 
is the third of a series illustrating the electron density varia- 
tions in the mean quiet ionosphere between latitudes 15°N 
and 50°N along the 75°W meridian. 


On the nature of equatorial spread -F,R. Cohen and K. L. 
ow J. Geophys. Research 66, No. 4, 1081-1106 (April 
961). 

Ionospheric propagation via scattering from the F region 
was sought at 50 Me/s over a transequatorial path (with 
midpoint near Huancayo, Peru) employing a 2580-km trans- 
mitter-receiver separation. Propagation via F scatter was 
present over this path about 10 per cent of the time, though 
only at night. A condition closely related to the occurrence 
of F scatter was the presence of equatorial spread-F configura- 
tions on the Huancayo ionograms. The height of the propa- 
gation medium supporting this F scatter was usually identifiable 
with the lowest height of the associated equatorial spread F 
on the ionograms. On the basis of an interpretation of ex- 
perimental results, equatorial spread F is shown to arise from 
scattering by relatively thin sheets of irregularities in electron 
density which occur at the bottom of the F layer or as much 
as 100 km lower. The thickness of these patches is estimated 
at the order of 50 km. The scattering irregularities comprising 
these patches are found to be elongated along the earth’s 
magnetic lines of force. Electron irregularities of scale 10 
meters or smaller measured in at least one dimension trans- 
verse to the magnetic field lines, and 1,000 meters or longer 
measured along the magnetic field lines, are shown to exist 
as high as 450 km in the nocturnal equatorial F region during 
equatorial spread-F conditions. Equatorial spread-F echoes 
observed at Huancayo are demonstrated to be arriving at all 
elevation angles in the magnetic east-west plane. The geo- 
graphical extension of a given scattering sheet in the magnetic 
east-west direction was at times as great as 1000 km. An 
observational procedure is suggested for distinguishing two 
fundamental varieties of spread-F echoes appearing on the 
equatorial ionograms. A necessary condition, that the con- 
tours of mean electron density be parallel to the magnetic 
lines of force, is proposed as a controlling factor for the oc- 
currence of spread F in the equatorial ionosphere. 


On the analysis of polarization rotation recordings of satellite 
radio signals, R. 8. Lawrence and C. G. Little, pp. 391-399 
of Some ionospheric results obtained during the International 
Geophysical Year; Proc. symp. organized by the URSI/AGI 
committee, Brussels, 1959, edited by W. J. G. Beynon (Elsevier 
Publ. Co., Amsterdam, 1960). 

At frequencies as low as 20 mc/s the rotation of the plane of 
polarization of satellite signals cannot be assumed to be pro- 
portional to f/ N B cos @ dl along the line of sight. A method 
of analysis is presented which eliminates this assumption 
and permits accurate estimates of total electron content below 
the satellite height to be obtained from 20 me observations. 
Some of the results obtained using the method to determine 
the sub-satellite electron content, throughout a satellite pass, 
are presented. These analyses permit the study of ionospheric 
tilts, and have revealed the presence of large-scale irregulari- 
ties in the ionospheric electron content. The deviation in 
sub-satellite content of these irregularities is of the order 





2% of the running mean value: the lateral extent of the irregu- 
larities is of the order of a few hundred kilometers. The ef- 
fects of such irregularities upon ionospheric refraction and 
satellite Doppler curves are briefly discussed, 


Prolonged space-wave fadeouts in tropospheric propagation, 
A. P. Barsis and M. E. Johnson, NBS Tech. Note 88& 
(PB16158$S) (1961) $2.00. 

This paper contains the results of studies performed during 
the last several years on the short-term variability of tropo- 
spheric signals received over within-the-horizon paths. Signal 
variations of this type have been termed “prolonged space- 
wave fadeouts,’”’ as they are mainly characterized by reduc- 
tions in signal level to many decibels below presumably 
constant values determined from geometrical optics methods. 
The data described here were obtained from measurements 
over propagation paths in the Pacific Coast region of Southern 
California, and the continental region of Eastern Colorado. 
Fadeouts are analyzed as a function of carrier frequency, 
path characteristics, and meteorological parameters. The 
study also includes an evaluation of fadeouts observed over 
a path using a mountain peak as a diffracting knife-edge like 
obstacle between transmitter and receiver. 


The relationship of low-height ionosonde echoes to auroral- 
zone absorption and VHF D scatter, J. K. Olesen and J. W. 
Wright, J. Geophys. Research 66, No. 4, 1127-1134 (April 
1961). 

Weak diffuse HF reflections at heights between 75 and 95 
km in the auroral zone are designated H,d. Diurnally, E,d 
occurs most frequently at noon at all seasons, and shows a 
pronounced minimum of occurrence around 2000 hours local 
time. There is little evidence for a seasonal variation from 
the year’s data presented here. It is suggested that E,d is 
related to the anomalies of auroral-zone absorption and that 
it is the layer responsible for VHF forward scatter. 


The Ephi system for VLF direction finding, G. Hefley, 
R. F. Linfield, and T. L. Davis, J. Research NBS 65C (Fng. 
and Instr.), No. 1, 43-49 (Jan.-Mar. 1961). 

A new system of VLF direction finding has been developed 
and tested by the National Bureau of Standards, Boulder, 
Colorado. The system has been named “Ephi’’ (E-¢) 
because the bearing of the transient signal is determined from 
the relative phase (¢) of the vertical electrical field (E) 
received at spaced antennas. The advantage of this scheme 
compared to conventional crossed-loop techniques is that it 
minimizes siting and polarization errors. A minimum of 
three antennas must be used to resolve directional ambiguity, 
and the preferred antenna location is at the vertices of an 
equilateral triangle with baselines equal to 4 to 40 wavelength 
(at 10 ke/s). Appropriate phase detectors, delay lines and 
coincidence circuits are used to obtain a directional code in 
preset sectors. Within practical instrumentation limits any 
number of sectors of variable widths can be operated simul- 
taneously and each can be rotated in azimuth. 


Integrated starlight over the sky, F. E. Roach and L. R. 
Megill, Astrophysical J. 133, No. 1, 228-242 (Jan. 1961). 
Calculations are made of the total integrated starlight over 
the entire sky based on the star counts in Groningen Pub., 
No. 43. The results are given in both the photographic and 
the visual magnitude scales in tabular and graphical form. 


Resonators for millimeter and submillimeter wavelengths, 
W. Culshaw, JRE Trans. Microwave Theory and Tech. MTT-9, 
No. 2, 185-144 (Mar. 1961). 

Further considerations on the mm-wave Fabry-Perot inter- 
ferometer are presented. Computed @ values for parallel 
metal plate resonators indicate that at spacings around 2.5 
cm, values ranging from 60,000 at 3 mm, to 300,000 at 0.1 mm 
wavelengths are possible. The plates must, however, be 
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quite flat. These results are important for many investiga- 
tions, and in particular for mm and sub-mm wave maser 
research. For the aperture per wavelength ratios possible 
here, diffraction effects should be small. Consideration is 
given to using curved reflectors or focused radiation in appli- 
cations where the fields must be concentrated. For this 
purpose, re-entrant conical spherical resonators are treated 
in detail, as regards operation in the TEM mode at high 
orders of interference. Expressions for the Q and shunt 
impedance are given, and high values are possible at mm and 
sub-mm wavelengths. Quasi-optical methods of coupling 
into and out of such a resonator are proposed, and the higher 
modes possible in such a resonator are considered. Results 
indicate that it could have application to the mm-wave 
generation problem, and that it represents a good resonant 
cavity for solid-state research at mm and sub-mm wave- 
lengths, and for maser applications in particular. 


The electromagnetic fields of a dipole in the presence of a 
thin plasma sheet, J. R. Wait, Appl. Sci. Research, Section B, 
8, 397-417 (1960). 

The problem of electric and magnetic dipoles located near a 
thin planar slab or sheet of ionized material is considered. 
A constant and uniform magnetic field is impressed on the 
slab. Under the assumption that the thickness of the slab 
is very small, expressions for the resultant fields are obtained. 
As a result of the anisotropy of the sheet it is indicated that 
the fields are elliptically polarized in general. On carryiag 
out a saddle-point evaluation of the integrals in the formal 
solution it is shown that the far fields may be split into 
“radiation” and ‘‘surface wave’’ components. The depend- 
ence of the radiation pattern and the surface wave character- 
istics on electron density, collision frequency and the impressed 
magnetic field is illustrated. 


The Nation’s electronic standards program: Where do we 
now stand? Harvey W. Lance, JRF Trans. Instrumentation 
I-9, No. 2, 94-100 (Sept. 1960). 

Recently increased standards activities throughout the nation 
have resulted in a large increase in the calibration services 
requested from the National Bureau of Standards and in a 
number of related problems. Several of these are discussed, 
including desirable criteria for the staffing and instrumenta- 
tion of electronic standards laboratories, the choice of inter- 
laboratory standards to be submitted to NBS, and the use 
of these standards before and after calibration. An attempt 
is made.to clarify the concept of ‘‘traceability of calibra- 
tions to NBS,” and conditions under which a calibration 
may be called traceable are discussed. Suggestions are 
made regarding the use of available calibration services and 
the action to be taken when needed calibration services are 
not available. The need for an association of standards 
laboratories is pointed out, and many areas are noted in 
which such an association could make valuable contributions. 


Propagation of error in a chain of standards, A. G. McNish 
and J. M. Cameron, JRE Trans. Instrumentation I-9, No. 2, 
101-104 (Sept. 1960). 

There is a widespread belief that in a chain of standards 
traceable back to a central laboratory, each standard must be 
calibrated in terms of a standard 10 times as accurate. 
Conversely, this implies a 10-fold degradation in accuracy in 
each successive calibration. If this were true, accuracy 
would be intolerably degraded at the working level. A more 
realistic view is that accuracy in a chain of standards is as 
good as its weakest link. This can be demonstrated by 
experimental data on the errors obtained in calibrations and 
can be explained by examining the factors which affect 
accuracy in calibration. 


An analysis of the accumulated error in a hierarchy of cali- 
brations, E. L. Crow, JRE Trans. Instrumentation I-9, No. 2, 
105-114 (Sept. 1960). 

Calibrations of many types are performed in a hierarchy of 
calibration laboratories fanning out from a national standard. 
Often the statement is made that the accuracy of each echelon 
of the hierarchy should be 10 times the accuracy of the immedi- 
ately following echelon. The validity of such statements is 


examined by deriving formulas for the total error accumulated 
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over the entire sequence when systematic and random errors 
may occur in each echelon, and by determining how a given 
total error may be achieved at minimum total cost under 
reasonable assumptions for the form of the cost-error functions. 


Atomic beam frequency standards, R. C. Mockler, R. E. 
Beehler, and C. S. Snider, JRE Trans. Instrumentation I-9, 
No. 2, 120-132 (Sept. 1960). 

A general qualitative description of atomic beam devices is 
given, with particular emphasis upon cesium atomic frequency 
standards. 

The various uncertainties in frequency measurements are dis- 
cussed including the effects of the uniform C field, the effect of 
a phase difference between the two oscillating fields that 
excite the atomic transition, and the effect of exciting an 
atomic transition with signals that are not monochromatic. 
The two cesium beam devices constructed at the National 
Bureau of Standards are described in some detail. The 
various tests are described for determining the uncertainties 
in the frequency measurements. The results of these tests 
are given. 

The accuracy of these rather dissimilar machines is considered 
to be +1.5X10-'!, and the precision of measurement is 
2x10-'*. The latter figure represents the standard deviation 
of the mean for measurement times of a few hours, or, about 
100 separate measurements. 


The power spectrum and its importance in precise frequency 
measurements, J. A. Barnes and R. C. Mockler, JRE Trans. 
Instrumentation I-9, No. 2, 149-155 (Sept. 1960). 

The power spectral density functions of a frequency multiplier 
chain, driven by several different crystal oscillators, were 
obtained by comparing the output with a second chain which 
was stabilized with an ammonia maser. The frequency of 
the maser stabilized chain was demonstrated to be relatively 
fixed; the power spectrum of the other chain was determined 
by two different methods. The results are compared. 
Possible errors and uncertainties introduced by the methods 
are discussed. An analysis is made that relates the instanta- 
neous frequency fluctuations of a signal with the power spectral 
density function. 

Analysis predicts that when frequency modulation occurs in 
the first stages of frequency multiplication or in the primary 
frequency oscillator, the output power spectrum is, in general, 
not symmetrical. Furthermore, the sidebands are increased 
in intensity by the multiplication process. This is, in fact, 
observed to be the case. It is shown that a frequency counter 
will measure the frequency of the center of gravity of the 
power spectrum. 

If signals having a complex power spectrum are used in precise 
frequency measurements, errors may result. 


Standards and measurements of microwave surface imped- 
ance, skin depth, conductivity and Q, H. E. Bussey, JRE 
Trans. Instrumentation I-9, No. 2, 171-175 (Sept. 1960). 

The metal walls of a cavity resonator become standards of 
skin depth when the internal Q is determined, if the walls are 
uniform and if contact loss is eliminated as in the circular 
TEo: mode. An end plate standardized in this way may be 
replaced by another flat piece of material and its unknown skin 
depth easily determined. A copper cavity in use has 94 
per cent of theoretical Q. Deoxidization of the surface raises 
the value to 97 per cent. Surface roughness and mode con- 
version may account for the remaining Q departure. Accuracy 
of Q measurement (+3 per cent) will be discussed and pre- 
liminary results on several metals will be given. 


A radio-frequency permittimeter, R. C. Powell and A. L. 
Rasmussen, [RE Trans. Instrumentation I-9, No. 2, 179-184 
(Sept. 1960). 

A coaxial RF impedance transformer in which the secondary is 
a single turn of the material to be measured is used with two- 
terminal impedance bridges to determine the complex per- 
mittivity or complex conductivity of low-impedance materials. 
No electrodes are needed, and many conductors, semiconduc- 
tors, electrolytes and high-permittivity materials can be 
evaluated to about 1 per cent, since errors due to electrode 
impedance and interaction as well as first-order series induct- 














ance are eliminated. The design, calibration, range, and 
accuracy are given along with measured values of such mate- 
rials as ferrites and strong electrolytic solutions, showing 
complex conductivities differing considerably from those 
previously observed by other methods. 


A transfer instrument for the intercomparison of microwave 
power meters, G. F. Engen, JRE Trans. Instrumentation 
I-9, 202-208 (Sept. 1960). 

In the intercomparison of two microwave power standards, 
or in the calibration of a microwave power meter by means 
of a second or “‘standard’’ power meter, the measurement of 
microwave impedance has played a major role. 

Through an extension of the reflectometer concept, it is now 
possible to devise a four-arm junction which, when properly 
adjusted, makes possible the intercomparison or calibration 
of such power meters with little or no regard for their im- 
pedance characteristics. In addition, the method is substan- 
tially independent of the impedance discontinuity which may 
be present at the input flange or connector. This latter result 
is of particular value in coaxial systems. 


Absolute measurement of temperatures of microwave noise 
sources, A. J. Estin, C. L. Trembath, J. S. Wells, and W. C. 
a” IRE Trans. Instrumentation 1-9, 209-213 (Sept. 
1960). 

In recent years, the requirement for noise-figure measure- 
ments of low-noise devices has led to a need for more precise 
knowledge of temperatures of noise sources. A system has 
been developed which is capable of measuring the absolute 
noise temperature of a gas discharge tube in WR-90 wave- 
guide to within 0.01 db. The comparison circuit is a modified 
Dicke radiometer. The standard noise source is a black 
body radiator at approximately 1000°C. Preliminary meas- 
urements of noise temperatures of a set of nine laboratory 
equivalents of a commonly used argon discharge tube indicate 
an excess noise temperature of 15.90 db at 200 ma discharge 
current with a variation between tubes of 0.03 db. 


A modulated subcarrier technique of measuring microwave 
phase shifts, G. E. Schafer, JRE Trans. Instrumentation I-9, 
217-219 (Sept. 1960). 

This paper describes a technique for comparing phases of two 
microwave signals which employs amplitude modulation in 
one channel of a two-channel system. This modulated sub- 
earrier technique produces a null response in a two-channel 
system for all ratios of the amplitudes of the waves traversing 
the separate channels. Therefore, high precision is main- 
tained for all ratios of amplitudes. A practical application 
of this technique using readily available components achieved 
a precision of phase measurements of a tenth of a degree for 
a change of attenuation of 50 db of component under test. 


Measurement of reflections and losses of waveguide joints 
and connectors using microwave reflectometer techniques. 
R. W. Beatty, G. F. Engen, and W. J. Anson, JRE Trans 
Instrumentation I-9, 219-226 (Sept. 1960). 

The reflection and loss of a waveguide joint are quantified 
by its VSWR and efficiency, respectively. These are con- 
veniently and accurately determined by techniques which 
employ a reflectometer with auxiliary tuners and an adjustable 
sliding termination. Depending upon the stability and gain 
of the associated apparatus, measurements of VSWR’s of 1.001 
and lower, and of efficiencies of 99.99 per cent and higher 
may be obtained. The attenuations of short lengths of wave- 
guide may be determined at the same time that the efficiency 
measurement is made. 


A precision RF attenuation calibration system, C. M. Allred 
and C. C. Cook, IRE Trans. Instrumentation I-9, 268-274 
(Sept. 1960). 

A new precision attenuation calibrating system with greatly 
increased sensitivity, stability, and measurement range has 
been completed at the National Bureau of Standards, Boulder 
Laboratories. The increased stability and sensitivity are 
achieved by the use of a highly accurate piston attenuator 
and precision phase shifter combined into a null system. 
The extended attenuation measurement range has been ob- 





tained by using a new mode launching system which is excited 
by essentially a constant current source of very high magni- 
tude. The padding necessary for correct impedance matching 
is kept at a minimum by use of a special noninteractive com- 
bining network. 


A precision RF power transfer standard, P. A. Hudson, 
IRE Trans. Instrumentation 1-9, 280-283 (Sept. 1960). 

A fast-reading, precision RF power transfer standard is 
described which is usable at fixed frequencies in the range 
10 Me to 1000 Me. The dynamic range is 1 watt to 1000 
watts with a resolution of 0.1 per cent to 0.01 per cent. 
Stability is 0.1 per cent per year. The standard consists of 
three directional couplers to cover the frequency range and 
a series of vacuum thermoelement detectors matched to 50 
ohms at each frequency interest. The thermoelements are 
connected to the secondary line of the appropriate coupler 
and provide a de output voltage which is calibrated in terms 
of RF power in the coupler main line. 

A series of small, portable wattmeters employing vacuum 
thermoelements are also described. 
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